[R] R Anova Analysis
Yanwei Tan
Tan at nbio.uni-heidelberg.de
Thu Apr 29 13:52:12 CEST 2010
Dear all,
I have a quite basic questions about anova analysis in R, sorry for
this, but I have no clue how to explain this result.
I have two datasets which are named: nmda123, nmda456. Each dataset has
three samples which were measured three times. And I would like to
compare means of them with Posthoc test using R, following please see
the output:
(CREB, mCREB and No virus are the name of samples)
> nmda123
Values ind
1 6.7171265 CREB
2 5.0343117 CREB
3 6.9000000 CREB
4 0.1195394 mCREB
5 0.1221876 mCREB
6 0.1900000 mCREB
7 1.0000000 No Virus
8 1.0000000 No Virus
9 1.0000000 No Virus
> nmda456
Values ind
1 6.4486940 CREB
2 6.2277490 CREB
3 6.5000000 CREB
4 0.2000000 mCREB
5 0.3766052 mCREB
6 0.4000000 mCREB
7 1.0000000 No Virus
8 1.0000000 No Virus
9 1.0000000 No Virus
> TukeyHSD(aov(Values ~ ind, data = nmda456))
Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = Values ~ ind, data = nmda456)
$ind
diff lwr upr p adj
mCREB-CREB -6.0666126 -6.3289033 -5.8043219 0.0000000
No Virus-CREB -5.3921477 -5.6544383 -5.1298570 0.0000000
No Virus-mCREB 0.6744649 0.4121743 0.9367556 0.0005382
> TukeyHSD(aov(Values ~ ind, data = nmda123))
Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = Values ~ ind, data = nmda123)
$ind
diff lwr upr p adj
mCREB-CREB -6.073237 -7.5618886 -4.584585 0.0000392
No Virus-CREB -5.217146 -6.7057976 -3.728495 0.0000943
No Virus-mCREB 0.856091 -0.6325606 2.344743 0.2588450
So my question is No virus-mCREB group. Even I looked at the data by
eyes, there is big difference between no virus and mCREB in data
nmda123, but why the pvalue is 0.2588450, but in nmda456 data, the
pvalue is 0.0005382. But I can see there is bigger difference in
nmda123 than nmda456, I do not know why. Sorry for my inexperiences in
statistics.
Thanks for your reply and time!
Cheers,
Wei
More information about the R-help
mailing list