[R-sig-Geo] can we use georob package to make a 3D nested variogram?
Bingwei Tian
bwtian at gmail.com
Fri Dec 5 08:56:53 CET 2014
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.
More information about the R-sig-Geo
mailing list