[R] manova with Error?
Peter Dalgaard BSA
p.dalgaard at biostat.ku.dk
Thu Oct 17 15:54:55 CEST 2002
"RenE J.V. Bertin" <rjvbertin at hotmail.com> writes:
> On 17 Oct 2002 13:50:30 +0200, Peter Dalgaard BSA <p.dalgaard at biostat.ku.dk> wrote regarding "Re:
> [R] manova with Error?"
> 8-) This is not my specialty, but could it be that you have an empty error
> 8-) stratum? I.e. if Subject:fac1:fac2 is the entire data set, then there
> 8-) is no DF for the residual error. (If so, try Error(Subject/(fac1+fac2)) )
> I'm not sure what you mean with "the entire data set". I have one value for obs1 and one for obs2 for each combination of Subject, fac1 and fac2.
That's what I meant.
> I tried your suggestion; got the same error message.
> BTW: Would that not also cause problems in a regular anova, with aov?
Not necessarily. (I knew that it works there.)
I think what you're trying to do makes sense, it's just that
summary.manova isn't quite up to the task (or perhaps aov() or
To let others chip in, let's try on data that are available to
example(summary.manova) # includes defn. of Y, rate, additive
fit <- manova(Y ~ rate + Error(additive))
of which only the raw print of "fit" and the very last line give
I would expect at least summary(fit[]) to work, but it is stumbling
ss <- list(0)
nmrows <- character(0)
and a little later
names(ss) <- nmrows
Anyone have further clues?
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