search for: errorsarlm

Hi all, I am trying to run an autologistic model using the function errorsarlm from spdep package. **I built an XY matrix extracting the two colums from matriz** coords1<-matriz[matriz$casos1==1, c(4,5)] coords1<-as.matrix(coords1) **I identify neighbours of region points** nb20<-dnearneigh(coords1,0,20,longlat=TRUE) ** I build a neighbours list with spatial we...

...I am trying to extract an estimator for the bread of the sandwich function. I used bread(fitted model) however it seems that I have missed something as an error message "no applicable method for "bread" appears. My fitted model is a Spatial simultaneous autoregressive error model.(errorsarlm in spdep package) Can anyone please tell me what I might be doing wrong? Your help is much appreciated. Thank you, Oshadhi Samarasinghe. RA Department of Economics University of Auckland New Zealand

Dear list, (also sent to Roger Bivand, but perhaps somebody of you can help me also) I am trying to use package spdep for fitting an SAR model with errorsarlm. However, I am not sure how to make a valid nb object out of my neighborhood. As far as I have seen, there is no documentation for nb.object. I have done the following: class(pschmid$nb) <- "nb" # pschmid is a prab object as generated from function prabinit, package # prabclus, and t...

...following code resulted in an error: > load("Panel.90s.ok.R") > attach(Panel.90s.ok) > neighs<-dnearneigh(cbind(x,y),0,50000) > help(nbdists) > dists<-nbdists(neighs,cbind(x,y)) > Weights<-nb2listw(neighs,zero.policy=T,glist=dists) > error.model.rprice<-errorsarlm(log(wells+1)~log(year)+log(area)+log(lag.rprice )+log(marketable+1)+log(prov.rds+1)+log(ind.rds+1)+log(seismic+1)+log(pipeline s+1)+log(min.year)+log(max.year)+log(avg.year),listw=Weights,zero.policy=T) Error in eigen(w, only.values = TRUE) : error code 2311 from Lapack routine dgeev The prob...

Greetings to the list, I was trying to estimate spatial error model in R, somehow I got the message below. Would you please help me with it? Many thanks in advance. Error in solve.default(asyvar, tol = tol.solve) : system is computationally singular: reciprocal condition number = 5.6964e-18 Regards, Dong

Hello, Has anyone modified the errorsarlm in the R package spdep to allow for more than a single spatial autocorrelation coefficient (i.e. 'lambda')? Or, if not, any initial suggestions on how to make that modification? I have looked at the source code for the function and realize that any attempt to do it on my own w...

...quot;SBA~", V) + rd <- round(nrow(data)*(2/3)) + d <- sample(seq(1:nrow(data)),rd) + dat1 <- data[d,] + dat2 <- data[-d,] + crd <- cbind(dat1$Longitude,dat1$Latitude) + dist80 <- dnearneigh(crd,0,100,longlat=F) + dist80sw <- nb2listw(dist80, style="B") + fm <- errorsarlm(as.formula(V), data=dat1, listw=dist60sw) + pred <- predict(fm,dat2) + C[i,1] <- cor(dat2$SBA,pred) + out <- cbind(C) + } + colMeans(out) + } > > sim.sp(df2007.5k.s2,CM,1,1000) Error in nb2listw(dist80, style = "B") : Empty neighbour sets found I guess I got error message...

...(data)*(2/3)) d <- sample(seq(1:nrow(data)),rd) dat1 <- data[d,] dat2 <- data[-d,] crd <- cbind(dat1$Longitude,dat1$Latitude) dist60 <- dnearneigh(crd,0,60,longlat=F) dist60sw <- nb2listw(dist60, style="B") fm <- errorsarlm(as.formula(V), data=dat1, listw=dist60sw) pred <- predict(fm,dat2) C[i,1] <- cor(dat2$SBA,pred) out <- cbind(C) } colMeans(out) } sim.sp(df2007.5k.s2,CM,1,1000) Error in nb2listw(dist60, style = "B") : Empty neighbour sets found And, here are my...

I'm trying to use the spdep package to calculate this: y = rho W y + e I don't want to use explanatory variables, just the lag from the dependent variable. How would I code this?

...95-3081833.post at n4.nabble.com> Content-Type: text/plain; charset=us-ascii This has been answered offlist (the poster also wrote directly to me as package maintainer, but did not post on the R-sig-geo list, as would have seemed natural). The resolution was to read ?formula, and to use either errorsarlm() or lagsarlm() in spdep with formula=y ~ 1. Apparently an insurance analyst in a hurry ... Roger Saswati Neogi wrote: > > > > I'm trying to use the spdep package to calculate this: > > y = rho W y + e > > I don't want to use explanatory variables, just the l...