copy from: http://www.reibang.com/p/3862faf22c14
# 清空環(huán)境變量
rm(list=ls())
# 安裝、載入包
install.packages("ggpubr")
library(ggpubr)
install.packages("ggsci")
library(ggsci)
# 加載數(shù)據(jù)集非凌，并查看
data("ToothGrowth")
str(ToothGrowth)
head(ToothGrowth)
image.png
image.png
ggboxplot函數(shù)
ggboxplot為ggpubr包的函數(shù)呛占。
ggboxplot(data, x, y, color = "black", bxp.errorbar = FALSE)
| data | a data frame |
|---|---|
| x | character string containing the name of x variable |
| y | character vector containing one or more variables to plot |
| color | outline color |
| palette | the color palette to be used for coloring or filling by groups. Allowed values include "grey" for grey color palettes; brewer palettes e.g. "RdBu", "Blues", ...; or custom color palette e.g. c("blue", "red"); and scientific journal palettes from ggsci R package, e.g.: "npg", "aaas", "lancet", "jco", "ucscgb", "uchicago", "simpsons" and "rickandmorty". |
| bxp.errorbar | logical value. If TRUE, shows error bars of box plots. |
#從lancet調(diào)色板提取9個(gè)顏色
pal = pal_lancet("lanonc")(9)
#選取前3個(gè)顏色
colors = pal_lancet("lanonc")(9)[1:3]
my_comparisons = list(c("0.5","1"),c("1","2"), c("0.5","2"))
p = ggboxplot( ToothGrowth , x = "dose", y = "len", color = "dose" , palette =colors, bxp.errorbar = T) +
stat_compare_means(comparisons =my_comparisons) +stat_compare_means(method = "anova" , label.x = 1.7, label.y = 45)
p
p + stat_boxplot(geom ="errorbar", color = colors, width =0.4, size =1) + geom_boxplot(aes(color = dose))
p
stat_compare_means
- ggpubr包的函數(shù)
- Add mean comparison p-values to a ggplot, such as box blots, dot plots and stripcharts.
- stat_compare_means(method = NULL, comparisons = NULL, label.x = NULL, label.y = NULL)
| method | a character string indicating which method to be used for comparing means. |
|---|---|
| comparisons | A list of length-2 vectors. The entries in the vector are either the names of 2 values on the x-axis or the 2 integers that correspond to the index of the groups of interest, to be compared. |
| label.x,label.y | numeric Coordinates (in data units) to be used for absolute positioning of the label. If too short they will be recycled. |
stat_boxplot
- ggplot2包函數(shù)晌姚，A box and whiskers plot (in the style of Tukey)
- The boxplot compactly displays the distribution of a continuous variable. It visualises five summary statistics (the median, two hinges and two whiskers), and all "outlying" points individually.
- stat_boxplot(geom ="errorbar", color = colors, width =0.4, size =1) + geom_boxplot(aes(color = dose))
結(jié)果展示：
image.png
