# NOT RUN {
require(spData)
data(boston)
id_obs <-sample(dim(boston.c)[1],400)
######################### data at observed sites
y <- log(boston.c[id_obs,"MEDV"])
x <- boston.c[id_obs, c("CRIM", "INDUS","ZN","NOX","AGE")]
coords <- boston.c[id_obs , c("LON", "LAT") ]
######################### data at predicted sites
x0 <- boston.c[-id_obs, c("CRIM", "INDUS","ZN","NOX", "AGE")]
coords0 <- boston.c[-id_obs , c("LON", "LAT") ]
mod <- scgwr( coords = coords, y = y, x = x )
pred0 <- predict0( mod=mod, coords0=coords0, x0=x0)
pred <- pred0$pred # predicted value
b <- pred0$b # spatially varying coefficients
b[1:5,]
bse <- pred0$bse # standard error of the coefficients
bt <- pred0$t # t-values
bp <- pred0$p # p-values
# }
Run the code above in your browser using DataLab