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
[EMAIL PROTECTED] [EMAIL PROTECTED]
303 472-3240 (cell) <-best choice
303 716-1617
