[R] formula for a mixed model ANOVA

Philippe Grosjean philippe.grosjean at ifremer.fr
Mon Sep 2 15:42:07 CEST 2002


Hi,

I am not very familiar with complex ANOVA, so, I am not sure about the
formula to use with aov() in this particular case. Could someone help me?

Here is my data. I have an unbalanced plan with 4 factors: "Site", "Sex",
"Age" and "Individual". "Site" is a 5-level fixed factor; "Sex" is of course
a 2-levels fixed factor, as is "Age". "Individual" is a nested random
factor: we sampled 4 to 10 individuals for each combinations Site*Sex*Age.
150 replicate measurements were made on each sampled individual ("OTM",
toxicity test; % of broken DNA for randomly selected cells).

if the formula "OTM ~ (Site*Sex*Age)/Individual" correct?

The purpose of this study is to choose a statistic that best summarizes the
150 measurements (mean, median, quantile, % of counts above a certain
level,...) in order to simplify the analysis to "stat(OTM) ~ Site*Sex*Age"
(where stat(OTM) is the statitstic applied by individual), in such a way
that intergroup variance is highest.

Any help or advise is welcome.
Thanks,

Philippe Grosjean


...........]<(({°<...............<°}))><...............................
( ( ( ( (
 ) ) ) ) )      Philippe Grosjean
( ( ( ( (
 ) ) ) ) )      IFREMER Nantes - DEL/AO
( ( ( ( (       rue de l'Ile d'Yeu, BP 21105, 44311 Nantes Cedex 3
 ) ) ) ) )      tel: (33) 02.40.37.42.29, fax: (33) 02.40.37.42.41
( ( ( ( (	e-mail: philippe.grosjean at ifremer.fr
 ) ) ) ) )
( ( ( ( (      "I'm 100% confident that p is between 0 and 1"
 ) ) ) ) )                                L. Gonick & W. Smith (1993)
.......................................................................



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