第 12 章 Tidy Statistics

本章介绍基本的方差分析内容

12.1 从一个案例开始

从这是一份1994年收集1379个对象关于收入、身高、教育水平等信息的数据集,数据在课件首页下载。

首先,我们下载后导入数据

earn height sex race ed age
79571 73.89 male white 16 49
96397 66.23 female white 16 62
48711 63.77 female white 16 33
80478 63.22 female other 16 95
82089 63.08 female white 17 43
15313 64.53 female white 15 30

我们的问题:男性是否就比女性挣的多?

12.2 单因素方差分析

## 
##  Welch Two Sample t-test
## 
## data:  earn by sex
## t = -12, df = 768, p-value <2e-16
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -25324 -18171
## sample estimates:
## mean in group female   mean in group male 
##                24246                45993
## 
## Call:
## lm(formula = earn ~ sex, data = wages)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -46092 -20516  -4639  11722 271956 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    24246       1004    24.1   <2e-16 ***
## sexmale        21748       1636    13.3   <2e-16 ***
## ---
## Signif. codes:  
## 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 29400 on 1377 degrees of freedom
## Multiple R-squared:  0.114,  Adjusted R-squared:  0.113 
## F-statistic:  177 on 1 and 1377 DF,  p-value: <2e-16
##               Df   Sum Sq  Mean Sq F value Pr(>F)    
## sex            1 1.53e+11 1.53e+11     177 <2e-16 ***
## Residuals   1377 1.19e+12 8.66e+08                   
## ---
## Signif. codes:  
## 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

12.3 双因素方差分析

我们采用ggpubr宏包下的ToothGrowth来说明,这个数据集包含60个样本,记录着每10只豚鼠在不同的喂食方法和不同的药物剂量下,牙齿的生长情况.

  • len : 牙齿长度
  • supp : 两种喂食方法 (橙汁和维生素C)
  • dose : 抗坏血酸剂量 (0.5, 1, and 2 mg)

问题:豚鼠牙齿的长度是否与药物的食用方法和剂量有关?

线性回归时,我们是通过独立变量来预测响应变量,但现在我们关注的重点会从预测转向不同组别差异之间的分析,这即为方差分析(ANOVA)。

这里是两个解释变量,所以问题需要双因素方差分析 (ANOVA)

检验表明不同类型之间存在显著差异,但是并没有告诉我们具体谁与谁之间的不同。需要多重比较帮助我们解决这个问题。使用TurkeyHSD函数

思考:交互效应是否显著?