Spatial Analysis of Multienvironment trials
Generating/Checking spatial coding
Consider the simultaneous analysis of a series of trials
with similar treatments and layout. For field trials,
the residuals are independent between trials but may be spatially
within trials. At the treatment level there is the
average effect of treatments across trials and the interaction
of treatment with trial.
For analysis in ASReml, the data for each trial needs to be
arranged conformably in a single or
sorted by trial. The data also needs to be sorted into plot order
but this can be done by ASReml if the rows and columns
are coded in the data file. Assume we have a data field
which identifies the trial, and fields
which code for the spatial arrangement within trials.
!SECTION site !ROWFAC row !COLFAC column
is specified, ASReml generates the R-structure lines to fit
a separable first-order autoregressive (
AR1 x AR1
) variance structure for each site.
is the name of a factor or variate containing column numbers (1 ...
nc where nc is the number of columns) on which the data is to
is the name of a factor or variate containing row numbers (1 ...
nr where nr is the number of rows) on which the data is to
name of a factor or variate containing site numbers (1 ...
ns where ns is the number of sites) on which the data is already sorted.
is specified ASReml will check that sections have been correctly dimensioned
but assumes the data records are presorted with respect to
are also specified, ASReml will construct the R-structure lines
required to model the spatial variance as AR1 x AR1
with initial correlation parameters of 0.1,
such that the data is sorted into plot order: columns within rows within sites.
The R structure lines that a user would normally be
required to work out and type into the
are written to the
file. The user may then cut and paste them into the
file for a later run if the structures need to be
The number of rows and columns within a trial does not
have to be the same across trials.
The following is a basic example assuming 5 sites (sections).
Syntax Basic multi-environment trial analysis
site 5 # sites coded 1 ... 5
column * # columns coded 1 ...
row * # rows coded 1 ...
variety !A # variety names
met.dat !SECTION site !ROWFAC row !COLFAC col
yield ~ site !r variety site.variety !f mv
site 2 0 # variance header line
# asreml inserts the 10 lines required to
# define the R structure lines for the five sites
Return to start