R/rbf.cv.R

In geospt: Geostatistical Analysis and Design of Optimal Spatial Sampling Networks

assign("rbf.cv",
       function(formula, data, eta, rho, n.neigh, func){
         z = extractFormula(formula, data, newdata=data)$z
         rbf.pred <- as.data.frame(matrix(NA,nrow= length(z), ncol=4))
         colnames(rbf.pred) <- c("x","y","var1.pred","var1.var")
         for(i in 1:length(z)){
           rbf.pred[i,3] <- rbf(formula, data[-i,], eta, rho, newdata=data[i,], n.neigh, func)[,3]
         }
         RMSPE <-  sqrt(sum((rbf.pred$var1.pred-z)^2)/length(z))
         RMSPE
       }
)