7 統計解析

制作中

inferパッケージのchisq_test()が便利。
tidy dataから、一発で検定できる。

library(infer)
mtcars %>%
  mutate(cyl = factor(cyl), am = factor(am)) %>%
  chisq_test(cyl ~ am)

7.1 検定

7.1.1 t検定

7.1.2 Wilcoxonの順位和検定

{coin}のwilcox_test()を使おう。
参考文献

wilcox_test(number ~ group, data = df, distribution = "exact")

7.1.3 分割表の検定

• カイ二乗検定, フィッシャーの正確検定

7.2 回帰分析

7.2.1 重回帰

7.2.2 ロジスティック回帰

7.2.2.1 broom

library(broom)
tidy(model, conf.int = TRUE, exponentiate = TRUE)

7.3 生存時間分析

• https://www.emilyzabor.com/tutorials/survival_analysis_in_r_tutorial.html

7.3.1 Cox回帰

{survival}と{survminer}を使う。

library(survival)
library(survminer)
cox <- coxph(Surv(time, event) ~ group + sex + age, method = "breslow", data = dat)
ggforest(cox)