Produces graphs and statistics for bivariate analysis of variance (R.F.A. Poultney).
|Controls printing of statistics from the bivariate analysis
||Controls output from the (univariate)
||Treatment terms to be fitted in the analysis of variance; this option must be set|
||Block model defining the error terms in the analysis of variance; if unset, the design is assumed to be unstratified (i.e. to have a single error term)|
||Single model term identifying the treatment term whose means are to be plotted|
||Stratum from which to extract treatment information; default is to take the bottom stratum|
||Limit on number of factors in a treatment term; default 3|
||Significance level to use in the calculation of the radius of the confidence region and the region of non-significance; default 0.95|
||Type of graphical output (
||controls the style of axes in a high-resolution graph (
||Plotting symbols for the means; default is to take the letters A to Z, then a to z|
||First variate for the bivariate analysis|
||Second variate for the bivariate analysis|
||Title for the graph|
ABIVARIATE produces a bivariate analysis of variance with a graphical representation of the results, as described by Dear & Mead (1983, 1984). The procedure was developed from a Genstat 4 macro, further information about which is given by Poultney & Riley (1986), and is intended primarily for data from intercropping experiments. The variates to be analysed (specified by parameters
Y2) are measurements, usually yields, taken on the two crops. The final parameter,
TITLE, defines a title for the graph.
The procedure will work for any of the designs that can be analysed by
ANOVA, except that there must be no pseudo-factors. Option
TREATMENTSTRUCTURE defines the treatment formulae for the analysis, and the block formula is defined by the
BLOCKSTRUCTURE can be omitted if there is a single error stratum (i.e. the analysis is of a completely randomized design). The
FACTORIAL option controls the number of factors in each treatment term, as in the
First of all,
ABIVARIATE calculates a univariate analysis of variance for each of the variates
Y2, with output controlled by the
APRINT option. The settings are the same as those in the
ANOVA directive; by default
Output from the bivariate analysis of variance, which follows, is controlled by the
error generates the error summary statistics from the bivariate analysis: Error Sum of Products, Variances after Adjustment for Covariance, Correlation Coefficient between
Y2, Radius of Standard Errors, Radius of Confidence Regions, and Radius of Non-Significance Regions. The setting
treatment produces the following statistics for each treatment term estimated within the specified error stratum: Treatment Sum of Products, Wilks’ Lambda, Bivariate F-Statistic.
The stratum from which the means (and other information) are to be taken is defined by
STRATUM option; if this is omitted, the lowest stratum is used. The significance level to use in the calculation of confidence regions is defined by the
PROBABILITY option; by default this is 0.95.
TERM option specifies a treatment term whose means are to be represented graphically. The means are plotted on axes transformed to allow for the variability in, and the correlation between, each crop variate. The plotting symbols can be defined as a factor or text using the option
LABELS. Alternatively they will be taken to be the first n values of the series A to Z, a to z where n is the number of means to be plotted. The graph can be either line printer or high resolution, the default being high resolution. The external axes of a high-resolution graph can be suppressed by setting
Problems arise in situations where the table of means to be plotted is incomplete; this can occur when a whole factor level is restricted out, or where the treatment structure is nested within a control. The length of the vector
LABELS is calculated as the number of cells in the table, including missing values. If
LABELS is declared, it must have length equal to the dimension of the table otherwise a fault will occur. Similarly, the calculation of the radius statistics is based on the assumption that the table of means is complete and has equal replication. These values, if printed, would be incorrect for a table with missing cells and so are suppressed. They can be calculated by hand as shown by Dear & Mead (1983).
(1) calculate the SSP matrix for all terms in the formula
(2) transform the variables such that the new set are uncorrelated and have unit error variance
(3) calculate new axes based on the maximum and minimum points of the transformed variables
(4) draw the graph of the transformed means with the axes rotated such that they are at the same angle to the vertical
Y2 can be restricted, however this restriction must be identical for the two variates. Some problems may occur when whole levels of factors are restricted out leaving empty cells in the table of means to be plotted (see above).
Dear, K.B.G. & Mead, R. (1983). The use of bivariate analysis techniques for the presentation, analysis and interpretation of data. Statistics in Intercropping Technical Report No. 1. Department of Applied Statistics, University of Reading, U.K.
Dear, K.B.G. & Mead, R. (1984). Testing assumptions and other topics in bivariate analysis. Statistics in Intercropping Technical Report No. 2. Department of Applied Statistics, University of Reading, U.K.
Poultney, R.F.A. & Riley, J. (1986). A Genstat Macro for the Bivariate Analysis of Intercropping Data. Genstat Newsletter, 17, 27-46
Commands for: Analysis of variance.
CAPTION 'ABIVARIATE example',\ !t('Data are from a melon and okra intercropping trial;',\ 'for further details see Poultney & Riley',\ '(1986, Genstat Newsletter 17).'); STYLE=meta,plain FACTOR [LABELS=!t(M1,M2); VALUES=8(1,2)3] Melon FACTOR [LABELS=!t(O1,O2); VALUES=4(1,2)6] Okra FACTOR [LABELS=!t(year1,year2,year3); VALUES=16(1,2,3)] Year FACTOR [LEVELS=4; VALUES=(1...4)12] Block FACTOR [LEVELS=12; VALUES=4(1...12)] Plot TEXT [VALUES='1','2','3','*','+','#'] Tsymbol VARIATE [NVALUES=48] MelonDat,OkraDat READ [SERIAL=yes] OkraDat,MelonDat 1.06 3.28 2.18 1.98 2.52 1.96 1.56 1.67 5.58 3.08 2.35 2.07 1.90 2.40 2.59 4.10 5.4 6.2 7.0 5.8 5.6 7.8 8.3 6.9 4.4 4.0 5.0 4.2 7.1 6.6 9.2 4.8 3.13 2.30 2.21 2.30 2.81 3.30 4.81 5.80 2.70 4.34 2.50 3.61 4.70 4.90 2.30 2.44 : 63.4 29.0 45.4 46.4 44.6 42.7 50.2 45.9 37.0 53.9 56.9 39.0 53.1 47.9 28.3 20.1 32.0 20.0 20.2 19.0 9.5 16.5 6.8 18.5 18.0 20.5 12.8 14.0 21.0 16.5 5.0 10.0 16.6 22.9 15.1 16.7 10.3 34.3 13.3 12.6 20.5 22.7 14.8 11.8 20.1 19.4 23.3 21.1 : ABIVARIATE [TREAT=Melon*Okra*Year; Block=Block/Plot; TERM=Okra.Year;\ LABEL=Tsymbol] OkraDat; MelonDat; TITLE='melon and okra data'