We investigated the stochastic response of a person sitting in a driving vehicle to quantify the impact of an uncertain parameter important in controlling defect reduction in terms of ride comfort. Using CarSim software and MATLAB/Simulink, we developed a fully coupled model that simulates a driving vehicle combined with an analytical nonlinear human model. Ride comfort was evaluated as a ride index considering the frequency weights defined in BS 6841. Additionally, to investigate the uncertainty of the ride index, a framework for calculating the ride index was proposed using the generalized polynomial (gPC) method. Further, sensitivity analysis of the ride index was performed for each uncertainty parameter, such as stiffness and damping. The results obtained through the gPC method were in good agreement with those obtained via Monte Carlo simulation (MCS) and were excellent in terms of computation time without a loss of numerical accuracy. Through in-depth investigation, we found that the stochastic distribution of the ride index varies differently for each uncertain parameter in the human model. By comparing linear and nonlinear human models, we also found that the nonlinearity of the human model is an important concern in the stochastic estimation of ride comfort.