[R-sig-Geo] xvalid or similar with krige.bayes()

[Paulo Justiniano Ribeiro Jr]

Well this is not that straightfoward (or uniquely defined)
since you need to specify what would 
be the "point prediction" in the Bayesian kriging.
The posterior mean is a candidate, but not uniquely defined.
Besides the analysis will be quite timing consuming
when running krige.bayes() for each of the location left aside

I personally would go for other measures of fit
but if you realy find CV within the Bayesian approach you can write
a simple loop call to the leaving-one-out



Paulo Justiniano Ribeiro Jr
LEG (Laboratorio de Estatistica e Geoinformacao)
Universidade Federal do Parana
Caixa Postal 19.081
CEP 81.531-990
Curitiba, PR  -  Brasil
Tel: (+55) 41 3361 3573
VOIP: (+55) (41) (3361 3600) 1053 1066
Fax: (+55) 41 3361 3141
e-mail: paulojus AT  ufpr  br
http://www.leg.ufpr.br/~paulojus

On Tue, 12 Mar 2013, Damian Maddalena wrote:

> Greetings:
>
> I need to perform LOOCV for krige.bayes().  I want to compare the output of 
> krige.bayes to the output of other interpolation methods. I can do LOOCV 
> easily for other interpolation methods using xvalid() in geoR or krige.cv() 
> in gstat but I cannot find any code examples that demonstrate how this can be 
> done for a krige.bayes() object.
>
> Any help would be greatly appreciated.
>
> Thank you.
>
> -Damian
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>