Me{spdep} error and results interpretation
Please do not post in HTML, and do follow the posting guide by providing a reproducible example from a standard data set. There is nothing anyone can do when the most likely cause of your problems is your own data and choice of methods.
On Thu, 25 Jul 2013, Amy Henry wrote:
Hi,
I fitted my spatial data (count data and environmental variables) to a
glm.nb model. I decided to detect and correct for spatial
autocorrelation using ME{spdep}. I used family poisson since with
negative binomial the model did not converge. I received this warning
after I added the eigenvectors into the model. Can you help me
understand what the error means? Also, should I infer from the example
below that spatial autocorrelation is negligible since the pr(ZI) values
are so high?
tu1<-glm.nb(total.aeg~Total.number.of.units, data=inhdb, offset(x))
tuME<-ME(total.aeg~Total.number.of.units, data=inhdb, family="poisson", offset(x), listw=sw.nb, alpha=0.5)
tuME
Eigenvector ZI pr(ZI) 0 NA NA 0.40 1 43 NA 0.39 2 19 NA 0.39 3 18 NA 0.36 4 34 NA 0.35 5 23 NA 0.41 6 25 NA 0.39 7 56 NA 0.36 8 51 NA 0.47 9 46 NA 0.41 10 36 NA 0.32 etc.....
inhdb$eigen_43<-fitted(tuME)[,1]
inhdb$eigen_19<-fitted(tuME)[,2]
tuglmME<-glm.nb(total.aeg~Total.number.of.units+fitted(tuME), data=inhdb, offset(x))
There were 27 warnings (use warnings() to see them) Warning messages: 1: In theta.ml(Y, mu, sum(w), w, limit = control$maxit, trace = control$trace > ... : iteration limit reached Thank you for your help! Best, Amy [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Roger Bivand Department of Economics, NHH Norwegian School of Economics, Helleveien 30, N-5045 Bergen, Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43 e-mail: Roger.Bivand at nhh.no