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

[Edzer Pebesma]

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