Etanercept pharmacokinetics in patients with rheumatoid arthritis (RA), ankylosing spondylitis (AS), and psoriasis were assessed separately with distinct models using population pharmacokinetics methods of limited precision. The different model structures and associated significant covariates identified by these earlier methods made it difficult to compare etanercept pharmacokinetics among disease groups. This integrated analysis aimed to establish a framework to evaluate previously established population pharmacokinetic models of etanercept, and to identify consistent and important demographic and disease factors that affected etanercept pharmacokinetics in a diverse population of healthy subjects and patients with RA and AS. In this integrated analysis, cumulative rich and sparse etanercept concentration data from 53 healthy volunteers, 212 patients with RA, and 346 patients with AS were examined and compared using nonlinear mixed effect methodology implemented the in NONMEM VI software package. A more precise estimation method (FOCEi) was employed and compared with the first-order method in population pharmacokinetics model building and evaluation. The integrated analysis found that an optimal population pharmacokinetics model with a 2-compartment structure adequately characterized etanercept pharmacokinetics in all subject groups. Health status or disease type did not significantly affect etanercept pharmacokinetics. In adult patients with RA and AS, age and body weight do not significantly affect etanercept pharmacokinetics.