❝VIP群里有观众老爷询问如何对箱线图进行几何形状填充，那么今天就来具体介绍一番；在原有的基础上做了一些小的改动也许恰好您正好有此特殊需求，需要着重体会八个字「变实为虚，变虚为实」，希望对各位观众老爷有所帮助；下面来看具体案例； ❞
加载R包
代码语言：javascript复制library(tidyverse)
library(ggsci)
library(ggprism)
library(rstatix)
library(ggpubr)
library(ggpmisc)
library(ggpattern)
加载数据
代码语言：javascript复制gapminder <- read_tsv("gapminder.xls")
数据清洗
代码语言：javascript复制df <- gapminder %>%
filter(year %in% c(1957,2002,2007),continent !="Oceania") %>%
select(country,year,lifeExp,continent)%>%
mutate(paired = rep(1:(n()/3),each=3),year=factor(year))
统计分析
代码语言：javascript复制df_p_val1 <- df %>% group_by(continent) %>%
wilcox_test(lifeExp ~ year) %>%
adjust_pvalue(p.col = "p", method = "bonferroni") %>%
add_significance(p.col = "p.adj") %>%
add_xy_position(x = "year", dodge = 0.8)
构建填充类型
代码语言：javascript复制❝下面函数来自「ggpattern」包官方文档 https://coolbutuseless.github.io/package/ggpattern/articles/developing-patterns.html ❞
tiling3_pattern <- function(params, boundary_df, aspect_ratio, legend = FALSE) {
args <- as.list(params)
args <- args[grep("^pattern_", names(args))]
# hexagonal tiling using "regular_polygon" pattern
args$pattern <- "polygon_tiling"
# three-color tiling using `fill`, `pattern_fill` and their "average"
avg_col <- gridpattern::mean_col(params$fill, params$pattern_fill)
args$pattern_fill <- c(params$fill, avg_col, args$pattern_fill)
args$x <- boundary_df$x
args$y <- boundary_df$y
args$id <- boundary_df$id
args$prefix <- ""
do.call(gridpattern::patternGrob, args)
}
options(ggpattern_geometry_funcs = list(tiling3 = tiling3_pattern))
数据可视化
代码语言：javascript复制df %>%
ggplot(aes(year,lifeExp))
stat_boxplot(aes(ymin = ..lower.., ymax = ..upper..),outlier.shape = NA,width=0.5)
stat_boxplot(geom = "errorbar", aes(ymin = ..ymax..),width=0.2,size=0.35)
stat_boxplot(geom = "errorbar", aes(ymax = ..ymin..),width=0.2,size=0.35)
geom_boxplot_pattern(aes(pattern = year,pattern_angle = year,fill=year),
pattern = 'tiling3', pattern_angle = 45,color="black",
pattern_spacing = 0.045,outlier.shape = NA,linetype = "dashed",width=0.5,size=0.35)
stat_summary(geom = "crossbar", fun = "median",width = 0.5,color="black",size=0.38)
stat_pvalue_manual(df_p_val1,label = "p.adj.signif",label.size=5,hide.ns = F)
scale_size_continuous(range=c(1,3))
# geom_smooth(method = "lm", formula = NULL,size=1,se=T,color="black",linetype="dashed",aes(group=1))
stat_cor(label.y = 25,aes(label = paste(..rr.label.., ..p.label.., sep = "~`,`~"),group=1),color="black",
label.x.npc = "left")
stat_regline_equation(label.y = 19,aes(group=1),color="black")
facet_wrap(.~continent,nrow=1)
scale_fill_npg()
scale_x_discrete(guide = "prism_bracket")
scale_y_continuous(limits = c(0,95),minor_breaks = seq(0,95,5),guide = "prism_offset_minor")
labs(x=NULL,y=NULL)
theme_prism(base_line_size =0.4)
theme(plot.margin=unit(c(0.5,0.5,0.5,0.5),units=,"cm"),
strip.text = element_text(size=12),
axis.line = element_line(color = "black",size = 0.4),
panel.grid.minor = element_blank(),
panel.grid.major = element_line(size = 0.2,color = "#e5e5e5"),
axis.text.y = element_text(color="black",size=10),
axis.text.x = element_text(margin = margin(t = -5),color="black",size=10),
legend.position = "none",
panel.spacing = unit(0,"lines"))
coord_cartesian()
