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

2014-12-09 Thread Edzer Pebesma
On 12/09/2014 02:13 PM, Bingwei Tian wrote: > Yes, you are right, only one nugget effect can be fitted. > I fitted the variogram with one nugget, but still did not better than eye > fitted. > > If your criterion for "better" is visual, then you must be right. If your criterion involves minimis

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

2014-12-09 Thread Bingwei Tian
Yes, you are right, only one nugget effect can be fitted. I fitted the variogram with one nugget, but still did not better than eye fitted. - Bingwei Ph.D. Student Kyoto University C-1-2-225, Katsura Campus, Kyoto University, Nishikyo-ku, 〒615-8530, Kyoto, Japan -- View this message in co

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

2014-12-05 Thread 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 Nu

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

2014-12-04 Thread Bingwei Tian
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=

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

2014-12-04 Thread Edzer Pebesma
On 12/05/2014 05:08 AM, Bingwei Tian wrote: > Dear list, > It seems only gstat can make a 3D nested variogram but can not auto fit it > until now, . ?fit.variogram gives me Description: Fit ranges and/or sills from a simple or nested variogram model to a sample variogram Or do you m

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

2014-12-04 Thread Bingwei Tian
Dear list, It seems only gstat can make a 3D nested variogram but can not auto fit it until now, . Recently I find georob is also a useful package to do 3D estimation with trend data, But I did not find the way to create a 3D nested variogram, Does anyone kown about this, or any other package for