rm.hotel: Repeated measures ANOVA (univariate data) using Hotelling's...

In mvhtests: Multivariate Hypothesis Tests

Repeated measures ANOVA (univariate data) using Hotelling's T^2 test

Description

Repeated measures ANOVA (univariate data) using Hotelling's T^2 test.

Usage

rm.hotel(x, a = 0.05)

Arguments

x	A numerical matrix with the repeated measurements. Each column contains the values of the repeated measurements.
a	The level of significance, default value is equal to 0.05.

Details

We now show how can one use Hotelling's T^2 test to analyse univariate repeated measures. Univariate analysis of variance for repeated measures is the classical way, but we can use this multivariate test as well. In the repeated measures ANOVA case, we have many repeated observations from the same n subjects, usually at different time points and the interest is to see whether the means of the samples are equal or not \mu_1=\mu_2=\ldots=\mu_k assuming k repeated measurements. We can of course change this null hypothesis and test many combinations of means. The idea in any case is to construct a matrix of contrasts. I will focus here in the first case only and in particular the null hypothesis and the matrix of contrasts \bf C are

\left( {\begin{array}{c} \mu_1=\mu_2 \\ \mu_2=\mu_3 \\ \vdots \\ \mu_{k-1}=\mu_k \end{array}} \right)= \left( {\begin{array}{ccccc} 1 & -1 & 0 & \ldots & 0 \\ 1 & 0 & -1 & \dots & 0 \\ \vdots & \vdots & \vdots & \vdots & \vdots \\ 1 & 0 & 0 & \ldots & -1 \\ \end{array}} \right)\pmb{\mu}={\bf C}\pmb{\mu}.

The contrast matrix \bf C has k-1 independent rows and if there is no treatment effect, {\bf C}\pmb{\mu}={\bf 0}.

The test statistic is

T_r^2=\frac{\left(n-k+1\right)}{\left(n-1\right)\left(k-1\right)}n\left({\bf C}\bar{\bf x}\right)^T \left({\bf CSC}^T\right)^{-1}\left({\bf C}\bar{\bf x}\right) \sim F_{k-1,n-k+1}.

Value

A list including:

m	The mean vector.
result	A vector with the test statistic value, it's associated p-value, the numerator and denominator degrees of freedom and the critical value.

Author(s)

Michail Tsagris.

R implementation and documentation: Michail Tsagris [e-mail].

See Also

maov, hotel2T2, james

Examples

x <- as.matrix(iris[, 1:4]) ## assume they are repeated measurements
rm.hotel(x)

mvhtests documentation built on April 4, 2025, 5:07 a.m.