pythonic生物人

复现Science正刊图,换数据即可用!

本次使用R语言的两种方法,复现一下Science正刊中的几张图!

science.adk0775 Fig. 1 E
science.adk0775 Fig. 1 E
science.adk0775 Fig. 2 H
science.adk0775 Fig. 2 H
science.adk0775 Fig. 5 A
science.adk0775 Fig. 5 A
❤️复现效果图-案例1.1❤️
❤️复现效果图-案例1.1❤️
❤️复现效果图-案例1.2❤️
❤️复现效果图-案例1.2❤️
❤️复现效果图-案例2.1❤️
❤️复现效果图-案例2.1❤️
❤️复现效果图-案例2.2❤️
❤️复现效果图-案例2.2❤️

✅图表解读

该图为气泡图(Dot Plot),常用于展示单细胞RNA测序数据中特定基因在不同细胞群体中表达情况:

图表元素
代表含义
解读
气泡颜色
基因在细胞群体中的平均表达水平(经Z-score标准化)
越偏红色表示表达量高于该基因在所有细胞中的平均水平;越偏蓝色表示表达量低于平均水平;白色接近平均水平。
气泡大小
表达该基因的细胞在该细胞群体中所占的比例(百分比)
点越大,表示该基因在该细胞群体中表达越普遍。
Y轴(基因)
你要展示的基因列表
每个点代表一个基因。
X轴(细胞群体)
根据celltype划分的不同细胞群体(如"B细胞Normal", "T细胞Tumor")
每个点代表一个细胞亚群。

注意,x和y轴有时候会交换。


✅读入测试数据

Image

✅方法1,关键代码

一行代码,很简洁,需要输入Seurat对象,

# 方法1,关键代码
scDotPlot(
    object = dotsingle_data,           
    features = gene_list, 
    group = "celltype_group"
) 
Image

这里添加了xy轴上基因和细胞群的聚类树。

另外,推荐几组经典配色,

"#F5F5F0", "#B7D1A3", "#35635B"
"#F5F5F0", "#B7D1A3", "#35635B"
"#4575B4", "#FFFFFF", "#D73027"
"#4575B4", "#FFFFFF", "#D73027"
"#808080", "#FFFFFF", "#000000"
"#808080", "#FFFFFF", "#000000"
"#F0E442", "#009E73", "#0072B2"
"#F0E442", "#009E73", "#0072B2"

✅方法2,关键代码

代码依旧很简单,

# 方法1,关键代码
jjDotPlot(object = dotsingle_data,
          gene = gene_list,
          id = 'celltype',
          split.by = 'groups'
         )
Image

点可以使用不同的形状,

shape = 23
shape = 23
shape = 16
shape = 16

本期结束!

测试数据+详细代码,后期会加入👉《保姆级R可视化教程》

图片
《👉保姆级R可视化教程》,28个章节,11w字,数百张图,旨在引导如何系统学习R语言可视化,部分展示,
图片
图片
图片
图片
复现效果图-abcd图❤️
复现效果图-b图❤️
图片
❤️复现效果图-ABCD图❤️
图片
❤️复现效果图-bcd图❤️
图片
图片
❤️复现效果图-b图❤️
图片
图片
图片
图片
图片
图片
图片
图片
图片
图片

图片

图片
图片
图片

(加入学习,收费,备注:299)

图片