AI-GEOSTATS: Cross Validation

Email to send questions: ai-geostats@jrc.it

AI-GEOSTATS: Cross Validation

Hello, I am trying to compare various kriging methods. I was hoping someone could explain what is wrong with my code to perform leave-one-out cross validation. I am iterating through precipitation time steps to get cross validation each time step. For ordinary kriging I am using the formula Pre

Re: AI-GEOSTATS: Cross Validation

Tobin, this is hard to find out without a reproducible example, and without knowing how the individual objects were created. A wild guess is the following for(i in 1:72){ DataCoord$y = Prec00[,i] crossval_ok[[i]] <- krige.cv(y~1, DataCoord, model = fitted_vario) } You do assume that the list un

Re: AI-GEOSTATS: Cross Validation

Hi Cara, You dont really give enough information here, but assuming that you want to build a formula from the column names of Prec00, I think the problem is your formula definition. This can for example be created as: as.formula(paste(names(Prec00)[i],"~1",sep="")) giving: crossval_ok[[i]] <-