R/idwST.cv1.R

In geosptdb: Spatio-Temporal Radial Basis Functions with Distance-Based Methods (Optimization, Prediction and Cross Validation)

assign("idwST.cv1",
function (param, formula, data, n.neigh, progress=FALSE)
{
    C <- param[1]
    factor.p <- param[2]
    z = extractFormula(formula, data, newdata = data)$z
    idw.pred <- as.data.frame(matrix(NA, nrow = length(z), ncol = 5))
    colnames(idw.pred) <- c("x", "y", "t", "var1.pred", "var1.var")
    if(progress)
      pb <- txtProgressBar(min = 0, max = length(z), char = "=", style = 3)
    for (i in 1:length(z)) {
      if(progress)
      setTxtProgressBar(pb, i)
      idw.pred[i, 4] <- idwST(formula, data[-i, ], newdata = data[i, ], n.neigh=n.neigh, C=C, factor.p=factor.p, progress=F)[, 4]
    }
    if(progress)
      close(pb)
    RMSPE <- sqrt(sum((idw.pred$var1.pred - z)^2)/length(z))
    RMSPE
}
)