Message-ID: <519E7A1B.4000207@uni-muenster.de>
Date: 2013-05-23T20:20:43Z
From: Edzer Pebesma
Subject: fit.variogram.reml
In-Reply-To: <CAOytA4ZetcNdio3hKW6WjsaqHtaWk-gvd8jBdfR5tZsbrMHkVg@mail.gmail.com>
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