Hello,

I am performing kriging with an external drift and ordinary kriging in a loop to have predicted values at each time step.

The predicted values are changing reasonably, however, the variance is not changing after the first time step.

Here is my simple for loop:

# KED
list_ked <- vector("list",72)

for(i in 1:72){
 list_ked[[i]] <- krige(Prec00[,i]~Z, locations=DataCoord, newdata=elev, model = fitted_vario)
}

# OK
list_ok <- vector("list",72)

for(i in 1:72){
     list_ok[[i]] <- krige(Prec00[,i]~1, DataCoord, newdata=coord_points, model = fitted_vario)
}

Would anyone know the solution?

Thank you very much,

Cara

	[[alternative HTML version deleted]]

