can we use georob package to make a 3D nested variogram?
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>