In order to understand the current status of the predictive performance of spallation models and issues towards further improvement, a comprehensive study of the proton-induced spallation product yields was conducted. In this study, four kinds of the latest spallation models (i.e. INCL++/ABLA07, INCL++/GEMINI++, INCL4.6/GEM, and CEM03.03), utilized in Monte Carlo particle transport calculation, are employed, and six pieces of cross section data measured at GSI, RIKEN, and J-PARC are used for comparison. As a result, the INCL++/ABLA07 calculations are in general agreement with the GSI and RIKEN experimental data, whereas some discrepancies are observed especially in the comparisons with the J-PARC experiments. Finally, several key issues inherent in each spallation model are described.