Displaying 1 result from an estimated 1 matches for "dist60".
2009 Jul 15
1
Error in simulation R-code
...x)
V <- paste(V, collapse="+")
V <- paste("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)
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)
}
s...