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@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
