泛癌分析数据库

❝ 已经付费加群的小伙伴无需二次付费,等待师兄后续更新即可! ❞ 封面 从这个系列开始,师兄就带着大家从各大顶级期刊中的Figuer入手,从仿照别人的作图风格到最后实现自己游刃有余的套用在自己的分析数据上!这一系列绝对是高质量!还不赶紧「点赞+在看」 ,

已经付费加群的小伙伴无需二次付费,等待师兄后续更新即可!

deb4235f41a9f4e4a78845494d3e838b.png

封面

从这个系列开始,师兄就带着大家从各大顶级期刊中的Figuer入手,从仿照别人的作图风格到最后实现自己游刃有余的套用在自己的分析数据上!这一系列绝对是高质量!还不赶紧「点赞+在看」,学起来!

本期分享的是期刊:「Nucleic Acids Research(IF = 19.16)」上面一篇文章中的一个「泛癌箱线图+显著性检验」

本系列所有代码和示例数据将会和生信常用图形系列绘图放在一起,扫描下方二维码添加师兄微信,「付费169元」,即可加入生信绘图交流群。「群内不仅提供生信常用图形系列的代码,还会提供本系列后续所有Figure的实例数据和代码,我会在文章更新后第一时间上传。」

9694bf6072c45fe609253d7e046e8ffd.png

师兄微信

当然了!如果你还想白嫖,师兄的文章中代码已经写的很清楚了!但是师兄还是希望你点个赞再走呗!

以上就是本期的全部内容啦!「欢迎点赞,点在看!」师兄会尽快更新哦!制作不易,你的打赏将成为师兄继续更新的十足动力!

「优惠方式:点赞+在看,并转发这两个系列任意一篇文章至朋友圈,集赞30个,即可享受¥139入群!」

d6e267d4c5f4e889cb4493aa1358a653.png

参考文献

话不多说,直接上图!

读图

307081df3731c8b903092bacb227ebcc.png

原图

本期图形来自于数据库CancerMiRNome数据库,数据均可在数据库中下载,图形难度不大,但是在箱线图外散点、以及显著性检验的细节上还是有些值得学习的地方!

效果展示

5d8d946dcc931dd1d6be29bd17d28fcd.png

效果展示

本次复现完全在R语言中进行,并且完美的解决了以上所有问题!请大家放心食用!

绘图群附加福利

凡是「已经加群」的小伙伴,你们在看文献的时候如果看到好看的Figure,可以发到群里!师兄会及时关注的,「如果被师兄选中,就会在推文中更新!」

Figure的要求如下(被选中的前提):

  • 首先肯定是要符合大众审美的,在图形样式上要过关。

  • 新颖独特,有与常见图形不一样的地方。

  • 有一定难度,太简单的大家都会做,没什么挑战性哈!

往期选中案例

0a1523672d4b0a0d082bcd389cf5150d.png

往期选中案例01

364854f0cf237b66d4099d98ec742bba.png

往期选中案例

R包载入及数据读取

# 载入R包:
library(ggplot2)
library(ggpubr)
library(tidyverse)
library(latex2exp)
library(rstatix)

# 读取数据:
data <- read.csv("hsa-let-7a-5p.TCGA_PanCancer_Expression_Data.csv", header = T)
head(data)
#             mir   project          sample group   expr
# 1 hsa-let-7a-5p TCGA-BRCA TCGA-3C-AAAU-01 Tumor 15.171
# 2 hsa-let-7a-5p TCGA-BRCA TCGA-3C-AALI-01 Tumor 15.494
# 3 hsa-let-7a-5p TCGA-BRCA TCGA-3C-AALJ-01 Tumor 15.453
# 4 hsa-let-7a-5p TCGA-BRCA TCGA-3C-AALK-01 Tumor 15.592
# 5 hsa-let-7a-5p TCGA-BRCA TCGA-4H-AAAK-01 Tumor 15.967
# 6 hsa-let-7a-5p TCGA-BRCA TCGA-5L-AAT0-01 Tumor 16.044

基础绘图

# 绘图:
ggplot(data)+
  # 基础图形:
  geom_boxplot(aes(project, expr, fill = group))+
  # 设置颜色:
  scale_fill_manual(values = c("Tumor" = "#d6503a", "Normal" = "#5488ef"))+
  # 设置主题:
  theme_bw()+
  theme(axis.text = element_text(face = "bold"),
        axis.text.x = element_text(angle = 45, hjust = 1))

43190fea97ca04a9ea25cece84465f14.png

基础绘图

完善美化

# 按中位数由高到低排列:
data_new <- data %>% 
  group_by(project) %>% 
  mutate(median = median(expr), group_max = max(expr)) %>% 
  arrange(desc(median))

# 调整因子顺序:
data_new$project <- factor(data_new$project, levels = unique(data_new$project))

data_new1 <- data_new %>% 
  group_by(project) %>% 
  mutate(group_number = length(unique(group)))

# 此数据用于添加显著性检验的label:
stat_data <- data_new1[which(data_new1$group_number == 2),] %>% 
  group_by(project) %>% 
  wilcox_test(expr ~ group) %>% 
  adjust_pvalue(method = "fdr") %>%
  add_significance('p')

data_new1 <- unique(data_new1[,c(2, 7,8)])

index <- which(data_new1$group_number == 2)

# 此数据用于画显著性检验的括号用:
bracket_data <- data.frame(x = index - 0.3,
                           xend = index + 0.3,
                           y = as.numeric(unlist(data_new1[index, 2]))+0.5,
                           annotation = "")

# 重新画图:
ggplot()+
  # 基础图形:
  geom_boxplot(data = data_new, aes(x = project, y =expr, fill = group),outlier.shape = 21, outlier.fill = "white")+
  # 设置颜色:
  scale_fill_manual(name = NULL, values = c("Tumor" = "#d6503a", "Normal" = "#5488ef"))+
  # 设置主题:
  theme_bw()+
  labs(y = TeX("miRNA Level ($log_{2}$CPM)", bold = T))+
  theme(axis.text = element_text(face = "bold"),
        axis.text.x = element_text(angle = 45, hjust = 1),
        axis.line = element_line(colour = "black"),
        panel.border = element_blank(),
        panel.background = element_blank(),
        legend.position = c(0.99, 0.99), 
        legend.justification = c(1,1))+
  # 显著性检验:
  geom_signif(stat = "identity",
              position = "identity",
              data = bracket_data,
              aes(x = x,
                  xend = xend,
                  y = y,
                  yend = y,
                  annotation = annotation))+
  annotate(geom = "text", x=index, y = bracket_data$y+0.25, label = as.character(stat_data$p.signif))

ggsave("box_plot.pdf", height = 6, width = 10)

54ec85dc423089860204d4bc780c684f.png

效果展示

往期优秀图形目录

bea954b575506093ff3fb35eb53d4d89.png

渐变火山图

e0c68e47fe7aaf7441e15613b659f905.png

气泡图+相关性热图

7cc99d890c9f90e653ebe7a615152065.png

复杂提琴图

2c92e2d23ee1176b058430754e087996.png

复杂热图

25e3a3f0e16d30849dcc0a68a84828bc.png

复杂散点图

75a2cf9dd62c402d1e2c88284cc3f2c8.png

复杂热图02

f8ba6165e8e112a59f6ae8ebf58e8301.png

甘特图

9e9ad056d7ffbc9b06e0c126427a6144.png

百分比柱状图

2ef3d8f6af52cc2fbb0764ddc3bbbbb5.png

箱线图美化

652c26323f50a32de5d7f23fc4420409.png

弦图

77c47ac1974932cc86ba8559869f8fe2.png

mantel test图

5611a6d5e3aec41c5277a6d012b259d0.png

瀑布图

2325edc9e31d300cd5fbf9dafa8918ba.png

曼哈顿图

a71ca26da70b206afa0dc632a08ec67c.png

KEGG富集图

648577518534b5d6b7f14aa1375ad7ea.png

哑铃图

以上内容仅为群内部分内容,不代表全部。详细目录请看下方列表!已经入群的小伙伴,无需付费购买推文,群内都会及时更新!

示例数据和代码获取

本系列所有代码和示例数据将会和生信常用图形系列绘图放在一起,公众号右下角添加师兄微信,「付费169元」,即可加入生信绘图交流群。「群内不仅提供生信常用图形系列的代码,还会提供本系列后续所有Figure的实例数据和代码,我会在文章更新后第一时间上传。」

当然了!如果你还想白嫖,师兄的文章中代码已经写的很清楚了!但是师兄还是希望你点个赞再走呗!

以上就是本期的全部内容啦!「欢迎点赞,点在看!」师兄会尽快更新哦!制作不易,你的打赏将成为师兄继续更新的十足动力!

「优惠方式:点赞+在看,并转发这两个系列任意一篇文章至朋友圈,集赞30个,即可享受¥139入群!」

往期文章

  1. 跟着Nature Medicine学作图 -- 箱线+散点图

  2. 跟着Nature Communications学作图 -- 渐变火山图

  3. 跟着Nature Communications学作图 -- 气泡图+相关性热图

  4. 跟着Nature Communications学作图 -- 复杂提琴图

  5. 跟着Nature Medicine学作图 -- 复杂热图

  6. 跟着Nature Com学作图 -- 复杂散点图

  7. 跟着高分SCI学作图 -- 气泡火山图

  8. 跟着高分SCI学作图--复杂热图02

  9. 跟着Nature Communications学作图 -- 甘特图

  10. 跟着Nature Communications学作图 -- 双向柱状图

  11. 跟着Nature Communications学做图 -- 批量散点图+拟合曲线

  12. 跟着Nature Communications学作图 -- 箱线图美化

  13. 跟着Nature Communications学作图 -- 复杂百分比柱状图

  14. 跟着Nature Communications学作图 -- 渐变柱状图

  15. 跟着Molecular Cancer学作图 -- 分半小提琴图

  16. 跟着高分SCI学作图 -- 复杂热图+渐变色连线

  17. 跟着Nature Communications学作图 -- 复杂热图+堆积柱状图注释

  18. 跟着Cell学作图 -- 单细胞分组火山图

  19. 跟着高分SCI学作图 -- 堆积柱状图+折线+热图注释

  20. 跟着Nucleic Acids Research学作图 -- 多组学九象限散点图

  21. 跟着Nature Communications作图 -- 分组散点图+直方图注释

  22. 跟着Nature Medicine学作图 -- 复杂环形热图

  23. 生信常用分析图形绘制01 -- 各种类型的热图你学会了吗?

  24. 生信常用分析图形绘制02 -- 解锁火山图真谛!

  25. 生信常用分析图形绘制03 -- 富集分析圈图

  26. 生信常用分析图形绘制04 -- 桑基图

  27. 生信常用分析图形绘制05 -- 弦图

  28. 生信常用分析图形绘制06 -- 富集分析分组柱状图+气泡图

  29. 生信常用分析图形绘制07 -- PCA图

  30. 生信常用分析图形绘制08--云雨图

  31. 生信常用分析图形绘制09 -- Venn图

  32. 生信常用分析图形绘制10 -- 环形互作网络图

  33. 生信常用分析图形绘制11 -- 相互作用网络图

  34. 生信常用分析图形绘制12 -- 聚类树美化

  35. 生信常用分析图形绘图13 -- 富集进阶气泡图

  36. 生信常用分析图形绘制14 -- mantel test相关性图

  37. 生信常用分析图形绘制15 -- 词云图

  38. 生信常用分析图形绘制16 -- 瀑布图

  39. 生信常用分析图形绘制17 -- 森林图

  40. 生信常用分析图形绘制18 -- 曼哈顿图

  41. 生信常用分析图形绘制19 -- KEGG分组柱状图

  42. 生信常用分析图形绘制20 -- 哑铃图

  43. 生信常用分析图形绘制21 -- 三线表绘制及美化

  44. 生信常用分析图形绘制22 -- 嵌套圈图

  45. 生信常用分析图形绘制23 -- 列线图

  46. 生信常用分析图形绘制24-- 蜂群图

  47. 生信常用分析图形绘制25 -- 箱线图+贝塞尔曲线

  48. 生信常用分析图形绘制26 -- 相关性矩阵图+矩阵散点图

  49. 生信常用分析图形绘制27 -- 交互式热图+局部放大

  50. R语言生信绘图群!轻松解决你的生信绘图问题!

知秋君
上一篇 2024-08-17 09:02
下一篇 2024-08-17 08:36

相关推荐