[R-sig-Geo] can we use georob package to make a 3D nested variogram?

[Edzer Pebesma]

Well, yes, you can't fit two nuggets independently, can you? Either fix
one of them, or remove one of them. Note that these two are different
models, when it comes to fitting parameters:

> vgm(1, "Sph", 900)
  model psill range
1   Sph     1   900
> vgm(1, "Sph", 900, 0)
  model psill range
1   Nug     0     0
2   Sph     1   900


On 12/05/2014 08:56 AM, Bingwei Tian wrote:
> Thanks edzer, As you said , I hope it can choose the number of
> nested structures, at least keep the eyefitted number and the corresponding
> models automatically.
> But, for my case , the fit seems do noting except change the second Nug to 
> 1.
> 
> R:> uk.eye1  <- vgm(psill = 0.155,  model = "Gau",  range=700,  nugget=0.01)
> R:> uk.eye   <- vgm(psill = 0.125,  model = "Sph",  range=35000,  nugget=0, 
> add.to=uk.eye1)
> R:> uk.eye
>   model psill range
> 1   Nug 0.010     0
> 2   Gau 0.155   700
> 3   Nug 0.000     0
> 4   Sph 0.125 35000
> 
> R:> uk.fit  <- fit.variogram(uk.vgm, uk.eye, fit.sills = TRUE, fit.ranges =
> TRUE,
> + fit.method = 7, debug.level = 1, warn.if.neg = FALSE)
> 
> Warning: singular model in variogram fit
> 
> R:> uk.fit
>   model psill range
> 1   Nug 0.010     0
> 2   Gau 0.155   700
> 3   Nug 1.000     0
> 4   Sph 0.125 35000
> 
> 
> 
> -----
> Bingwei
> 
> Ph.D. Student 
> Kyoto University
> C-1-2-225, Katsura Campus, Kyoto University, 
> Nishikyo-ku, 〒615-8530, Kyoto, Japan
> --
> View this message in context: http://r-sig-geo.2731867.n2.nabble.com/can-we-use-georob-package-to-make-a-3D-nested-variogram-tp7587522p7587524.html
> Sent from the R-sig-geo mailing list archive at Nabble.com.
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> 

-- 
Edzer Pebesma,  Co-Editor-in-Chief Computers & Geosciences
Institute for Geoinformatics (ifgi), University of Münster
Heisenbergstraße 2, 48149 Münster, Germany. Phone: +49 251
83 33081 http://ifgi.uni-muenster.de GPG key ID 0xAC227795

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 490 bytes
Desc: OpenPGP digital signature
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20141205/c32c18a0/attachment.bin>