naginterfaces.library.anova.random¶
- naginterfaces.library.anova.random(y, iblock, nt, it, tol, irdf)[source]¶
random
computes the analysis of variance and treatment means and standard errors for a randomized block or completely randomized design.For full information please refer to the NAG Library document for g04bb
https://www.nag.com/numeric/nl/nagdoc_29.2/flhtml/g04/g04bbf.html
- Parameters
- yfloat, array-like, shape
The observations in the order as described by and .
- iblockint
Indicates the block structure.
There are no blocks, i.e., it is a completely randomized design.
There are blocks and the data should be input by blocks, i.e., must contain the observations for block followed by the observations for block , etc.
There are blocks and the data is input in parallel with respect to blocks, i.e., must contain the first observation for block , must contain the first observation for block must contain the first observation for block must contain the second observation for block , etc.
- ntint
The number of treatments. If only blocks are required in the analysis then set .
- itint, array-like, shape
Note: the required length for this argument is determined as follows: if : ; otherwise: .
indicates which of the treatments plot received, for .
If , is not referenced.
- tolfloat
The tolerance value used to check for zero eigenvalues of the matrix . If a default value of is used.
- irdfint
An adjustment to the degrees of freedom for the residual and total.
The degrees of freedom for the total is set to and the residual degrees of freedom adjusted accordingly.
The total degrees of freedom for the total is set to , as usual.
- Returns
- gmeanfloat
The grand mean, .
- bmeanfloat, ndarray, shape
If , contains the mean for the th block, , for .
- tmeanfloat, ndarray, shape
If , contains the (adjusted) mean for the th treatment, , for , where is the mean of the treatment adjusted observations, .
- tablfloat, ndarray, shape
The analysis of variance table. Column 1 contains the degrees of freedom, column 2 the sum of squares, and where appropriate, column 3 the mean squares, column 4 the -statistic and column 5 the significance level of the -statistic. Row 1 is for Blocks, row 2 for Treatments, row 3 for Residual and row 4 for Total. Mean squares are computed for all but the Total row; -statistics and significance are computed for Treatments and Blocks, if present. Any unfilled cells are set to zero.
- cfloat, ndarray, shape
If , the upper triangular part of contains the variance-covariance matrix of the treatment effects, the strictly lower triangular part contains the standard errors of the difference between two treatment effects (means), i.e., contains the covariance of treatment and if and the standard error of the difference between treatment and if , for , for .
- irepint, ndarray, shape
If , the treatment replications, , for .
- rfloat, ndarray, shape
The residuals, , for .
- effloat, ndarray, shape
If , the canonical efficiency factors.
- Raises
- NagValueError
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, no blocks or treatments in model.
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, .
Constraint: .
- (errno )
On entry, and .
Constraint: and if , must be a multiple of .
- (errno )
On entry, at least one treatment is not present. Treatment does not appear in .
- (errno )
On entry, , and .
Constraint: .
- (errno )
On entry, the values in are constant.
- (errno )
The computation of the eigenvalues has failed to converge.
- (errno )
A computed standard error is zero.
- Warns
- NagAlgorithmicWarning
- (errno )
The treatments are totally confounded with blocks.
- (errno )
The residual degrees of freedom is zero.
- (errno )
The residual mean square is zero.
- (errno )
The design is disconnected.
- Notes
In the NAG Library the traditional C interface for this routine uses a different algorithmic base. Please contact NAG if you have any questions about compatibility.
In a completely randomized design, experimental material is divided into a number of units, or plots, to which a treatment can be applied. In a randomized block design the units are grouped into blocks so that the variation within blocks is less than the variation between blocks. If every treatment is applied to one plot in each block it is a complete block design. If there are fewer plots per block than treatments then the design will be an incomplete block design and may be balanced or partially balanced.
For a completely randomized design, with treatments and plots per treatment, the linear model is
where is the th observation for the th treatment, is the overall mean, is the effect of the th treatment and is the random error term. For a randomized block design, with treatments and blocks of plots, the linear model is
where is the effect of the th block and the notation indicates that the th treatment is applied to the th plot in the th block.
The completely randomized design gives rise to a one-way analysis of variance. The treatments do not have to be equally replicated, i.e., do not have to occur the same number of times. First the overall mean, , is computed and subtracted from the observations to give . The estimated treatment effects, are then computed as the treatment means of the mean adjusted observations, , and the treatment sum of squares can be computed from the sum of squares of the treatment totals of the divided by the number of observations per treatment total, . The final residuals are computed as , and, from the residuals, the residual sum of squares is calculated.
For the randomized block design the mean is computed and subtracted from the observations to give . The estimated block effects, ignoring treatment effects, , are then computed using the block means of the and the unadjusted sum of squares computed as the sum of squared block totals for the divided by number of plots per block, . The block adjusted observations are then computed as . In the case of the complete block design, with the same replication for each treatment within each block, the blocks and treatments are orthogonal, and so the treatment effects are estimated as the treatment means of the block adjusted observations, . The treatment sum of squares is computed as the sum of squared treatment totals of the divided by the number of replicates to the treatments, . Finally the residuals, and hence the residual sum of squares, are given by .
For a design without the same replication for each treatment within each block the treatments and the blocks will not be orthogonal, so the treatments adjusted for blocks need to be computed. The adjusted treatment effects are found as the solution to the equations
where is the vector of the treatment totals for block adjusted observations, , is a diagonal matrix with equal to the number of times the th treatment is replicated, and is the incidence matrix, with equal to the number of times treatment occurs in block . The solution to the equations can be written as
where is a generalized inverse of . The solution is found from the eigenvalue decomposition of . The residuals are first calculated by subtracting the estimated treatment effects from the block adjusted observations to give . However, since only the unadjusted block effects have been removed and blocks and treatments are not orthogonal, the block means of the have to be subtracted to give the correct residuals, and residual sum of squares.
The mean squares are computed as the sum of squares divided by the degrees of freedom. The degrees of freedom for the unadjusted blocks is , for the completely randomized and the complete block designs the degrees of freedom for the treatments is . In the general case the degrees of freedom for treatments is the rank of the matrix . The -statistic given by the ratio of the treatment mean square to the residual mean square tests the hypothesis
The standard errors for the difference in treatment effects, or treatment means, for the completely randomized or the complete block designs, are given by:
where is the residual mean square and in the complete block design. In the general case the variances of the treatment effects are given by
from which the appropriate standard errors of the difference between treatment effects or the difference between adjusted means can be calculated.
In the complete block design all the information on the treatment effects is given by the within block analysis. In other designs there may be a loss of information due to the non-orthogonality of treatments and blocks. The efficiency of the within block analysis in these cases is given by the (canonical) efficiency factors, these are the nonzero eigenvalues of the matrix , divided by the number of replicates in the case of equal replication, or by the mean of the number of replicates in the unequally replicated case, see John (1987). If more than one eigenvalue is zero then the design is said to be disconnected and some treatments can only be compared using a between block analysis.
- References
Cochran, W G and Cox, G M, 1957, Experimental Designs, Wiley
Davis, O L, 1978, The Design and Analysis of Industrial Experiments, Longman
John, J A, 1987, Cyclic Designs, Chapman and Hall
John, J A and Quenouille, M H, 1977, Experiments: Design and Analysis, Griffin
Searle, S R, 1971, Linear Models, Wiley