Re: [R] Unbalanced Anova: What is the best approach?

From: peter dalgaard <pdalgd_at_gmail.com>
Date: Sun, 03 Apr 2011 11:10:34 +0200

On Apr 3, 2011, at 09:24 , Krishna Kirti Das wrote:

> I have a three-way unbalanced ANOVA that I need to calculate (fixed effects
> plus interactions, no random effects). But word has it that aov() is good
> only for balanced designs. I have seen a number of different recommendations
> for working with unbalanced designs, but they seem to differ widely (car,
> nlme, lme4, etc.). So I would like to know what is the best or most usual
> way to go about working with unbalanced designs and extracting a reliable
> ANOVA table from them in R?

Actually, without random effects, aov() is not too crazy, but you might as well use plain lm(). In both cases, the main point is that you need to be aware that there is no such thing as "the" ANOVA table: Sums of squares will depend on the order of testing, and there is nothing to do about that (except getting balanced data).

Pragmatically, I'd test the three-factor interaction, then use drop1() on a model with two-factor interactions, if nothing glaringly obvious pops up, try reduction to additive model and then use drop1() again. Obviously, if significant interactions appear, you cannot just remove them and need to investigate what they mean.

-- 
Peter Dalgaard
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Email: pd.mes_at_cbs.dk  Priv: PDalgd_at_gmail.com

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Sun 03 Apr 2011 - 09:14:36 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sun 03 Apr 2011 - 15:00:26 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive