Therapeutic drug monitoring for busulfan is currently performed by multiple plasma sampling. Saliva is considered a noninvasive therapeutic drug monitoring matrix. This study aimed to investigate intravenous busulfan pharmacokinetics (PK) in plasma and saliva, and establish a limited sampling strategy (LSS) for predicting the area under the concentration–time curve from time zero to infinity in plasma (AUC0‐∞,p) by using saliva samples. Therefore, the PK of busulfan was studied in 37 Chinese patients. Pearson correlation analysis was used to evaluate the correlation between the AUC of busulfan in plasma and saliva. LSS models were established by the multiple linear regression analysis. The prediction error, the mean prediction error, and the root mean square error were used to evaluate the predictive accuracy. The agreement between the predicted and observed AUC0‐∞ in saliva was investigated by the intraclass correlation coefficient and Bland–Altman analysis. The accuracy and robustness of the models were evaluated by using the bootstrap procedure. The result of PK analysis 62.2% of patients (23/37) was within the target range of AUC0‐∞,p. A good correlation between saliva and plasma busulfan AUC0‐∞ was observed (r = 0.63, p < .01). The bias and precision of the models 7 and 13 were less than 15%. The intraclass correlation coefficient exceeded 0.9, and the limits of agreement were within ±15%. The 2‐point LSS model in saliva is a convenient and desirable approach to predict the AUC0‐∞ of 4 times daily intravenous busulfan in plasma, which can be used to design personalized dosing for busulfan.