We investigate earthquake-induced landslides using a geostatistical model featuring a latent spatial effect (LSE). The LSE represents the spatially structured residuals in the data, which remain after adjusting for covariate effects. To determine whether the LSE captures the residual signal from a given trigger, we test the LSE in reproducing the pattern of seismic shaking from the distribution of seismically induced landslides, without prior knowledge of the earthquake being included in the model. We assessed the landslide intensity, i.e., the expected number of landslides per mapping unit, for the area in which landslides triggered by the Wenchuan and Lushan earthquakes overlap. We examined this area to test our method on landslide inventories located in near and far fields of the earthquake. We generated three models for both earthquakes: i) seismic parameters only (proxy for the trigger); ii}) the LSE only; and iii) both seismic parameters and the LSE. The three configurations share the same morphometric covariates. This allowed us to study the LSE pattern and assess whether it approximated the seismic effects. Our results show that the LSE reproduced the shaking patterns for both earthquakes. In addition, the models including the LSE perform better than conventional models featuring seismic parameters only. Due to computational limitations we carried out a detailed analysis for a relatively small area (2112 km2), using a dataset with higher spatial resolution. Results were consistent with those of a subsequent analysis for a larger area (14648 km2) using coarser resolution data.