Displaying 1 result from an estimated 1 matches for "cv69".
Did you mean:
c769
2011 Jan 06
1
Cross validation for Ordinary Kriging
...ifficulty 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 error saying:
Error in SpatialPointsDataFrame(coordinates(data),
data.frame(matrix(as.numeric(NA), :
row.names of data and coords do not match
I don't know what is wrong. The x2 data is a SpatialPointsdatafram...