幹細胞治療の有効性を評価する(続き)
対照群と治療群のヒツジの心臓における平均の心ポンプ機能に差があるかを評価するため、仮説検定を実施し、p値を計算します。
差を検定するので、差をとる順番が重要です。calculate() 関数の order 引数で指定できます。たとえば、group1 - group2 を得るには order = c("group1", "group2") を使います。
この演習はコースの一部です
R による数値データの推測
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
n_replicates <- 1000
# Generate 1000 differences in means via randomization
diff_mean_ht <- stem.cell %>%
# y ~ x
specify(___ ~ ___) %>%
# Null = no difference between means
hypothesize(null = "___") %>%
# Shuffle labels 1000 times
generate(reps = ___, type = "___") %>%
# Calculate test statistic
calculate(stat = "___", order = c("esc", "___")