[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