[R] MANOVA extension of paired t-test?

[Stewart Schultz]

I would like to test the hypothesis that the difference between pairs, for
several variables, is zero.  This is easily done separately for each
variable with:

lm(Y ~ rep(0, nrow(Y)))

where Y is a matrix whose columns are the differences for each variable
between pair members.

However, I would like to get an overall probability across all variables
from a Wilks or Pillai-Bartlett statistic as in MANOVA.  However, the
analogous command:

manova(Y ~ rep(0, nrow(Y)))

does not work presumably because the predictor variable is not a factor.

This would be a repeated-measures MANOVA, except that it is a trivial case
of only one repeat, so the difference between pair members can simply be
used as the variable of interest.


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._