Estimación del tamaño de la muestra usando lenguaje R (2)

Este artículo tiene por defecto α = 0,05 (nivel de señal) y β = 0,2 (potencia), que se pueden ajustar según las necesidades de la investigación.

Paquete de importación

library(pwr)

(1) Desviación estándar conocida y diferencia esperada

1. Prueba t de una muestra: se espera que una determinada medida de tratamiento aumente el nivel de una determinada sustancia en 8 mg/l, con una desviación estándar de 10 mg/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. Prueba t de dos muestras: (1) Se sabe que la diferencia esperada entre el tratamiento A y el tratamiento B es del 30% de la desviación estándar y los tamaños de muestra de los dos grupos son iguales. Calcule el tamaño de la muestra (2) Se sabe que el tratamiento esperado A es peor que el tratamiento B. La desviación estándar es 60% y se sabe que 90 personas reciben el tratamiento A. Calcule el tamaño de la muestra de 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) Tasa de muestra conocida: (1) Se sabe que la tasa de complicaciones esperada del tratamiento A es del 7% y la tasa de complicaciones del tratamiento general es del 12%. Calcule el tamaño de la muestra de A. (2) Se sabe que la tasa de complicaciones esperada del tratamiento A es del 7% y la del tratamiento B es del 12%. Calcule los tamaños de muestra de A y B. (3) Se sabe que la tasa de complicaciones esperada del tratamiento A es del 7%, el tamaño de la muestra es 164 y la tasa de complicaciones del tratamiento B es del 12%. Calcule el tamaño de la muestra de B. (4) Se sabe que la efectividad esperada del método A es del 90%, el método B es del 70%, el método C es del 80% y el método D es del 60%; calcule el tamaño total de la muestra.

#（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
)