R 中的 A/B 测试
Lauryn Burleigh
Data Scientist
library(ggplot2)
ggplot(pizza, aes(x = Time,
fill = Topping)) +
geom_histogram() +
facet_grid(Topping~.)

library(pwr)
pwr.t.test(d = 0.73, power = 0.80,
sig.level = 0.05,
type = "two.sample",
alternative = "two.sided")
双样本 t 检验功效计算
n = 30.44799
d = 0.73
sig.level = 0.05
power = 0.8
alternative = two.sided
注意:n 为每组样本量
不显著(p > 0.05)= 方差相等
若显著(p < 0.05)= 方差不等
library(car)
leveneTest(Time ~ Topping,
data = Pizza)
Levene 方差齐性检验
Df F value Pr(>F)
group 1 0.1457 0.7031
t.test(Time ~ Topping, data = Pizza,
paired = FALSE,
alternative = "two.sided",
var.equal = TRUE)
双样本 t 检验
data: Time 按 Topping 分组
t = 2.3811, df = 198, p-value = 0.01821
备择假设:Pepperoni 与 Cheese 组
均值之差不等于 0
95% 置信区间:
0.0599370 0.6377601
library(effectsize)
cohens_d(Time ~ Topping, data = Pizza)
Cohen's d | 95% CI
<----------------------
0.34 | [0.06, 0.62]
library(pwr)
pwr.t.test(n = 1000,
sig.level = 0.0182,
d = 0.34,
type = "two.sample")
双样本 t 检验功效计算
n = 100
d = 0.34
sig.level = 0.0182
power = 0.510256
alternative = two.sided
注意:n 为每组样本量
R 中的 A/B 测试