Displaying 2 results from an estimated 2 matches for "tsport_ac".
Did you mean:
tsport_ace
2011 Jan 05
1
Prediction error for Ordinary Kriging
...lt;- ceiling(diff(t(bb))/cs)
rsa2_grd <- GridTopology(cellcentre.offset = cc,cellsize = cs, cells.dim = cd)
getClass("SpatialGrid")
p4s <- CRS(proj4string(rsa2))
x2_SG <- SpatialGrid(rsa2_grd, proj4string = p4s)
x2_SP <- SpatialPoints(cbind(x2$X, x2$Y))
v <- variogram(log1p(tsport_ace) ~ 1, x2, cutoff=100, width=9)
te<- fit.variogram(v,vgm(0.0437, "Exp", 26, 0))
y <- krige(tsport_ace~1, x2, x2_SG, model = ve.fit)
spplot(y, 1, col.regions = bpy.colors(100), sp.layout = list("sp.lines",as(rsa2,
"SpatialLines"),no.clip = TRUE))
I'm lookin...
2011 Jan 06
1
Cross validation for Ordinary Kriging
ear ALL,
The last part of my thesis analysis is the cross validation. Right now I am
having difficulty using the cross validation of gstat. Below are my commands
with the tsport_ace as the variable:
nfold <- 3
part <- sample(1:nfold, 69, replace = TRUE)
sel <- (part != 1)
m.model <- x2[sel, ]
m.valid <- x2[-sel, ]
t<- fit.variogram(v,vgm(0.0437, "Exp", 26, 0))
cv69 <- krige.cv(tsport_ace ~ 1, x2, t, nfold = nrow(x2))
The last line gives an err...