[R] GepR - bayesian modelling - please help!
Monica Palaseanu-Lovejoy
monica.palaseanu-lovejoy at stud.man.ac.uk
Fri Apr 23 11:45:41 CEST 2004
Hi,
I am coming back with a question about GeoR hopping somebody
can help.
I am doing a bayesian kriging and the prediction coordinates are
shifted 0.1 m from the data coordinates.
I am using the following command and i am getting the following
error:
zn.bayes <- krige.bayes(zn.gdata, loc = loc1, model =
model.control(cov.model = "exponential", lambda = 0),
prior = prior.control(phi.prior ="exponential", phi =
89.1948), output=output.control(n.predictive=2,
mean.var = TRUE, quantile = c(0.025, 0.25, 0.5, 0.75,
0.975), threshold = c(300)))
Error: Error in cond.sim(env.loc = base.env, env.iter = iter.env,
loc.coincide = get("loc.coincide", :
chol: matrix not pos def, diag[13]= -3.487549e-019
What i am doing wrong? The location loc1 is a data frame with the
shifted coord. values. Any suggestion is very well appreciated.
Thanks,
Monica
More information about the R-help
mailing list