Invasive fungal disease (IFD) is a significant cause of morbidity and mortality in hospitalized patients. To maximize the efficacy of voriconazole treatment, the study established the relationship between voriconazole pharmacokinetic/pharmacodynamic (PK/PD) and probability of response and optimized voriconazole dosage regimen in patients with IFD based on Monte Carlo simulation. Forty-four patients proven with IFD were involved in this study. Among them, the overall cure rate was 75% (33/44) and there was a significant difference between Cmin /MIC values in patients with lack of response (n = 11) and those with successful response (n = 33) (mean value: 1.91 vs. 11.33; P < 0.05). Logistic regression model showed a high correlation between voriconazole Cmin /MIC ratio and clinical response (P = 0.044, OR = 1.349). According to Monte Carlo simulation results under different voriconazole dosing regimens, we could draw a conclusion that 200 mg voriconazole administered intravenously or orally twice daily for Candida infections and 300 mg administered orally or with 200 mg administered intravenously twice daily for Aspergillus infections were rational, which could achieve a value of the cumulative fraction of response >90%. This study built the relationship between voriconazole PK/PD and clinical response and obtained the reasonable empirical dosage regimen, which can be used to customize individual dosage regimen and improve the efficacy of voriconazole treatment.