R语言实现样本量的估算(2)

本文默认α=0.05(sig.level),β=0.2(power),根据研究需要可调整。

导入包

复制代码
library(pwr)

(1)已知标准差和预期差异

1、单样本t检验:某治疗措施预期提高某物质水平8mg/L,标准差为10mg/L。

复制代码
#单样本t检验
pwr.t.test(d = 8/10, #期望的平均差异/总体标准差
           sig.level = 0.05,
           power = 1-0.2,
           type = "one.sample",
           alternative = "greater")#two.sided两边;greater:表示单侧检验中的大尾检验;less: 表示单侧检验中的小尾检验

2、两样本t检验:(1)已知预期A治疗措施与B治疗措施相差30%的标准差,两组样本量相同,计算样本量(2)已知预期A治疗措施差于B治疗措施60%的标准差,已知A治疗措施90人,计算B样本量。

复制代码
#(1)两样本t检验(样本数量相同)
pwr.t.test(d = 0.3,#标准差的30%
           sig.level = 0.05,
           power = 1 - 0.2,
           type = "two.sample",
           alternative = "two.sided"
)
#(2)两样本t检验(样本数量不同)
pwr.t2n.test(d=0.6,
             n1=90,
             sig.level = 0.05, 
             power =0.8,
             alternative="greater")

(2)已知样本率:(1)已知A治疗措施预期并发症率7%,一般治疗措施并发症率12%,计算A的样本量。(2)已知A治疗措施预期并发症率7%,B治疗措施并发症率12%,计算A、B的样本量。(3)已知A治疗措施预期并发症率7%,样本量164,B治疗措施并发症率12%,计算B的样本量。(4)已知预期A法有效率为90%,B法70%,C法80%,D法60%,计算总的样本量。

复制代码
#(1)一组样本
pwr.p.test(h=ES.h(0.07,0.12),power=0.8,sig.level=0.05,alternative="less")
#(2)两组相同样本
pwr.2p.test(h=ES.h(0.07,0.12),power=0.8,sig.level=0.05,alternative="less")
#(3)两组不同样本量
pwr.2p2n.test(h = ES.h(0.07,0.12),
              n1=164,
              n2=90,
              sig.level = 0.05,
              alternative = "less"
)
#(4)多样本率
prob <- rbind(c(0.9, 0.7, 0.8,0.6), # 有效率
              c(0.1, 0.3, 0.2,0.4)) # 无效率
pwr.chisq.test(w = ES.w2(prob/4), # 效应大小
               df = (4-1)*(3-1), #自由度
               sig.level = 0.05,
               power = 1-0.2
)
相关推荐
烟锁池塘柳01 小时前
【R语言】R 语言中打印含有双引号的字符串时会出现 “\” 的原因解析
r语言
全栈开发圈4 天前
干货分享|如何从0到1掌握R语言数据分析
开发语言·数据分析·r语言
小杜的生信筆記7 天前
基于R语言,“上百种机器学习模型”学习教程 | Mime包
开发语言·学习·机器学习·r语言·sci
在打豆豆的小潘学长7 天前
【R语言】多样本单细胞分析_SCTransform+Harmony方案(2)
开发语言·r语言
TS的美梦7 天前
ROGUE: 【张院士团队R包】一种基于熵的用于评估单细胞群体纯度的度量标准
开发语言·r语言
weixin_493202639 天前
R语言代码加密(1)
r语言
Tiger Z9 天前
《R for Data Science (2e)》免费中文翻译 (第3章) --- Data transformation(2)
r语言·数据科学·中文翻译
星石传说11 天前
使用R将nc文件转换为asc文件或者tif文件
r语言·生信
Mister Leon12 天前
机器学习Adaboost算法----SAMME算法和SAMME.R算法
算法·机器学习·r语言
Tiger Z13 天前
R 语言科研绘图第 67 期 --- 箱线图-显著性
r语言·论文·科研·绘图·研究生