From: Sumitrajit Dhar <s-dhar_at_northwestern.edu>

Date: Sat, 29 Nov 2008 21:40:54 -0600

[...]

Date: Sat, 29 Nov 2008 21:40:54 -0600

I am trying to figure out how run a repeated measures ANOVA on the following data set.

1 FSI052A A 1 NA 2 FSI052B B 1 NA 3 FSI053A A 1 NA 4 FSI055A A 1 NA 5 FSI055B B 1 NA 6 FSI057A A 1 NA 7 FSI057B B 1 NA 8 FSI058A A 1 NA 9 FSI058B B 1 12.93 10 FSI060A A 1 NA 11 FSI060B B 1 NA 12 FSI067A A 1 NA 13 FSI070A A 1 NA 14 FSI070B B 1 NA 15 FSI075A A 1 NA 16 FSI075B B 1 NA 17 FSI078A A 1 NA 18 FSI078B B 1 NA 19 FSI052A A 2 13.65 20 FSI052B B 2 11.71 21 FSI053A A 2 20.91 22 FSI055A A 2 9.52 23 FSI055B B 2 13.43 24 FSI057A A 2 14.57 25 FSI057B B 2 10.00 26 FSI058A A 2 14.95 27 FSI058B B 2 16.97 28 FSI060A A 2 13.06 29 FSI060B B 2 7.96

[...]

Both "frequency" and "trial" are factorized.

I know that I should be using 'lme' rather than aov, but....

> summary(aov(dplvl~trial*frequency+Error(subject/

(trial*frequency)),data=tenB))

Df Sum Sq Mean Sq F value Pr(>F) trial 1 7.15 7.15 0.2566 0.6203 frequency 1 40.84 40.84 1.4657 0.2461 trial:frequency 1 0.03 0.03 0.0013 0.9723 Residuals 14 390.11 27.86 Error: subject:frequency Df Sum Sq Mean Sq F value Pr(>F) frequency 1 183.893 183.893 11.925 0.003271 ** trial:frequency 1 16.393 16.393 1.063 0.317866 Residuals 16 246.738 15.421

--- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Error: Within Df Sum Sq Mean Sq F value Pr(>F) Residuals 35 269.13 7.69 What I do not understand is why I am getting df=1 for frequency when there are 5 levels of that factor. Thanks for all your help. Regards, Sumit [[alternative HTML version deleted]]Received on Sun 30 Nov 2008 - 03:46:56 GMT______________________________________________ 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.

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 30 Nov 2008 - 11:30:28 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.
*