title: “功效分析-方差分析”
output: html_document
此为宫廷土豆豆《R语言实战》学习笔记,感谢作者Kabacoff以及译者高涛等人。——2019.9.27
**pwr.anova.test()**可以对平衡单因素方差分析进行统计功效的计算。格式为:
pwr.anova.test(k=, n=, f=, sig.level=, power= )
,其中k为组的个数,n是各组中的样本大小。
对于单因素方差分析,效应值可以通过f来衡量:
现在举例说明函数的用法。现对5个组做单因素方差分析,要达到0.8的功效,效应值为0.25(如果面对全新的研究情况,没有任何过去的经验可借鉴时, 可以选择Cohen提出的基准值,中等的效应值是0.25)。选择0.05的显著性水平,计算各组需要的样本大小。代码如下:
library(pwr)
pwr.anova.test(k=5, f=0.25, sig.level=0.05, power=0.8)
结果为:
Balanced one-way analysis of variance pow