A/B测试之t-test

贝叶斯的t检验可以参考John K. Kruschke的著作Doing Bayesian Data Analysis的第16章。

本文的例子使用的数据集可以从github取得：

Doing-Bayesian-Data-Analysis-in-brms-and-the-tidyverse/TwoGroupIQ.csv at master · ASKurz/Doing-Bayesian-Data-Analysis-in-brms-and-the-tidyverse · GitHub

这是关于某种药物能否影响IQ的实验。

首先，导入数据并拟合模型

library(tidyverse)
library(tidybayes)
library(brms)

my_data <- read_csv("TwoGroupIQ.csv")
(mean_y <- mean(my_data$Score))
(sd_y <- sd(my_data$Score))
stanvars <- 
  stanvar(mean_y, name = "mean_y") + 
  stanvar(sd_y,   name = "sd_y") + 
  stanvar(1/29,   name = "one_over_twentynine")
fi