方差分析多重比较可视化

本文主要是演示，看一下方差分析中，如果进行多重比较，如何对多重比较的结果可视化 1。

1. 数据设置

首先，载入butron.maize这个数据 2

library(pacman)
p_load(agridat)
data("butron.maize")
dat = butron.maize
head(dat)

##          gen male female  env yield
## 1  A509xA637 A509   A637 pv95  6.81
## 2  A509xA661 A509   A661 pv95  5.02
## 3 A509xCM105 A509  CM105 pv95  7.39
## 4  A509xEP28 A509   EP28 pv95  6.84
## 5  A509xEP31 A509   EP31 pv95  6.04
## 6  A509xEP42 A509   EP42 pv95  7.10

看一下数据结构，确保因子和数字类型。3

str(dat)

## 'data.frame':    245 obs. of  5 variables:
##  $ gen   : Factor w/ 49 levels "A509xA637","A509xA661",..: 1 2 3 4 5 6 7 8 9 10 ...
##