PROC VARIOGRAM features automated semivariogram fitting. In particular, the procedure selects a theoretical semivariogram model to fit the empirical semivariance and produces estimates of the model parameters in addition to a fit plot. You have the option to save these estimates in an item store, which is a binary file format that is defined by the SAS System and that you cannot modify. Then, you can retrieve this information at a later point from the item store for future analysis with PROC KRIGE2D or PROC SIM2D.

The coal seam thickness empirical semivariogram in Figure 109.7 shows first a slow, then rapid, rise from the origin. This behavior suggests that you can approximate the empirical semivariance with a Gaussian-type form

as shown in the section Theoretical Semivariogram Models. Based on this remark, you choose to fit a Gaussian model to your classical semivariogram. Run PROC VARIOGRAM again and specify
the MODEL
statement with the FORM=
GAU option. By default, PROC VARIOGRAM uses the weighted least squares (WLS) method to fit the specified model, although you
can explicitly specify the METHOD=
option to request the fitting method. You want additional information about the estimated parameters, so you specify the
CL
option in the MODEL
statement to compute their 95% confidence limits and the COVB
option of the MODEL
statement to produce a table with their approximate covariances. You also specify the STORE
statement to save the fitting outcome into an item store file with the name `SemivStoreGau`

and a desired label. You run the following statements:

proc variogram data=sashelp.thick outv=outv; store out=SemivStoreGau / label='Thickness Gaussian WLS Fit'; compute lagd=7 maxlag=10; coordinates xc=East yc=North; model form=gau cl / covb; var Thick; run; ods graphics off;

After you run the procedure you get a series of output objects from the fitting analysis. In particular, Figure 109.10 shows first a model fitting table with the name and a short label of the model that you requested to use for the fit. The table also displays the name and label of the specified item store.

If you specify no parameters, as in the current example, then PROC VARIOGRAM initializes the model parameters for you with default values based on the empirical semivariance; for more details, see the section Theoretical Semivariogram Model Fitting. The initial values provided by the VARIOGRAM procedure for the Gaussian model are displayed in the table in Figure 109.11.

Otherwise, in PROC VARIOGRAM you can specify initial values for parameters with the PARMS statement. Alternatively, you can specify fixed values for the model scale and range with the SCALE= and RANGE= options, respectively, in the MODEL statement. A nugget effect is always used in model fitting. Unless you explicitly specify a fixed nugget effect with the NUGGET= option in the MODEL statement or initialize the nugget parameter in the PARMS statement, the nugget effect is automatically initialized to zero. See the section Syntax: VARIOGRAM Procedure for more details about how the MODEL statement and the PARMS statement handle model parameters.

The output in Figure 109.12 comes from the optimization process that takes place during the model parameter estimation. The optimizer produces an optimization information table, information about the optimization technique that is used, optimization-related results, and notification about the optimization convergence.

Figure 109.12: Fitting Optimization Information

Spatial Correlation Analysis with PROC VARIOGRAM |

The VARIOGRAM Procedure

Dependent Variable: Thick

Angle: Omnidirectional

Current Model: Gaussian

Dual Quasi-Newton Optimization |

Dual Broyden - Fletcher - Goldfarb - Shanno Update (DBFGS) |

Hessian Computed by Finite Differences (Using Analytic Gradient) |

Optimization Results | |||
---|---|---|---|

Iterations | 12 | Function Calls | 45 |

Gradient Calls | 0 | Active Constraints | 1 |

Objective Function | 11.433894152 | Max Abs Gradient Element | 3.0128744E-8 |

Slope of Search Direction | -3.986332E-8 |

The fitting process is successful, and the parameters converge to the estimated values shown in Figure 109.13. For each parameter, the same table also displays the approximate standard error, the degrees of freedom, the t value, the approximate p-value, and the requested 95% confidence limits.

The approximate covariance matrix of the estimated parameters is displayed in Figure 109.14.

The fitting summary table in Figure 109.15 displays statistics about the quality of the fitting process. In particular, the table shows the weighted error sum of squares
in the `Weighted SSE`

column and the Akaike information criterion in the `AIC`

column. See more information about the fitting criteria in section Quality of Fit.

Figure 109.16 demonstrates the fitted theoretical semivariogram against the empirical semivariance estimates with the weighted least squares method. The fit seems to be more accurate closer to the origin , and this is explained as follows: A smaller corresponds to smaller semivariance; in turn, this corresponds to smaller semivariance variance, as shown in the section Theoretical and Computational Details of the Semivariogram. By definition, the WLS optimization weights increase with decreasing variance, which leads to a more accurate fit for smaller distances in the WLS fitting results.