[R-sig-Geo] fit.variogram.reml

[Paulo Justiniano Ribeiro Jr]

I'm unsure the REML's in gstat and geoR reffer to the same method.
In the former the REML is based on the sample variogram
whereas for the latter the REML is applied to the assumed
(Gaussian) distribution of the data, in a modification of
the original likelihood function.



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 Thu, 23 May 2013, Edzer Pebesma wrote:

> Saman, package geoR contains functions for (RE)ML fitting that also fit
> the range parameter.
>
> Please note that in your output below the nugget is negative, which
> can't be the case -- one could argue that this is a software error.
>
> On 05/23/2013 10:02 PM, Saman Monfared wrote:
>> Hi,
>> In gstat and  package the fit.variogram.reml function provided for
>> RMLE method to  estimation of variogram parameters. This function only
>> estimate nugget and sill and range parameter is fix.
>>
>> for example:
>>
>> v.rml<-fit.variogram.reml(rate~1,data,model =vgm(200000,"Sph",95241.99,0))
>> v.rml
>>   model     psill    range
>> 1   Nug -1711.462     0.00
>> 2   Sph  3431.108 95241.99
>>
>> How can I set the best range value?
>> If  exist a way to estimate range like( fit.ranges = TRUE) in
>> fit.variogram function?
>>
>> Best,
>>
>
> -- 
> Edzer Pebesma
> Institute for Geoinformatics (ifgi), University of Münster
> Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
> 83 33081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>