Re: minor bug in conditional simulation on 3d grid

["Edzer J. Pebesma" <e.pebesma@geog.uu.nl>]

The question is: should we consider a stratified approach with only
one stratum as stratified? If you leave out the "s=4" part in data(),
the result will be correct; if you specify, besides MN, more
variables, the result will also be correct.

Anyway, I corrected this; see source code.

Best regards,
--
Edzer

bob sandefur wrote:

> hi-
> The cmd file:
>
> data(MN): 'jnkxyzvinmn.eas', x=1, y=2, z=3, v=6, min=3, max=12,radius=120;
> data(): 'jnkpau.eas', x=1, y=2, ,z=3, s=4 ;
> variogram(MN): 0.066 Nug() + 0.22 Sph(80) + 0.27 Sph(278) + 0.348057
> Sph(929.455);
> method: gs;
> set output='jnksimpau.out';
> set nsim=1;
>
> Gives the following in jnksimpau.out
> $ head jnksimpau.out
> prediction results from file: jnkpau.eas variable: col[0]
> 5
> x
> y
> z
> s
> pred(MN)
> var(MN)
> 6275 12781 4017.5 1 0.231736 NA
> 7595 11551 4182.5 1 -0.733279 NA
>
> So the field count is off by 1 (should be 6)
>
> Regards
>
> Robert (Bob) L. Sandefur PE
> Geostatistical Associate
> CAM
> 200 Union Suite G-13
> Lakewood, Co
> 80228
>
> rsandefur@cam-llc.com
> rlsandefur@aol.com
>
> 303 472-3240 (cell) <-best  choice
>
> 303 716-1617
>
>
>