If Month is nested within Quadrat I think you want: aov(ProportioninTreatment ~ Treatment*Month +Error(Quadrat/Month), RM) If Treatment is also nested within Quadrat, you want: aov(ProportioninTreatment ~ Treatment*Month +Error(Quadrat/(Treatment*Month)), RM) On Wed, Apr 1, 2009 at 12:42 AM, Jessica L Hite/hitejl/O/VCU
<hitejl at vcu.edu> wrote:
I have data on the proportion of clutches experiencing different fates (e.g., 4 different sources of mortality) for 5 months . I need to test 1) if the overall proportion of these different fates is different over the entire study and 2) to see if there are monthly differences within (and among) fate types. Thus, I am pretty sure this is an RM analysis -( I measure the same quadrats each month). I am fine running the analysis in R - with the code below, however, there is no output for the among group variation...this is an important component - any ideas on how to solve this problem? I have included code and sample data below. Many thanks in advance for help and suggestions. J ?both.aov <- aov(ProportioninTreatment ~ factor(Treatment)*factor(Month) + Error(factor(Quadrat)), RM) Error: factor(id) ? ? ? ? ?Df ?Sum Sq Mean Sq F value Pr(>F) Residuals ?3 0.51619 0.17206 ? ? ? ? ? ? ? #####why only partial output here? ####### Error: Within ? ? ? ? ? ? ? ? ? Df Sum Sq Mean Sq F value ? Pr(>F) factor(Fate1) ? ? ? 3 1.2453 ?0.4151 ?3.5899 0.017907 * time ? ? ? ? ? ? ? ?1 0.9324 ?0.9324 ?8.0637 0.005929 ** factor(Fate1):time ?3 0.9978 ?0.3326 ?2.8763 0.042272 * Residuals ? ? ? ? ?69 7.9783 ?0.1156 Fate1 Proportion in Fate ? ? ?ASIN ?Month Quadrat 1 ? ? 0.117647059 0.350105778 1 ? ? 1 1 ? ? 0 ? ? 0 ? ? 2 ? ? 1 1 ? ? 0.111111111 0.339836909 3 ? ? 1 1 ? ? 0 ? ? 0 ? ? 4 ? ? 1 1 ? ? 0 ? ? 0 ? ? 5 ? ? 1 1 ? ? 0 ? ? 0 ? ? 1 ? ? 2 1 ? ? 0 ? ? 0 ? ? 2 ? ? 2 1 ? ? 0.2 ? 0.463647609 3 ? ? 2 1 ? ? 0.25 ?0.523598776 4 ? ? 2 1 ? ? 0.111111111 0.339836909 5 ? ? 2 1 ? ? 0 ? ? 0 ? ? 1 ? ? 3 1 ? ? 0 ? ? 0 ? ? 2 ? ? 3 1 ? ? 0 ? ? 0 ? ? 3 ? ? 3 1 ? ? 0.384615385 0.668964075 4 ? ? 3 1 ? ? 0 ? ? 0 ? ? 5 ? ? 3 1 ? ? 0 ? ? 0 ? ? 1 ? ? 4 1 ? ? 0 ? ? 0 ? ? 2 ? ? 4 1 ? ? 0 ? ? 0 ? ? 3 ? ? 4 1 ? ? 0.166666667 0.420534336 4 ? ? 4 1 ? ? 0 ? ? 0 ? ? 5 ? ? 4 2 ? ? 0.352941176 0.636132062 1 ? ? 1 2 ? ? 0.2 ? 0.463647609 2 ? ? 1 2 ? ? 0.333333333 0.615479708 3 ? ? 1 2 ? ? 1 ? ? 1.570796327 4 ? ? 1 2 ? ? 0 ? ? 0 ? ? 5 ? ? 1 2 ? ? 0.5 ? 0.785398163 1 ? ? 2 2 ? ? 0 ? ? 0 ? ? 2 ? ? 2 2 ? ? 0.6 ? 0.886077124 3 ? ? 2 2 ? ? 0.416666667 0.701674124 4 ? ? 2 2 ? ? 0.222222222 0.490882678 5 ? ? 2 2 ? ? 0 ? ? 0 ? ? 1 ? ? 3 2 ? ? 0.2 ? 0.463647609 2 ? ? 3 2 ? ? 0 ? ? 0 ? ? 3 ? ? 3 2 ? ? 0.461538462 0.746898594 4 ? ? 3 2 ? ? 0 ? ? 0 ? ? 5 ? ? 3 2 ? ? 0 ? ? 0 ? ? 1 ? ? 4 2 ? ? 0 ? ? 0 ? ? 2 ? ? 4 2 ? ? 0.307692308 0.588002604 3 ? ? 4 2 ? ? 0.666666667 0.955316618 4 ? ? 4 2 ? ? 0 ? ? 0 ? ? 5 ? ? 4 3 ? ? 0 ? ? 0 ? ? 1 ? ? 1 3 ? ? 0 ? ? 0 ? ? 2 ? ? 1 3 ? ? 0.444444444 0.729727656 3 ? ? 1 3 ? ? 0 ? ? 0 ? ? 4 ? ? 1 3 ? ? 1 ? ? 1.570796327 5 ? ? 1 3 ? ? 0.5 ? 0.785398163 1 ? ? 2 3 ? ? 0 ? ? 0 ? ? 2 ? ? 2 3 ? ? 0 ? ? 0 ? ? 3 ? ? 2 3 ? ? 0.25 ?0.523598776 4 ? ? 2 3 ? ? 0.555555556 0.841068671 5 ? ? 2 3 ? ? 0 ? ? 0 ? ? 1 ? ? 3 3 ? ? 0 ? ? 0 ? ? 2 ? ? 3 3 ? ? 0 ? ? 0 ? ? 3 ? ? 3 3 ? ? 0.153846154 0.403057075 4 ? ? 3 3 ? ? 0.666666667 0.955316618 5 ? ? 3 3 ? ? 0 ? ? 0 ? ? 1 ? ? 4 3 ? ? 0 ? ? 0 ? ? 2 ? ? 4 3 ? ? 0 ? ? 0 ? ? 3 ? ? 4 3 ? ? 0 ? ? 0 ? ? 4 ? ? 4 3 ? ? 0.875 1.209429203 5 ? ? 4 4 ? ? 0.294117647 0.573203309 1 ? ? 1 4 ? ? 0.2 ? 0.463647609 2 ? ? 1 4 ? ? 0 ? ? 0 ? ? 3 ? ? 1 4 ? ? 0 ? ? 0 ? ? 4 ? ? 1 4 ? ? 0 ? ? 0 ? ? 5 ? ? 1 4 ? ? 0 ? ? 0 ? ? 1 ? ? 2 4 ? ? 0 ? ? 0 ? ? 2 ? ? 2 4 ? ? 0 ? ? 0 ? ? 3 ? ? 2 4 ? ? 0.083333333 0.292842771 4 ? ? 2 4 ? ? 0.111111111 0.339836909 5 ? ? 2 4 ? ? 0 ? ? 0 ? ? 1 ? ? 3 4 ? ? 0 ? ? 0 ? ? 2 ? ? 3 4 ? ? 0 ? ? 0 ? ? 3 ? ? 3 4 ? ? 0 ? ? 0 ? ? 4 ? ? 3 4 ? ? 0.166666667 0.420534336 5 ? ? 3 4 ? ? 0 ? ? 0 ? ? 1 ? ? 4 4 ? ? 0 ? ? 0 ? ? 2 ? ? 4 4 ? ? 0.461538462 0.746898594 3 ? ? 4 4 ? ? 0 ? ? 0 ? ? 4 ? ? 4 4 ? ? 0.125 0.361367124 5 ? ? 4 ? ? ? ?[[alternative HTML version deleted]]
______________________________________________ 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.
Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar: http://tinyurl.com/mikes-public-calendar ~ Certainty is folly... I think. ~