Susceptibility analysis for predicting earthquake-induced landslides has most frequently been done using deterministic methods; multivariate statistical methods have not previously been applied. In this study, however, we introduce a statistical methodology that uses the intensity of earthquake shaking as a landslide triggering factor. This methodology is applied in a study of shallow earthquake-induced  landslides in central western Taiwan. The results show that we can accurately interpret landslide distribution in the study area and predict the occurrence of landslides in neighboring regions. This susceptibility model is capable of predicting shallow landslides induced during an earthquake scenario with similar range of ground shaking, without requiring the use of geotechnical, groundwater or failure depth data.