Saves experiment × treatment means from analysis of a series of trials by
VASERIES (R.W. Payne).
||Limit on the number of factors in the terms generated from the
||Saves residual variances from the experiments|
||Definitions of the models used by
||Terms for which means are to be saved|
||Experiment × term tables of means|
||Experiment × term tables of standard errors of means|
||Average standard errors of differences for the experiments|
VASMEANS saves experiment-by-term tables of means, and associated information, from a series of experiments analysed by the
VASERIES procedure. The model definitions and
REML save structures from the analysis can be specified by the
SAVE options. If either of these is not specified, the output is taken from the most recent
The parameters of
VASMEANS save the information about the fixed terms. The
TERMS parameter specifies a model formula, which Genstat expands to form the series of terms about which you wish to save information. The
FACTORIAL option sets a limit on the number of factors in each term. Any term containing more than that limit is deleted.
The subsequent parameters allow you to specify identifiers of data structures to store the information about each of the terms that you have specified. The
MEANS parameter saves tables of predicted means, classified by the experiment factor and the factors of the term. Similarly the
SEMEANS parameter saves tables of standard errors. The
AVESEDMEANS parameter saves tables of average standard errors for the means, classified by the experiment factor (as required for the
VMETA procedure). If you have a single term, you can supply a table for each of these parameters, as appropriate. However, if you have several terms, you must supply a pointer which will then be set up to contain as many tables as there are terms.
RESIDUALVARIANCES option saves residual variances from the experiments, in a table classified by the experiment factor (also as required for the
CAPTION 'VASMEANS example'; STYLE=meta SPLOAD [PRINT=*] '%gendir%/examples/WheatTrials.gsh' " find best random model for each trial " VASERIES [PRINT=best,devi,aic,sic,dfrandom,summary; FIXED=Entry;\ EXPERIMENTS=Trial; ROWS=Row; COLUMNS=Column; BLOCKS=Block;\ MVINCLUDE=yvariate; TRYSPATIAL=ifregular; TRYTRENDS=yes;\ RSTRATEGY=fastoptimal; VCONSTRAINTS=positive]\ Yield; MODELDEFINITIONS=modeldefs; SAVE=save " save trial x entry means, residual variances and average SEDs " VASMEANS [RESIDUALVARIANCES=Sigma2] Entry; MEANS=Means;\ SEMEANS=SE; AVESEDMEAN=AverageSEDs PRINT Means,SE PRINT Sigma2,AverageSEDs