The VARIOGRAM Procedure

Suppose you specify NDIRECTIONS=3 in the COMPUTE statement in PROC VARIOGRAM. This results in three angle classes defined by midpoint angles between and : , , and , where is the angle tolerance. If you do not specify an angle tolerance by using the ANGLETOLERANCE= option in the COMPUTE statement, the following default value is used:

For example, if NDIRECTIONS=3, the default angle tolerance is . When the directed line segment in Figure 19 is superimposed on the coordinate system that shows the angle classes, its angle is approximately , measured clockwise from north. In particular, it falls within , the second angle class (Figure 20).

Note: If the designated points and are labeled in the opposite order, the orientation is in the opposite direction—that is, approximately instead of approximately . This does not affect angle class selection; the angle classes and are the same.

Figure 20: Selected Pair Falls within the Second Angle Class

If you specify an angle tolerance less than the default, such as ATOL=15, some point pairs might be excluded. For example, the selected point pair in Figure 20, while closest to the axis, might lie outside . In this case, the point pair would be excluded from the semivariance computation. This setting can be desirable if you want to reduce interference between neighboring angles. An angle tolerance that is too small might result in too few point pairs in some distance classes for the empirical semivariance estimation (see also the discussion in the section Choosing the Size of Classes).

On the other hand, you can specify an angle tolerance greater than the default. This can result in a point pair being counted in more than one angle classes. This has a smoothing effect on the variogram and is useful when only a small amount of data is present or the available data are sparsely located. However, in cases of anisotropy the smoothing effect might have the side effect of amplifying weaker anisotropy in some direction and weakening stronger anisotropy in another (Deutsch and Journel 1992, p. 59).

Changes in the values of the BANDWIDTH= option have a similar effect. See the section Bandwidth Restriction for an explanation of how BANDWIDTH= functions.

An alternative way to specify angle classes and angle tolerances is with the DIRECTIONS statement. The DIRECTIONS statement is useful when angle classes are not equally spaced. When you use the DIRECTIONS statement, consider specifying the angle tolerance too. The default value of the angle tolerance is when a DIRECTIONS statement is used instead of the NDIRECTIONS= option in the COMPUTE statement. This might not be appropriate for a particular set of angle classes. See the section DIRECTIONS Statement for more details.