[R] MANOVA / Hotelling's Test

Peter Dalgaard BSA p.dalgaard at biostat.ku.dk
Wed Nov 11 15:58:14 CET 1998


"Uli Flenker; Raum 704" <uli at biochem.dshs-koeln.de> writes:

> I have to perform Hotteling's T^2 Test - more generally a MANOVA - on a
> set of data.  Is there a (simple?) possibility to do it in R? 
> 
> 
> The somewhat obvious way to do it would be
> 
> > summary(aov(Y ~ x1 + x2 ... + xn)) 
> 
> where Y would be a two-column matrix. 
> 
> 
> But this does not work in any possible combination of matrix or factor
> dimensions! Is it principally not possible to perform modelling of
> multivariate data using the "~"-operator? 
> 
> I'm using R 0.62.2 on Linux 2.0.25. 

Let's put it this way: It isn't, but it should be...

0.63 will contain some of the pieces, but not all. There's an extended
aov() and objects of class "maov", and in general you can use
matrix-responses. However, all you get from summary and friends are
lists of univariate results. Splus has manova(), which is essentially
just a wrapper for aov() which attaches class "manova" to the result,
but also summary.manova() computing the relevant test statistics. R
doesn't have this yet.

Probably the most practical way to proceed is to calculate generalised
sums of squares and products with 

> r<-resid(lm(y~x))
> t(r)%*%r

and then put TW Anderson's book in your lap and calculate Wilk's
lambda and asymptotic approximations of it's distibution. (Hotelling's
T^2 is some trivial transformation of W's lambda and can be converted
to an exact F distributed variable)

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list