It is recognized that there is a dearth of studies that provide a comprehensive understanding of vehicle-occupant system dynamics for various road conditions, sitting occupancies and vehicle velocities. In the current work, an in-house-developed 50 degree-of-freedom (DOF) multi-occupant vehicle model is employed to obtain the vehicle and occupant biodynamic responses for various cases of vehicle velocities and road roughness. The model is solved using MATLAB scripts and library functions. Random road profiles of Classes A, B, C and D are generated based on PSDs (Power Spectral Densities) of spatial and angular frequencies given in the manual ISO 8608. A study is then performed on vehicle and occupant dynamic responses for various combinations of sitting occupancies, velocities and road profiles. The results obtained underscore the need for considering sitting occupancies in addition to velocity and road profile for assessment of ride comfort for a vehicle.