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