[SOLVED] Kribayes function issues

[MaximeBeauchamp]

Hi everyone,

I am running a cokriging with external drift (2 variables) according to the following command:
grid.db = kriging(data.db, grid.db, v_fit, neigh = neigh.create(ndim = 2, type = 0),
uc = c("1", "f1", "f2"), radix = "CK")

No problem with that. I wanted to run additional tests with the kribayes function, just changing the previous kriging call by:

grid.db = kribayes(data.db, grid.db, v_fit, neigh = neigh.create(ndim = 2, type = 0),
uc = c("1", "f1", "f2"),
dmean = c(0,1,0),
dcov = matrix(c(0.001, 0, 0, 0, 0.001, 0, 0, 0, 0.001), nrow = 3),
radix = "CK")

but I got this error :
Abort : Wrong number of arguments in RC_get_value_double (3 > 1)
Error in kribayes(data.db, grid.db, v_fit, neigh = neigh, uc = uc, radix = "CK") :
Abort caught by C and return to R
Calls: call_kriging -> kribayes -> .Call
Exécution arrêtée

Do you know what is this error?

Thanks for your help,

Maxime

[Fabien Ors]

Thank you for reporting.
Currently bayesian kriging is not available in multivariate context.
This feature will be added to the next RGeostats version.