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

[Edzer Pebesma]

Paulo, they do refer to the same method. gstat::fit.variogram.reml only
fits sill parameters (nugget, partial sill), but doesn't compute nor use
sample variograms.

On 05/24/2013 04:14 AM, Paulo Justiniano Ribeiro Jr wrote:
> 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
>>

-- 
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