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 <milansharma2012@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]]

