[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