I have a
SpatialPolygonsDataFrame object in R with about 3000 polygons and associated data. I have columns of expected and observed cases. For the moment, I don't have any covariates.
I'd like to perform spatial smoothing. I have used functions such as
EBlocal() from the
DCluster package but have not been happy with the results.
Therefore I would like to do a conditional autoregression using
spautolm() - I would like to model a Poisson model of observed versus expected with neighbour weighting.
I cannot figure out how to specify such a model. How do I specify merely a Poisson model for the
spautolm command with no covariates?
Or am I barking up the wrong tree altogether?