[R-sig-Geo] Me{spdep} error and results interpretation

[Roger Bivand]

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