naginterfaces.library.anova.contrasts

naginterfaces.library.anova.contrasts(tmean, irep, rms, rdf, ct, tabl, tol, usetx, tx)

contrasts computes sum of squares for a user-defined contrast between means.

For full information please refer to the NAG Library document for g04da

https://support.nag.com/numeric/nl/nagdoc_30.3/flhtml/g04/g04daf.html

Parameters

tmeanfloat, array-like, shape $\left(\text{nt}\right)$

The treatment means, ${\hat{\tau }}_{\text{i}}$, for $\text{i}=1,2,\dots ,t$.

irepint, array-like, shape $\left(\text{nt}\right)$

The replication for each treatment mean, $n_{\text{i}}$, for $\text{i}=1,2,\dots ,t$.

rmsfloat

The residual mean square, ${\hat{\sigma }}^2$.

rdffloat

The residual degrees of freedom.

ctfloat, array-like, shape $(\text{nt},\text{nc})$

The columns of $ct$ must contain the $\text{nc}$ contrasts, that is $ct[\text{i}-1,\text{j}-1]$ must contain ${\lambda }_{\text{i}\text{j}}$, for $\text{j}=1,2,\dots ,\text{nc}$, for $\text{i}=1,2,\dots ,t$.

tablfloat, array-like, shape $(\text{nc},5)$

The elements of $tabl$ that are not referenced as described below remain unchanged.

tolfloat

The tolerance, $ϵ$ used to check if the contrasts are orthogonal and if they are orthogonal to the mean. If $tol\le 0.0$ the value machine precision is used.

usetxbool

If $usetx=True$ the means ${\tau }_i^{\ast }$ are provided in $tx$ and the formula (2) is used instead of formula (1).

If $usetx=False$ formula (1) is used and $tx$ is not referenced.

txfloat, array-like, shape $\left(\text{nt}\right)$

If $usetx=True$ $tx$ must contain the means ${\tau }_{\text{i}}^{\ast }$, for $\text{i}=1,2,\dots ,t$.

Returns

estfloat, ndarray, shape $\left(\text{nc}\right)$

The estimates of the contrast, ${\Lambda }_{\text{j}}$, for $\text{j}=1,2,\dots ,\text{nc}$.

tablfloat, ndarray, shape $(\text{nc},5)$

The rows of the analysis of variance table for the contrasts. For each row column 1 contains the degrees of freedom, column 2 contains the sum of squares, column 3 contains the mean square, column 4 the $F$-statistic and column 5 the significance level for the contrast. Note that the degrees of freedom are always one and so the mean square equals the sum of squares.

Raises

NagValueError

(errno $1$)

On entry, $rdf=⟨value⟩$.

Constraint: $rdf\ge 1.0$.

(errno $1$)

On entry, $rms=⟨value⟩$.

Constraint: $rms>0.0$.

(errno $1$)

On entry, $\text{nt}=⟨value⟩$.

Constraint: $\text{nt}\ge 2$.

(errno $1$)

On entry, $\text{nc}=⟨value⟩$.

Constraint: $\text{nc}\ge 1$.

Warns

NagAlgorithmicWarning

(errno $2$)

The $⟨value⟩$$⟨value⟩$ and $⟨value⟩$$⟨value⟩$ contrasts are not orthogonal.

(errno $2$)

The $⟨value⟩$$⟨value⟩$ contrast is not orthogonal to the mean.

Notes

No equivalent traditional C interface for this routine exists in the NAG Library.

In the analysis of designed experiments the first stage is to compute the basic analysis of variance table, the estimate of the error variance (the residual or error mean square), ${\hat{\sigma }}^2$, and the (variance ratio) $F$-statistic for the $t$ treatments. If this $F$-test is significant then the second stage of the analysis is to explore which treatments are significantly different.

If there is a structure to the treatments then this may lead to hypotheses that can be defined before the analysis and tested using linear contrasts. For example, if the treatments were three different fixed temperatures, say $18$, $20$ and $22$, and an uncontrolled temperature (denoted by $N$) then the following contrasts might be of interest.

$$\begin{array}{c}\begin{array}{ccccc}& 18& 20& 22& N\\ & & & & \\ \left(a\right)& \frac{1}{3}& \frac{1}{3}& \frac{1}{3}& -1\\ & & & & \\ \left(b\right)& -1& 0& 1& 0\end{array}\end{array}$$

The first represents the average difference between the controlled temperatures and the uncontrolled temperature. The second represents the linear effect of an increasing fixed temperature.

For a randomized complete block design or a completely randomized design, let the treatment means be ${\hat{\tau }}_i$, $i=1,2,\dots ,t$, and let the $j$th contrast be defined by ${\lambda }_{ij}$, $i=1,2,\dots ,t$, then the estimate of the contrast is simply:

$${\Lambda }_j=\sum _{i=1}^t{\hat{\tau }}_i{\lambda }_{ij}$$

and the sum of squares for the contrast is:

$${SS}_j=\frac{{\Lambda }_j^2}{{\sum }_{i=1}^t{\lambda }_{ij}^2/n_i}$$

where $n_i$ is the number of observations for the $i$th treatment. Such a contrast has one degree of freedom so that the appropriate $F$-statistic is ${SS}_j/{\hat{\sigma }}^2$.

The two contrasts ${\lambda }_{ij}$ and ${\lambda }_{i{j}^{\prime }}$ are orthogonal if ${\sum }_{i=1}^t{\lambda }_{ij}{\lambda }_{i{j}^{\prime }}=0$ and the contrast ${\lambda }_{ij}$ is orthogonal to the overall mean if ${\sum }_{i=1}^t{\lambda }_{ij}=0$. In practice these sums will be tested against a small quantity, $ϵ$. If each of a set of contrasts is orthogonal to the mean and they are all mutually orthogonal then the contrasts provide a partition of the treatment sum of squares into independent components. Hence the resulting $F$-tests are independent.

If the treatments come from a design in which treatments are not orthogonal to blocks then the sum of squares for a contrast is given by:

$${SS}_j=\frac{{\Lambda }_j{\Lambda }_j^{\ast }}{{\sum }_{i=1}^t{\lambda }_{ij}^2/n_i}$$

where

$${\Lambda }_j^{\ast }=\sum _{i=1}^t{\tau }_i^{\ast }{\lambda }_{ij}$$

with ${\tau }_{\text{i}}^{\ast }$, for $\text{i}=1,2,\dots ,t$, being adjusted treatment means computed by first eliminating blocks then computing the treatment means from the block adjusted observations without taking into account the non-orthogonality between treatments and blocks. For further details see John (1987).

References

Cochran, W G and Cox, G M, 1957, Experimental Designs, Wiley

John, J A, 1987, Cyclic Designs, Chapman and Hall

Winer, B J, 1970, Statistical Principles in Experimental Design, McGraw–Hill