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]