需完成作业
image.png
原始数据
image.png
data <- read.delim("D:/R/R 学习/rwadata/week1-3.txt")
p1 <- ggplot(data,aes(x=X,y=Y,color=Group))+
geom_line(size = .8)+
geom_point(shape = 21,color = 'black',size=3,fill = 'white')+
scale_color_d3()+
theme_classic()+
scale_color_npg()+
theme(legend.position = c(.1,.8))
p2 <- ggplot(data = data) +
geom_bar( mapping = aes(x = Group, fill= Group),width = 1)+
theme(plot.title = element_text(color="grey",size=12,hjust = 0.5,vjust = 1),
legend.position ="top",
axis.title.y = element_blank())+#+
scale_color_npg()+
labs(title="Isis by group",x = NULL, y = NULL)+
coord_polar(theta="x")+
theme_light()+
theme(plot.title = element_text(color="grey",size=12,hjust = 0.5,vjust = 1),
legend.position ="top",
axis.title.y = element_blank())+
scale_color_npg()
p2
#ggbetweenstats(),用于绘制组间比较的箱线图、小提琴图及小提琴箱线组合图，并提供详细的参数检验和统计信息
p3 <- ggbetweenstats(
data = data,
x = Group,
y = Y,
plot.type = "violin"#图标类型，不加参数则默认为boxviolin，其它可选为：box、violin
# title = "Distribution of sepal length across Iris species"
)+
scale_color_npg()
p3
#提取统计数据
ggbetweenstats(data, Group, Y) %>%
extract_stats()
# $subtitle_data
# NULL
#
# $caption_data
# NULL
#
# $pairwise_comparisons_data
# # A tibble: 6 x 9
# group1 group2 statistic
# <chr> <chr> <dbl>
# 1 Group1 Group2 -13.4
# 2 Group1 Group3 -7.19
# 3 Group1 Group4 -8.00
# 4 Group2 Group3 2.11
# 5 Group2 Group4 4.59
# 6 Group3 Group4 1.16
# p.value alternative
# <dbl> <chr>
# 1 0.00000000283 two.sided
# 2 0.000733 two.sided
# 3 0.00000665 two.sided
# 4 0.917 two.sided
# 5 0.0279 two.sided
# 6 0.917 two.sided
# distribution p.adjust.method
# <chr> <chr>
# 1 q Holm
# 2 q Holm
# 3 q Holm
# 4 q Holm
# 5 q Holm
# 6 q Holm
# test expression
# <chr> <list>
# 1 Games-Howell <language>
# 2 Games-Howell <language>
# 3 Games-Howell <language>
# 4 Games-Howell <language>
# 5 Games-Howell <language>
# 6 Games-Howell <language>
#
# $descriptive_data
# NULL
#
# $one_sample_data
# NULL
#
# $tidy_data
# NULL
#
# $glance_data
# NULL
p4 <- ggbetweenstats(
data = data,
x = Group,
y = Y,
plot.type = "violin",
title = "Distribution of sepal length across Iris species",
results.subtitle = TRUE,#决定是否将统计检验的结果显示为副标题（默认TRUE）;如果设置为FALSE,则仅返回绘图
subtitle = TRUE,#副标题,默认显示统计检验结果,自定义则results.subtitle=FALSE
outlier.tagging = TRUE,#是否标记离群异常值，默认FALSE
outlier.shape = 19,#异常值形状,可设置为NA将其隐藏（不是删除，因此不会影响统计检验结果）
outlier.color = "pink",#异常值颜色
outlier.label.args = list(size = 4),#异常值标签大小
point.args = list(position = ggplot2::position_jitterdodge(dodge.width = 0.6),
alpha = 0.5, size = 3, stroke = 0),#传递给geom_point的参数设置
violin.args = list(width = 0.4, alpha = 0.2),#传递给geom_violin的参数设置
ggtheme = theme_bw(),#主题修改，可直接调用ggplot2的主题，默认主题为ggstatsplot::theme_ggstatsplot()
package = "ggsci",#提取调色板所需的包
palette = "uniform_startrek")#选择提取包中的调色板
p4
pdf("result/homework2.pdf",width=10,height = 12)
plot_grid(p1,p2,p3,ncol=2,
labels = LETTERS[1:3],
# labels = "AUTO",
align = 'h')
dev.off()
pdf("result/homework2-2.pdf")
p1+p2+p3
dev.off()
p4 <- plot_grid(p1,p2,p3,ncol=2,
labels = LETTERS[1:3],
# labels = "AUTO",
align = 'h')
p4
ggsave("result/homework2-2.png",p4,dpi=600,width=8,height = 12)
完成情况
image.png
仍需改进，欢迎提供更改意见
参考代码
[如何一键绘制具有显著性检验信息的箱线图/小提琴图？](如何一键绘制具有显著性检验信息的箱线图/小提琴图？ (qq.com)
（[Polar coordinates — coord_polar • ggplot2 (tidyverse.org)]
https://ggplot2.tidyverse.org/reference/coord_polar.html
