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

[Damian Maddalena]

Thank you for your reply, Paulo.

As long I can find a statistically valid method to compare the Bayesian 
kriging surfaces to the others I will be happy. I don't necessarily need 
CV. It was just what I knew for the other interpolation methods.

Please forgive my simple question, but can you please suggest a few 
terms that I might use to search for "others measures of fit?" Something 
that I could search to find some good code examples.  I am fairly new to 
coding in R and having something to examine and modify would be a big help.

Thank you.

-Damian




On 03/13/2013 01:37 PM, Paulo Justiniano Ribeiro Jr wrote:
> 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
>>