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