R structures

Dimensions

In a traditional analysis, the residuals are IID (Independently and Identically Distributed) in which case the order does not matter. However, ASReml is designed to handle cases where residuals are correlated. Sometimes, these correlations can be defined as a direct product of variance structures, assuming the residuals are in an appropriate order.

IID residuals

The variance structure can be written as 
residual units
or 
residual idv(units)

This analysis will estimate any variance components as ratios relative to the implicit residual variance unless !SIGMAP is specified.

Typical field trial - regular rectangular grid

The variance structure can be written as 
 ar1(row).ar1(column)
This particular coding assumes there is a factor in the data called row which indexes the rows of the field layout, and another factor column which indexes the columns of the field layout, and that that there are data records for all cells in the complete row x column grid. (It will be necessary to estimate any missing values in the data.)
AR1 is the VCODE for a first order Auroregressive matrix. It requires one parameter (correlation) given an initial value of 0.1.

If the data file does not contain row and column information, but a correlated structure is required, the data records need to be sorted on row and column (even though the coding is not in the file) and the structures (assuming 12 rows and 36 columns) would be written as 
residual ar1(12).ar1(36)

Spatial data - irregular pattern - Matern function

The variance structure can be written as 
residual mat1(fac(X,Y))
MAT1 is the VCODE for a MATERN correlation structure with 1 parameter estimated. It requires one parameter (range) given an initial value of 10. The fac(X,Y) model function collects the plot coordinates for variables X and Y for the matern function to use.

Multivariate data

The R structure (assuming 3 traits) can be written as 
residual units.us(Trait)
The data is sorted traits within records so the structure for traits is specified second. Again we let ASReml supply the number of records.
US is the VCODE for an unstructed variance matrix. It requires six initial values: 3 variances and 3 covariances in the order V11 C21 V22 C31 C32 V22 (lower triangle rowwise). However, it is generally difficult to guess suitable values and ASReml will obtain initial values as a proportion of the simple variances and covariances of the residual.

See Also