Objectives To characterize the pharmacokinetics (PK) of polymyxin B in Chinese critically ill patients. The factors significantly affecting PK parameters are identified, and a limited sampling strategy for therapeutic drug monitoring of polymyxin B is explored. Methods Thirty patients (212 samples) were included in a population PK analysis. A limited sampling strategy was developed using Bayesian estimation, multiple linear regression and modified integral equations. Non-linear mixed-effects models were developed using Phoenix NLME software. Results A two-compartment population PK model was used to describe polymyxin B PK. Population estimates of the volumes of central compartment distribution (V) and peripheral compartment distribution (V2), central compartment clearance (CL) and intercompartmental clearance (Q) were 7.857 L, 12.668 L, 1.672 L/h and 7.009 L/h. Continuous renal replacement therapy (CRRT) significantly affected CL, and body weight significantly affected CL and Q. The AUC0-12h of polymyxin B in patients with CRRT was significantly lower than in patients without CRRT. CL and Q increased with increasing body weight. A limited sampling strategy was suggested using a two-sample scheme with plasma at 0.5h and 8h after the end of infusion (C0.5 and C8) for therapeutic drug monitoring in the clinic. Conclusions A dosing regimen should be based on body weight and the application of CRRT. A two-sample strategy for therapeutic drug monitoring could facilitate individualized treatment with polymyxin B in critically ill patients.