Hi Edzer,
Thanks for the clue.
However, I am still confused and probably doing something wrong. Using the
dX-option, I expect a reduction in the number of pairs for the semivariogram
fitting. In my case with n=120 and three categories a reduction from
0.5*120*(120-1)= 7140 to 3*{0.5*40*(40-1)}=2340. Any suggestions?
points(fractie): 'result.dat', x=1, y=2, v=3, X=4&5&6,dX=0.1, mv=-99;
# X=1: x, X=2:y, X=3: fraction, X=4: binary cat.1, X=5: binary cat 2, X=6:
binary cat.3 , [X=7: number, X=8: lost]
variogram(fractie): 2 Nug(0) + 1 Sph(60); # initial spherical variogram
set fit=3;
set cutoff=90;
set width=12;
set output='result1.out';
set debug=5;
set logfile='result1.log';
Paula Westerman
Wageningen Agricultural University
Group Crop and Weed Ecology
[The C.T. de Wit Graduate School]
Haarweg 333
Postbus 430
6700 AK Wageningen
tel: +31 (0)317 484769 /482536
fax: +31-(0)317-484892 /485572
- - - - - - - - - - - - - - Original Message - - - - - - - - - - - - - -
check out the dX=... option, as in
data(x): 'file', ... , X=-1&3&4&5, dX=0.5;
--
Edzer
- - - - - - - - - - - - End of Original Message - - - - - - - - - - - -
