Plots models fitted to 2-dimensional auto- and cross-variograms (D.A. Murray).
Options
PLOT = string token |
Controls how to display the plotted variograms (separate, scattermatrix); default scat |
|---|---|
ESTIMATES = pointer |
Pointer containing model estimates saved from MCOVARIOGRAM |
Parameter
COVARIOGRAM = pointer |
Pointer to supply the semi-variances, distances and associated information as saved from FCOVARIOGRAM |
|---|
Description
DCOVARIOGRAM plots 2-dimensional auto- and cross-variograms using data generated by FCOVARIOGRAM. DCOVARIOGRAM can also be used to display the fitted model for isotropic models using estimates generated from MCOVARIOGRAM. The data should be supplied in a pointer that has been saved using the COVARIOGRAM option from FCOVARIOGRAM. This pointer provides the auto-variograms, cross-variograms and associated information required for the plots. Its elements contain:
1 a matrix with columns of variograms and cross-variograms and rows indexed by lags within directions;
2 a variate of counts at the lags in each direction;
3 distances of the lags in each direction;
4 horizontal angles;
5 vertical angles;
6 variances;
7 distance classes;
8 method;
9 pointer containing identifiers of the DATA variates;
10 number of dimensions.
The ESTIMATES option can be used to plot an isotropic fitted model of coregionalization where the estimates are taken directly from MCOVARIOGRAM. Graphical output is controlled using the PLOT option. The setting separate produces each auto- and cross-variogram on a separate plot. Alternatively, they can be combined onto a single scatter matrix using the scattermatrix setting.
Options: PLOT, ESTIMATES.
Parameter: COVARIOGRAM.
Method
DCOVARIOGRAM uses Genstat’s standard graphics and calculation commands. See MCOVARIOGRAM directive for details on the models and estimates used for plotting the model of coregionalization.
Reference
Webster, R. & Oliver, M.A. (2001). Geostatistics for Environmental Scientists. Wiley, Chichester.
See also
Directives: FCOVARIOGRAM, MCOVARIOGRAM, COKRIGE.
Procedure: DVARIOGRAM.
Commands for: Spatial statistics, Graphics.
Example
CAPTION 'DCOVARIOGRAM example',\
!t('Data from survey at CEDAR farm - south-east England.',\
'Exchangeable potassium, available phosphorus and yield',\
'of barley measured from topsoil of a 6.4 ha field.',\
'Data supplied with kind permission from Dr Z.L.Carroll');\
STYLE=minor,plain
VARIATE [NVALUES=160] k; !(138,159,175,203,154,149,130,119,123,101,113,\
146,145,164,187,167,142,144,113,108,149,119,135,137,124,143,153,140,\
140,135,135,131,120,117,117,129,152,151,137,147,132,133,143,127,137,\
125,137,135,122,160,122,138,124,131,123,157,140,152,136,145,182,195,\
207,147,163,162,165,161,146,144,239,230,215,177,181,152,140,123,137,\
119,222,243,207,204,188,186,189,155,158,137,197,230,223,172,177,187,\
174,134,146,127,151,160,179,164,151,146,135,138,130,131,160,168,160,\
170,169,145,143,120,138,139,151,149,139,181,154,140,153,158,155,147,\
160,176,173,186,164,158,161,139,125,142,110,161,170,168,152,151,129,\
180,162,168,168,161,190,194,206,222,239,190,183,152)
VARIATE [NVALUES=160] p;!(33.8,41,53.2,45.8,42.8,44.6,45.6,45.6,42,\
42.6,35.4,35.8,30.4,41.4,42.6,49.8,44.4,43.6,37.8,39,48.8,31.6,41,\
42.6,35,39.2,42.8,58.2,48,46.6,34.6,33.2,21,20.4,18.8,21.8,20,19.8,\
23.2,27.8,30,29.4,30,22.8,24.4,16.8,23.4,21.8,24.2,22.4,34,43.8,40.2,\
35.6,27.6,35,32.4,39.2,45,43.6,41.8,41.2,41.4,50.4,54,51.8,49.8,49.6,\
38.8,37,63.8,67.4,63.6,52.2,46.2,41.6,39.2,39.8,48.2,40.8,63.2,67.8,\
71.8,69,58,55.8,55.6,50.2,56.2,51.8,77.4,89,76.6,66.4,60.4,60,60.2,\
54.8,51.2,58.8,70.4,61.4,61.6,51.6,57,51,47.2,50.8,48.2,49.6,66,72.8,\
60,68.8,75.6,49.6,48.8,51.4,58.4,56,62.4,64,58.6,51,56.2,59.2,59.6,\
59.6,63.2,62.8,71.6,64.4,54.6,63.8,57,67.6,58.6,56.8,68.2,62.6,54.8,\
62,61.2,56,46,53.6,57.8,52.4,55.2,52.8,62.2,65.2,60,59.8,73.4,63.4,\
66,57.4,58.2,57)
VARIATE [NVALUES=160] yield; !(3.69,3.39,3.35,3.26,3.57,4.43,3.18,2.43,\
3,3.19,3.7,3.72,3.8,2.4,2.58,2.94,*,3.79,3.52,2.9,2.28,3.15,3.3,3.58,\
2.96,3.28,3.16,2.58,1.86,3.08,*,3.43,3.47,3.41,3.43,3.86,3.69,3.28,3.13,\
2.18,3.38,3.49,3.49,3.5,3.04,3.41,3.25,3.28,3.32,3.06,2.97,3.05,3.75,\
3.2,3.46,3.2,3.15,3.22,3.62,3.13,2.31,3.03,3.75,3.2,3.46,3.2,3.15,3.22,\
3.62,3.13,2.91,3.13,2.34,2.61,3.11,3.23,3,2.85,2.37,2.9,2.19,1.96,2.61,\
2.56,1.95,2.8,3.04,2.94,3.38,3.65,1.81,1.39,1.28,2.33,2.34,2.95,2.73,\
2.47,2.67,3.31,3.67,3.14,2.96,2.33,1.75,2.92,2.6,2.59,2.9,2.89,3.78,\
3.6,2.84,3.24,2.85,2.88,2.95,3.11,3.36,3.11,3,2.96,2.52,2.79,2.87,3,\
2.96,2.52,2.68,2.36,3.68,3.35,3.4,3.2,3.02,3.14,2.23,2.74,2.6,3.1,3.25,\
3.43,2.83,3.06,3.14,2.82,2.77,2.58,3.03,2.14,3.36,3.39,3.29,3.36,3.21,\
3.03,3.38,3.14,2.81,3.87)
VARIATE [NVALUES=160] x; VALUES=!(10(100,120...400))
VARIATE [NVALUES=160] y; VALUES=!((100,120...280)16)
VARIATE lags; VALUES=!(40,60...200)
FCOVARIOGRAM [PRINT=statistics,variograms; MAXLAG=200; STEPLENGTH=lags;\
DIRECTIONS=!(90); SEGMENTS=180; MAXCONE=500; MINCOUNT=1;\
METHOD=common; COVARIOGRAM=save_cov] DATA=k,p,yield; X1=x; X2=y
DCOVARIOGRAM save_cov
MCOVARIOGRAM [PRINT=estimates; MAXLAG=200; COVARIOGRAM=save_cov;\
MINCOUNT=2; WEIGHT=equal] MODEL=nugget,spherical;\
INITIAL=*,130; ESTIMATES=save_est
DCOVARIOGRAM [ESTIMATES=save_est] save_cov