Diagnostic Residual statistics and plots from spatial analyses
increases the amount of information reported on the residuals
obtained from the analysis of a two dimensional regular grid field trial.
is written to the
The variogram has been suggested as a useful diagnostic for assisting
with the identification of appropriate variance models for spatial
data (Cressie, 1991). Gilmour et al. (1997) demonstrate its
usefulness for the identification of the sources of variation in the
analysis of field experiments. If the elements of the data vector (and
hence the residual vector) are indexed by a vector of spatial
coordinates, then the ordinates of the
sample variogram are given by the half the variance of the difference
between the residuals.
The sample variogram reported by ASReml has two forms depending on
whether the spatial coordinates represent a complete rectangular
lattice (as typical of a field trial) or not.
In the lattice case,
the sample variogram is calculated from the triple (Li1,Lli2,vi) where
and Li2 are the displacements in the two directions. As there
will be many vi with the same displacements,
ASReml calculates the means for each displacement pair Li1, Li2
either ignoring the signs (default) or separately
for same sign and opposite sign ( !TWOWAY),
after grouping the larger displacements: 9-10, 11-14, 15-20, ....
The result is displayed
as a perspective plot of the one or two surfaces
by absolute displacement group.
this case, the two directions may be on different scales.
Otherwise ASReml forms a variogram based on radial coordinates.
It calculates the distance between points di=
sqrt(Li12+ Li22) and an angle
thetai subtended by the line from (0,0) to
with the x-axis (-180
|lt thetai ||lt 180).
The angle can be calculated as thetai =
||lt thetai||lt 180) since there is radial
symmetry (the numbers are the same for ij
The variogram presented averages the vi within
12 distance classes and 4, 6 or 8
sectors (selected using a
centred on an angle of (i-1)*180/s (i=1,...s).
A figure is produced which
reports the trends in average(vi) with increasing distance for each sector.
computes the variogram from predictors of random effects which appear to have a
variance structures defined in terms of distance.
The variogram details are reported in the .res file.
modifies the appearance of the variogram calculated from
the residuals obtained from a two dimensional field trial.
The default form is based on absolute 'distance' in each direction.
This form distinguishes same sign and different sign distances
and plots the variances separately as two layers in the same figure.
requests that the variogram formed from the BLUPs
predicted for a
term be based on s (=4, 6 or 8)
sectors of size 180/s. The appropriate sectors are
centred on the X and Y directions. The details are written
file as well as creating a figure.
Return to start