fit.variogram.reml
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