RCTグループ間の統計的検出力計算

install.packages("pwr", repos = "https://cran.rstudio.com/")

##  パッケージを 'C:/Users/miyamoto ayano/AppData/Local/R/win-library/4.3' 中にインストールします 
##  ('lib' が指定されていないため)

##  パッケージ 'pwr' は無事に展開され、MD5 サムもチェックされました 
## 
##  ダウンロードされたパッケージは、以下にあります 
##      C:\Users\miyamoto ayano\AppData\Local\Temp\Rtmp6P9Wt3\downloaded_packages

library(pwr)
# 比較する回数（統制群 vs 3つの介入群）
n_comparisons <- 3

# 補正後の有意水準
corrected_sig_level <- 0.05 / n_comparisons

# 補正後の有意水準を使って、必要なサンプルサイズを計算
power_result_bonferroni <- pwr.t.test(
  d = 0.2,         # 検出したい効果量 (小)
  sig.level = corrected_sig_level, # Bonferroni補正後の有意水準
  power = 0.8,       # 検出力 80%
  type = "two.sample",
  alternative = "two.sided"
)

print(power_result_bonferroni)

## 
##      Two-sample t test power calculation 
## 
##               n = 524.8907
##               d = 0.2
##       sig.level = 0.01666667
##           power = 0.8
##     alternative = two.sided
## 
## NOTE: n is number in *each* group