功效分析-代码笔记

RStudio面板展示

配套代码

#功效分析 pwr包

install.packages("pwr")
library(pwr)

pwr.f2.test(u=3,f2=0.0769,sig.level=0.05,power=0.90)    
#sig.leve为显著性水平，power为功效，u为分母自由度，v为分子自由度，f2为效应值
pwr.f2.test(u=3,sig.level = 0.05,power=0.9,f2=0.0769)
#假设显著性水平为0.05，在百分之九十的置信度情况下，至少需要185个受试者才可以

#方差功效分析案案例 两组样品做单因素方差分析
pwr.anova.test(k=2,f=0.25,sig.level = 0.05,power=0.9)
#功效fower为0.9，f效应值为0.25，k为组的个数，n为组的大小为要求解的量，
#最终求得n为85.03，每组需要86个样本，总样本数需要172个