[R] manova with Error?

[Peter Dalgaard BSA]

"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?"
> 
> Hi,
> 
> 8-) 
> 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.

Too bad.
 
> 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
manova() isn't).

To let others chip in, let's try on data that are available to
everyone:

example(summary.manova) # includes defn. of Y, rate, additive
fit <- manova(Y ~ rate + Error(additive))
fit
summary(fit)
summary.aov(fit)
summary(fit[[1]])
summary(fit[[2]])
summary(fit[[3]])

of which only the raw print of "fit" and the very last line give
sensible results. 

I would expect at least summary(fit[[2]]) to work, but it is stumbling
upon

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