使用ggboxplot函数绘制箱线图
加载所需R包
library(ggpubr)
基本用法：
Usage
ggboxplot(data, x, y, combine = FALSE, merge = FALSE, color = "black",
fill = "white", palette = NULL, title = NULL, xlab = NULL,
ylab = NULL, facet.by = NULL, panel.labs = NULL,
short.panel.labs = TRUE, linetype = "solid", size = NULL, width = 0.7,
notch = FALSE, select = NULL, remove = NULL, order = NULL,
add = "none", add.params = list(), error.plot = "pointrange",
label = NULL, font.label = list(size = 11, color = "black"),
label.select = NULL, repel = FALSE, label.rectangle = FALSE,
ggtheme = theme_pubr(), ...)
常用参数
Arguments
data #a data frame
x #character string containing the name of x variable.
y #character vector containing one or more variables to plot
combine #logical value. Default is FALSE. Used only when y is a vector containing multiple variables to plot. If TRUE, create a multi-panel plot by combining the plot of y variables.
merge #logical or character value. Default is FALSE. Used only when y is a vector containing multiple variables to plot. If TRUE, merge multiple y variables in the same plotting area. Allowed values include also "asis" (TRUE) and "flip". If merge = "flip", then y variables are used as x tick labels and the x variable is used as grouping variable.
color #outline color.
fill #fill 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".
title #plot main title.
xlab #character vector specifying x axis labels. Use xlab = FALSE to hide xlab.
ylab #character vector specifying y axis labels. Use ylab = FALSE to hide ylab.
facet.by #character vector, of length 1 or 2, specifying grouping variables for faceting the plot into multiple panels. Should be in the data.
panel.labs #a list of one or two character vectors to modify facet panel labels. For example, panel.labs = list(sex = c("Male", "Female")) specifies the labels for the "sex" variable. For two grouping variables, you can use for example panel.labs = list(sex = c("Male", "Female"), rx = c("Obs", "Lev", "Lev2") ).
short.panel.labs #logical value. Default is TRUE. If TRUE, create short labels for panels by omitting variable names; in other words panels will be labelled only by variable grouping levels.
linetype #line types.
size #Numeric value (e.g.: size = 1). change the size of points and outlines.
width #numeric value between 0 and 1 specifying box width.
notch #是否添加缺口 if FALSE (default) make a standard box plot. If TRUE, make a notched box plot. Notches are used to compare groups; if the notches of two boxes do not overlap, this suggests that the medians are significantly different.
select #character vector specifying which items to display.
remove #character vector specifying which items to remove from the plot.
order #character vector specifying the order of items.
add #character vector for adding another plot element (e.g.: dot plot or error bars). Allowed values are one or the combination of: "none", "dotplot", "jitter", "boxplot", "point", "mean", "mean_se", "mean_sd", "mean_ci", "mean_range", "median", "median_iqr", "median_mad", "median_range"; see ?desc_statby for more details.
add.params #parameters (color, shape, size, fill, linetype) for the argument 'add'; e.g.: add.params = list(color = "red").
error.plot #plot type used to visualize error. Allowed values are one of c("pointrange", "linerange", "crossbar", "errorbar", "upper_errorbar", "lower_errorbar", "upper_pointrange", "lower_pointrange", "upper_linerange", "lower_linerange"). Default value is "pointrange" or "errorbar". Used only when add != "none" and add contains one "mean_*" or "med_*" where "*" = sd, se, ....
label #the name of the column containing point labels. Can be also a character vector with length = nrow(data).
font.label #a list which can contain the combination of the following elements: the size (e.g.: 14), the style (e.g.: "plain", "bold", "italic", "bold.italic") and the color (e.g.: "red") of labels. For example font.label = list(size = 14, face = "bold", color ="red"). To specify only the size and the style, use font.label = list(size = 14, face = "plain").
repel #a logical value, whether to use ggrepel to avoid overplotting text labels or not.
label.rectangle #logical value. If TRUE, add rectangle underneath the text, making it easier to read.
ggtheme #function, ggplot2 theme name. Default value is theme_pubr(). Allowed values include ggplot2 official themes: theme_gray(), theme_bw(), theme_minimal(), theme_classic(), theme_void()
... #other arguments to be passed to geom_boxplot, ggpar and facet.
使用示例
Examples
# Load data
data("ToothGrowth")
df <- ToothGrowth
head(df)
## len supp dose
## 1 4.2 VC 0.5
## 2 11.5 VC 0.5
## 3 7.3 VC 0.5
## 4 5.8 VC 0.5
## 5 6.4 VC 0.5
## 6 10.0 VC 0.5
# Basic plot
# width: change box plots width
p1 <- ggboxplot(df, x = "dose", y = "len", width = 0.8)
# Change orientation: horizontal
p2 <- ggboxplot(df, "dose", "len", orientation = "horizontal")
p2
p2
# Notched box plot 添加缺口
p3 <- ggboxplot(df, x = "dose", y = "len",
notch = TRUE)
p3
p3
# Add dots
p4 <- ggboxplot(df, x = "dose", y = "len",
add = "dotplot")
p4
p4
# Add jitter points and change the shape by groups
p5 <- ggboxplot(df, x = "dose", y = "len",
add = "jitter", shape = "dose")
p5
p5
# Select and order items 选择特定列和指定特定排列顺序
# Select which items to display: "0.5" and "2"
p6 <- ggboxplot(df, "dose", "len",
select = c("0.5", "2"))
p6
p6
# Change the default order of items
p7 <- ggboxplot(df, "dose", "len",
order = c("2", "1", "0.5"))
p7
p7
# Change colors 更改填充色和边框色
# Change outline and fill colors
p8 <- ggboxplot(df, "dose", "len",
color = "black", fill = "gray")
p8
p8
# Change outline colors by groups: dose
# Use custom color palette
# Add jitter points and change the shape by groups
p9 <- ggboxplot(df, "dose", "len",
color = "dose", palette =c("#00AFBB", "#E7B800", "#FC4E07"),
add = "jitter", shape = "dose")
p9
p9
# Change fill color by groups: dose
p10 <- ggboxplot(df, "dose", "len",
fill = "dose", palette = c("#00AFBB", "#E7B800", "#FC4E07"))
p10
p10
# Box plot with multiple groups
# fill or color box plot by a second group : "supp"
p11 <- ggboxplot(df, "dose", "len", color = "supp",
palette = c("#00AFBB", "#E7B800"))
p11
p11
p12 <- ggboxplot(df, "dose", "len", color= "supp", facet.by = "supp",
palette = c("#00AFBB", "#E7B800"))
p12
p12
参考来源：
https://www.rdocumentation.org/packages/ggpubr/versions/0.1.4/topics/ggboxplot
sessionInfo()
## R version 3.5.1 (2018-07-02)
## Platform: x86_64-apple-darwin15.6.0 (64-bit)
## Running under: OS X El Capitan 10.11.3
##
## Matrix products: default
## BLAS: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRblas.0.dylib
## LAPACK: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRlapack.dylib
##
## locale:
## [1] zh_CN.UTF-8/zh_CN.UTF-8/zh_CN.UTF-8/C/zh_CN.UTF-8/zh_CN.UTF-8
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] ggpubr_0.1.7.999 magrittr_1.5 ggplot2_3.0.0
##
## loaded via a namespace (and not attached):
## [1] Rcpp_0.12.18 rstudioapi_0.7 bindr_0.1.1 knitr_1.20
## [5] tidyselect_0.2.4 munsell_0.5.0 colorspace_1.3-2 R6_2.2.2
## [9] rlang_0.2.2 stringr_1.3.1 plyr_1.8.4 dplyr_0.7.6
## [13] tools_3.5.1 grid_3.5.1 gtable_0.2.0 withr_2.1.2
## [17] htmltools_0.3.6 assertthat_0.2.0 yaml_2.2.0 lazyeval_0.2.1
## [21] rprojroot_1.3-2 digest_0.6.16 tibble_1.4.2 crayon_1.3.4
## [25] bindrcpp_0.2.2 purrr_0.2.5 glue_1.3.0 evaluate_0.11
## [29] rmarkdown_1.10 labeling_0.3 stringi_1.2.4 compiler_3.5.1
## [33] pillar_1.3.0 scales_1.0.0 backports_1.1.2 pkgconfig_2.0.2
