Defines the residual term for a
REML analysis, or the residual term for an experiment within a meta-analysis (combined analysis of several experiments).
||Level of the
||Model term to be used as the residual|
||Whether the structure is formed by direct product construction or by definition of the whole matrix (
||Allows an initial estimate to be provided for the residual variance of the experiment|
||Allows the residual variance to be fixed at its initial value (
||Coordinates of the data points to be used in calculating distance-based models|
||Type of covariance model associated with the term(s), or with individual factors in the term(s) if
||Order of model|
||Heterogeneity for correlation matrices (
||How to calculate distances when
||Factors over which to form direct products|
||To define matrix values for the term or the factors when
||To define values for matrix inverses (instead of the fixed matrices themselves) when
||Initial parameter values for each correlation matrix|
||Texts containing strings
||Non-zero values in the variate indicate groups of parameters whose values are to be constrained to be equal|
VRESIDUAL is used to define the residual term for a
REML analysis or to define separate residual terms for different experiments within a multi-experiment (or meta-) analysis. The
TERM option is used to specify the formula for the residual term. This term need not have been specified previously by the
For a single experiment,
VRESIDUAL can be used to impose a covariance structure on the residual term. This could also be done by specifying the covariance structure using
VRESIDUAL has the advantage that the algorithm then checks that the correct residual term is used.
In a multi-site experiment,
VRESIDUAL can be used to specify a different residual model for each separate experiment. The
EXPERIMENT option is used to define the experiment(s) for which the model is to be used.
VARIANCE option is used to give an initial value for the residual variance in the current experiment(s). You can set option
CONSTRAINT=fix to fix the residual variance at the initial value rather than estimating it (as a positive value).
The definition of the residual terms then follows mainly as for the definition of correlated error terms through
VSTRUCTURE The exception is that power models can be defined only in terms of the coordinates of the data points, not by specifying coordinates for the factor levels. (So the
COORDINATES parameters of
VSTRUCTURE are not present in
For a multi-experiment analysis, the factors and variates for the separate experiments should be concatenated into structures which run over all the experiments. For example, consider an experiment set up at two sites to compare a set of 24 varieties in four replicates. In one site the experiment was laid out as a grid of eight rows by 12 columns, in the other a grid of 16 rows by six columns was used. In these circumstances, a single set of factors (of length 192) can be used to specify the design, using factors to describe variety, rows and columns, plus a factor
expt defining the allocation of units to experiments.
VCOMPONENTS [FIXED=Variety; EXPERIMENTS=Expt]
VRESIDUAL [EXPERIMENT=1; TERM=Row.Col] MODEL=AR,AR;\
VRESIDUAL [EXPERIMENT=2; TERM=Row.Col] MODEL=AR,AR;\
Where some factors differ between experiments, these should be defined on the units relevant to the appropriate experiment(s) and missing elsewhere. When an
EXPERIMENTS factor has been defined, the default action of the
MVINCLUDE option of
REML is changed to include units with missing y-values and missing factor levels.
Commands for: REML analysis of linear mixed models.