R/rbf.tcv.R

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

assign("rbf.tcv",
   function(formula, data, eta, rho, n.neigh, func){
   z = extractFormula(formula, data, newdata=data)$z
   s = coordinates(data)
   rbf.pred <- as.data.frame(matrix(NA,nrow= length(z), ncol=8))
   colnames(rbf.pred) <- c("var1.pred","var1.var","observed","residual","zscore","fold","x","y")
   for(i in 1:(length(z))){
   rbf.pred[i,1] <- rbf(formula, data[-i,], eta, rho, newdata=data[i,], n.neigh, func)[,3]
   rbf.pred[i,6] <- i
}
rbf.pred[,3]<- z
rbf.pred[,7:8]<-s
rbf.pred[,4]<- rbf.pred[,3]-rbf.pred[,1]
rbf.pred
}
)