We develope statistical methods for the analysis of multivariate abundance data in ecology. We extend the generalized linear latent variable models to handle spatial and temporal correlations inherent in data and they will be applied to two existing datasets from community ecology, e.g. to study the effect of restoration on Finnish peatland vegetation.