Voriconazole is the first-line therapy for invasive aspergillosis (IA). To determine the minimum inhibitory concentration of Aspergillus, a voriconazole pharmacokinetic-pharmacodynamic (PK-PD) model linked to galactomannan response was developed and evaluated, and its clinical correlation for IA treatment was elucidated. Adult patients with probable or definite IA and at least one serum voriconazole measurement were included. A two-compartment voriconazole PK model was linked to a previously described PD model of galactomannan response. PK and PD parameters were estimated using a nonparametric adaptive grid technique. The relationship between the ratio of voriconazole exposure that induced half-maximum galactomannan response (EC50) and the observed terminal galactomannan concentration was evaluated. The factors associated with the PK-PD parameters and mortality were also determined. Between January 2013 and December 2022, 41 patients were prescribed voriconazole for IA. The 30-day mortality rate was 17%. A high correlation was found for the observed-predicted Bayesian posterior estimates of voriconazole and galactomannan levels. Moreover, a nonlinear relationship was identified between AUC:EC50 and terminal galactomannan. The factors associated with higher AUC:EC50 were intravenous administration and intubation. In the survival analysis, higher EC50 tended to be associated with mortality, higher AUC was significantly associated with increased mortality, and higher AUC:EC50 tended to be associated with higher mortality. After adjusting for the intravenous route, higher AUC and AUC:EC50 were not associated with mortality. Individual EC50 estimation can provide insights into in vivo host and organism responses. Elevated EC50 showed comparable and unfavorable trends to higher minimum inhibitory concentration. Thus, determining EC50 might help guide individualized target serum voriconazole levels.