naginterfaces.library.anova.confidence

naginterfaces.library.anova.confidence(typ, tmean, rdf, c, clevel)

confidence computes simultaneous confidence intervals for the differences between means. It is intended for use after random() or rowcol().

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

https://support.nag.com/numeric/nl/nagdoc_30.2/flhtml/g04/g04dbf.html

Parameters

typstr, length 1

Indicates which method is to be used.

$typ=\text{'T'}$

The Tukey–Kramer method is used.

$typ=\text{'B'}$

The Bonferroni method is used.

$typ=\text{'D'}$

The Dunn–Sidak method is used.

$typ=\text{'L'}$

The Fisher LSD method is used.

$typ=\text{'S'}$

The Scheffe’s method is used.

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

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

rdffloat

$\nu$, the residual degrees of freedom.

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

The strictly lower triangular part of $c$ must contain the standard errors of the differences between the means as returned by random() and rowcol(). That is $c[i,j]$, $i>j$, contains the standard error of the difference between the $i$th and $j$th mean in $tmean$.

clevelfloat

The required confidence level for the computed intervals, ($1-\alpha$).

Returns

cilfloat, ndarray, shape $(\text{nt}\times (\text{nt}-1)/2)$

The $((\text{i}-1)(\text{i}-2)/2+\text{j})$th element contains the lower limit to the confidence interval for the difference between $\text{i}$th and $\text{j}$th means in $tmean$, for $\text{j}=1,2,\dots ,\text{i}-1$, for $\text{i}=2,3,\dots ,t$.

ciufloat, ndarray, shape $(\text{nt}\times (\text{nt}-1)/2)$

The $((\text{i}-1)(\text{i}-2)/2+\text{j})$th element contains the upper limit to the confidence interval for the difference between $\text{i}$th and $\text{j}$th means in $tmean$, for $\text{j}=1,2,\dots ,\text{i}-1$, for $\text{i}=2,3,\dots ,t$.

isigint, ndarray, shape $(\text{nt}\times (\text{nt}-1)/2)$

The $((\text{i}-1)(\text{i}-2)/2+\text{j})$th element indicates if the difference between $\text{i}$th and $\text{j}$th means in $tmean$ is significant, for $\text{j}=1,2,\dots ,\text{i}-1$, for $\text{i}=2,3,\dots ,t$. If the difference is significant then the returned value is $1$; otherwise the returned value is $0$.

Raises

NagValueError

(errno $1$)

On entry, $typ=⟨value⟩$.

Constraint: $typ=\text{'T'}$, $\text{'B'}$, $\text{'D'}$, $\text{'L'}$ or $\text{'S'}$.

(errno $1$)

On entry, $clevel=⟨value⟩$.

Constraint: $0.0<clevel<1.0$.

(errno $1$)

On entry, $rdf=⟨value⟩$.

Constraint: $rdf\ge 1.0$.

(errno $1$)

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

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

(errno $2$)

On entry, $i=⟨value⟩$ and $j=⟨value⟩$.

Constraint: $c[i-1,j-1]>0.0$.

(errno $3$)

There has been a failure in the computation of the studentized range statistic.

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 the computation of analysis of a designed experiment 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$, the residual degrees of freedom, $\nu$, and the (variance ratio) $F$-statistic for the $t$ treatments. The second stage of the analysis is to compare the treatment means. If the treatments have no structure, for example the treatments are different varieties, rather than being structured, for example a set of different temperatures, then a multiple comparison procedure can be used.

A multiple comparison procedure looks at all possible pairs of means and either computes confidence intervals for the difference in means or performs a suitable test on the difference. If there are $t$ treatments then there are $t(t-1)/2$ comparisons to be considered. In tests the type $1$ error or significance level is the probability that the result is considered to be significant when there is no difference in the means. If the usual $t$-test is used with, say, a $6\%$ significance level then the type $1$ error for all $k=t(t-1)/2$ tests will be much higher. If the tests were independent then if each test is carried out at the $100\alpha$ percent level then the overall type $1$ error would be ${\alpha }^{\ast }=1-{(1-\alpha )}^k\simeq k\alpha$. In order to provide an overall protection the individual tests, or confidence intervals, would have to be carried out at a value of $\alpha$ such that ${\alpha }^{\ast }$ is the required significance level, e.g., five percent.

The $100(1-\alpha )$ percent confidence interval for the difference in two treatment means, ${\hat{\tau }}_i$ and ${\hat{\tau }}_j$ is given by

$$({\hat{\tau }}_i-{\hat{\tau }}_j)\pm T_{(\alpha ,\nu ,t)}^{\ast }se({\hat{\tau }}_i-{\hat{\tau }}_j)\text{,}$$

where $se\left(\right)$ denotes the standard error of the difference in means and $T_{(\alpha ,\nu ,t)}^{\ast }$ is an appropriate percentage point from a distribution. There are several possible choices for $T_{(\alpha ,\nu ,t)}^{\ast }$. These are:

1. $\frac{1}{2}{q}_{(1-\alpha ,\nu ,t)}$, the studentized range statistic, see stat.inv_cdf_studentized_range. It is the appropriate statistic to compare the largest mean with the smallest mean. This is known as Tukey–Kramer method.

2. $t_{(\alpha /k,\nu )}$, this is the Bonferroni method.

3. $t_{({\alpha }_0,\nu )}$, where ${\alpha }_0=1-{(1-\alpha )}^{1/k}$, this is known as the Dunn–Sidak method.

4. $t_{(\alpha ,\nu )}$, this is known as Fisher’s LSD (least significant difference) method. It should only be used if the overall $F$-test is significant, the number of treatment comparisons is small and were planned before the analysis.

5. $\sqrt{(k-1)F_{1-\alpha ,k-1,\nu }}$ where $F_{1-\alpha ,k-1,\nu }$ is the deviate corresponding to a lower tail probability of $1-\alpha$ from an $F$-distribution with $k-1$ and $\nu$ degrees of freedom. This is Scheffe’s method.

In cases (b), (c) and (d), $t_{(\alpha ,\nu )}$ denotes the $\alpha$ two tail significance level for the Student’s $t$-distribution with $\nu$ degrees of freedom, see stat.inv_cdf_students_t.

The Scheffe method is the most conservative, followed closely by the Dunn–Sidak and Tukey–Kramer methods.

To compute a test for the difference between two means the statistic,

$$\frac{{\hat{\tau }}_i-{\hat{\tau }}_j}{se({\hat{\tau }}_i-{\hat{\tau }}_j)}$$

is compared with the appropriate value of $T_{(\alpha ,\nu ,t)}^{\ast }$.

References

Kotz, S and Johnson, N L (ed.), 1985, Multiple range and associated test procedures, Encyclopedia of Statistical Sciences (5), Wiley, New York

Kotz, S and Johnson, N L (ed.), 1985, Multiple comparison, Encyclopedia of Statistical Sciences (5), Wiley, New York

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