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 <hit...@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@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. ~ ______________________________________________ R-help@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.
