Hello,

I am trying to compare the zscores of using various kriging methods with inverse distance weighting.

When I run

crossval_idw <- krige.cv(Prec00hr~1, DataCoordhr, nfold=2)

I get all NA's for the zscore. As the residual is a valid number, I am assuming that the standard deviation is not being calculated properly. However, I can calculate the standard deviation of the residual:

sd(crossval_idw$residual)

As an aside, I also get all NA's for the variances even though, I also can calculate the variances separately.

var(crossval_idw$residual)

Any ideas?

Thank you for your help,

Cara

	[[alternative HTML version deleted]]

