First thing that comes to mind would be to suggest you to use a different method inside the function lagsarlm. In the spdep manual it says
"spam" or "Matrix_J" for strictly symmetric weights lists of styles "B" and "C" and since you have a weight matrix of style B. That could fix your problem. Juan On Sat, Dec 29, 2012 at 1:43 PM, Milan Sharma <milansharma2...@yahoo.com>wrote: > > > Hi all, > can any one of you see whats wrong in the code below? > > model<-lm(Tcrash~X100MVMT+Alcohol+Alignment+PVMT, data=data)#works > > sar<-spautolm(Tcrash~X100MVMT+Alcohol+Alignment+PVMT,listw=WI_queenW, > data=data)#works > > lag <- lagsarlm(Tcrash~X100MVMT+Alcohol+Alignment+PVMT, data=data, > listw=listwW)#shows error > Error in solve.default(inf, tol = tol.solve) : > system is computationally singular: reciprocal condition number = > 2.44638e-12 > > lag <- lagsarlm(Tcrash~X100MVMT+Alcohol+Alignment+PVMT, data=data, > listw=WI_queenW) > Error in solve.default(inf, tol = tol.solve) : > system is computationally singular: reciprocal condition number = > 1.51647e-13 > > WI_queenW > Characteristics of weights list object: > Neighbour list object: > Number of regions: 71 > Number of nonzero links: 364 > Percentage nonzero weights: 7.22079 > Average number of links: 5.126761 > > Weights style: B > Weights constants summary: > n nn S0 S1 S2 > B 71 5041 364 728 8160 > ............................. > Milan > [[alternative HTML version deleted]] > > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > > -- Juan Tomás Sayago Gómez Graduate Research Assistant West Virginia University - RRI 886 Chestnut Ridge Road, Room 520 P.O. Box 6825 Morgantown, WV 26506-6825 [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo