This study attempts to investigate the influence of the dynamic soil property uncertainties on the ground response analysis via a case study. For this purpose, nonlinear time-domain ground response analysis and uncertainties in soil parameters are coupled by a MATLAB code. To take full advantage of the real data, two investigation boreholes are drilled in the site. The analysis was performed deterministically and then extended to the stochastic context to consider the variability of the plastic index, shear wave velocity, and unit weight of the soil. Furthermore, the capability of the three different methods for predicting the stochastic fundamental period, including modal analysis, approximate method, and nonlinear method, is investigated. The comparison reveals that the mean value of the approximate method provides closer predictions regarding the fundamental periods obtained through the nonlinear method. In the stochastic analysis, the maximum Coefficient of Variation (COV) of the peak ground motion parameters, fundamental period, response spectrum, and amplification factor is calculated. The results demonstrate that the soil parameters heterogeneity has a significant effect on the variation of the surface Peak Ground Displacement (PGD). Among the other stochastic responses, the fundamental period has received the least effect from soil parameters' uncertainty.