Calculates the percentage variance and sum of squares accounted for in the strata of an ANOVA analysis (R.W. Payne).
Options
PRINT = string tokens |
Controls printed output (%variance, %ss); default %var |
METHOD = string token |
Which terms to take into account (treatments, blocksandtreatments); default trea |
%VARIANCE = pointer |
Saves the percentage variance accounted for |
%SS = pointer |
Saves the percentage sum of squares accounted for |
SAVE = ANOVA save structure |
Save structure for the analysis; by default this will be the most recent ANOVA |
No parameters
Description
Procedure A%VARIANCE calculates the percentage of the variance and of the sum of squares that are accounted for in the strata of an analysis of variance performed by the ANOVA directive. By default, this will be the most recent ANOVA analysis. However, you calculate these for an earlier analysis, by using the SAVE option of A%VARIANCE to specify its save structure (saved using the SAVE parameter of the earlier ANOVA command).
By default only the treatment terms estimated in the stratum concerned are included when accounting for its variance and sum of squares. For example, suppose we have a randomized-block design with block structure Blocks / Plots = Blocks + Blocks.Plots and a treatment term Treat estimated within blocks (i.e. in the Blocks.Plots stratum). The default, METHOD=treatments, assesses only how Treat accounts for the Blocks.Plots variation. Alternatively, if you set METHOD=blocksandtreatments, any block terms marginal to each stratum are also included. In this example the term Blocks is marginal to Blocks.Plots, and these blocks will have been defined in order to decrease the amount of within-block variation. So we may be want to use this setting to see how effective that has been.
Output is controlled by the PRINT option, with settings: %variance prints the percentage variance accounted for, and %ss prints the percentage sum of squares accounted for. By default PRINT=%variance. Options: PRINT, METHOD, %VARIANCE, %SS, SAVE. Parameters: none.
Method
The percentage variance accounted for is calculated as
100 × (total mean square – residual mean square) / total mean square
and the percentage sum of squares accounted for is similarly calculated as
100 × (total sum of squares – residual sum of squares) / total sum of squares
For both methods, the residual mean squares and residual sums of squares are taken from the residual lines of the relevant strata in the original analysis of variance. To account only for treatments, A%VARIANCE does an analysis with the original block structure but no treatments. The totals are given by the sums of squares and mean squares of the relevant strata in that analysis. To account for marginal block terms as well as treatments, A%VARIANCE takes the sums of squares and mean squares for each stratum from an analysis containing only that stratum. (So no other block terms will have been fitted.)
See also
Directive: ANOVA.
Commands for: Analysis of variance.
Example
CAPTION 'A%VARIANCE example',!t('Split plot design, see the',\
'Guide to Genstat, Part 2, Section 4.2.1.'); STYLE=meta,plain
FACTOR [NVALUES=72; LEVELS=6] Blocks
& [LEVELS=3] Wplots
& [LEVELS=4] Subplots
GENERATE Blocks,Wplots,Subplots
FACTOR [LABELS=!T('0 cwt','0.2 cwt','0.4 cwt','0.6 cwt')] Nitrogen
& [LABELS=!T(Victory,'Golden rain',Marvellous)] Variety
VARIATE Yield; DECIMALS=2; EXTRA=' of oats in cwt. per acre'
READ [SERIAL=yes] Nitrogen,Variety,Yield
4 3 2 1 1 2 4 3 1 2 3 4 3 1 2 4 4 1 2 3 2 1 3 4
2 3 4 1 4 2 3 1 1 4 2 3 3 4 1 2 1 3 4 2 2 3 4 1
4 1 3 2 3 4 1 2 3 4 2 1 3 1 4 2 4 3 1 2 1 2 3 4 :
3 3 3 3 1 1 1 1 2 2 2 2 3 3 3 3 1 1 1 1 2 2 2 2
2 2 2 2 3 3 3 3 1 1 1 1 3 3 3 3 2 2 2 2 1 1 1 1
2 2 2 2 1 1 1 1 3 3 3 3 1 1 1 1 2 2 2 2 3 3 3 3 :
156 118 140 105 111 130 174 157 117 114 161 141
104 70 89 117 122 74 89 81 103 64 132 133
108 126 149 70 144 124 121 96 61 100 91 97
109 99 63 70 80 94 126 82 90 100 116 62
96 60 89 102 112 86 68 64 132 124 129 89
118 53 113 74 104 86 89 82 97 99 119 121 :
" Convert yields to cwt per acre."
CALCULATE Yield=(Yield*80)/(112*4)
" Define the treatment structure: factorial effects of V and N."
TREATMENTS Variety*Nitrogen
" Subplots nested within whole-plots nested within blocks."
BLOCK Blocks/Wplots/Subplots
ANOVA [PRINT=means] Yield
A%VARIANCE [PRINT=%variance,%ss]
& [METHOD=blocksandtreatments]