In this paper we deal with the problem of variable selection in spatiotemporal autoregressive (STAR) models with neighbourhood effects. We propose a procedure to carry out the selection process, taking into account the uncertainty associated with estimation of the parameters and the predictive behaviour of the compared models, in order to give more realism to the analysis. We set up a multi-objective programming problem that combines the use of different criteria to measure both these aspects. We use genetic algorithms which are very flexible and suitable for our multicriteria decision problem. In particular, the procedure allows us to estimate the number of spatial and temporal nearest neighbours as well as their relative effects. The methodology is illustrated through an application to the real estate market of Zaragoza. Copyright (C) 2010 John Wiley & Son, Ltd.STAR , variable selection , genetic algorithms , neighbourhood effects ,