【生信分析】棒棒糖也能成图?R语言绘制棒棒糖图

棒棒糖图介绍

简单老来说,就是条形图的替代方法,但是比柱状图能展示更加丰富的信息

直接上图,给大家展示一下我画的棒棒糖图:

棒棒糖图基本画法

由分组变量“cyl”着色的棒棒糖图表:

ggdotchart(dfm, x = "name", y = "mpg",
           color = "cyl",                                # 按组显示颜色
           palette = c("#00AFBB", "#E7B800", "#FC4E07"), # 自定义调色板
           sorting = "ascending",                        # 按降序对值排序
           add = "segments",                             # 添加从y=0到点的线段
           ggtheme = theme_pubr()                        # 主题
           )

添加 排序,旋转,点大小,数字

  • 按降序排序:sorting = "descending"。
  • 垂直旋转绘图: rotate = = TRUE 
  • 每个组内的 mpg 值进行排序:group = "cyl"
  • 点大小设置: dot.size=6
  • 添加 mpg 值作为标签。 label = "mpg" 或者 label = round(dfm$mpg)
ggdotchart(dfm, x = "name", y = "mpg",
           color = "cyl",                                # 按组显示颜色
           palette = c("#00AFBB", "#E7B800", "#FC4E07"), # 自定义调色板
           sorting = "descending",                       # 按降序对值排序
           add = "segments",                             # 添加从y=0到点的线段
           rotate = TRUE,                                # 垂直旋转
           group = "cyl",                                # 按组排序
           dot.size = 6,                                 # 点大小
           label = round(dfm$mpg),                       # 将mpg值添加为点标签
           font.label = list(color = "white", size = 9, 
                             vjust = 0.5),               # 调整标签参数
           ggtheme = theme_pubr()                        # ggplot2主题
           )

添加 偏差

  • 更改分段颜色和大小: add.params = list(color = “lightgray”, size = 2)
ggdotchart(dfm, x = "name", y = "mpg_z",
           color = "cyl",                                # 按组显示颜色
           palette = c("#00AFBB", "#E7B800", "#FC4E07"), # 自定义调色板
           sorting = "descending",                       # 按降序排序值
           add = "segments",                             # 添加从y=0到点的分段
           add.params = list(color = "lightgray", size = 2), # 添加从y=0到点的分段
           group = "cyl",                                # 按组排序
           dot.size = 6,                                 # 点大小
           label = round(dfm$mpg_z,1),                   # 将mpg值添加为点标签
           font.label = list(color = "white", size = 9, 
                             vjust = 0.5),               # 调整标签参数
           ggtheme = theme_pubr()                        # ggplot2主题
           )+
  geom_hline(yintercept = 0, linetype = 2, color = "lightgray")

分组配色

  • 按组对 y 文本进行着色(使用 y.text.col = TRUE):
 
 
ggdotchart(dfm, x = "name", y = "mpg",
           color = "cyl",                                # 分组颜色
           palette = c("#00AFBB", "#E7B800", "#FC4E07"), # 自定义调色板
           sorting = "descending",                       # 按降序对值进行排序
           rotate = TRUE,                                # 垂直旋转
           dot.size = 2,                                 # 点尺寸
           y.text.col = TRUE,                            # 按组为y文本上色
           ggtheme = theme_pubr()                        # ggplot2主题
           )+
  theme_cleveland()                                      # 添加虚线网格

想要了解更多可以访问 :

kassambara/ggpubr: 'ggplot2' Based Publication Ready Plots (github.com) 

  • 12
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
### 回答1: GSEA(基因集富集分析)是一种常用的生物信息学分析方法,用于研究基因集在基因表达谱中的富集情况。下面是使用R语言进行GSEA生信分析的代码示例: 1. 首先,需要安装和加载必要的R包,例如GSEA包和其他必要的依赖包。 ```R install.packages("GSEA") library(GSEA) ``` 2. 加载基因表达数据集,通常是一个包含基因表达矩阵的数据文件。假设文件名为"expression_data.txt",其中包含基因表达矩阵和对应的样本信息。 ```R expression_matrix <- read.table("expression_data.txt", header = TRUE) ``` 3. 定义基因集,可以是预定义的基因集数据库(例如MSigDB)中的基因集,也可以是自定义的基因集。 ```R gene_sets <- c("GO_Biological_Process", "KEGG_Pathways", "Custom_Gene_Set") ``` 4. 进行GSEA分析,使用`gsea()`函数。其中,`gene_expr_matrix`参数为基因表达矩阵,`gene_sets`参数为基因集,`class_vector`参数为样本类别信息向量。 ```R gsea_results <- gsea(gene_expr_matrix = expression_matrix, gene_sets = gene_sets, class_vector = sample_classes) ``` 5. 分析结果包括富集分数(Enrichment Score)、正负富集基因集和富集谱等。可以通过可视化方法进一步探索和解释这些结果。 ```R enrichment_score <- gsea_results$es positive_sets <- gsea_results$pos_sets negative_sets <- gsea_results$neg_sets gene_set_plot <- plot(gsea_results) ``` 以上是使用R语言进行GSEA生信分析的基本代码示例。根据具体的研究问题和分析目标,还可以进行更多的数据预处理和可视化分析。 ### 回答2: GSEA(Gene Set Enrichment Analysis)是一种生物信息学分析工具,可用于确定基因集在给定基因表达数据中的富集程度。下面是R语言中实现GSEA分析的示例代码。 首先,需要安装并加载GSEABase、clusterProfiler和enrichplot等相关的R包。 ```R install.packages("GSEABase") install.packages("clusterProfiler") install.packages("enrichplot") library(GSEABase) library(clusterProfiler) library(enrichplot) ``` 接下来,准备基因表达数据和基因集数据。假设基因表达数据保存在一个矩阵中,行表示基因,列表示样本;基因集数据保存在GMT格式文件中,每行包含一个基因集的名称、描述和基因列表。 ```R expression_data <- read.table("expression_data.txt", header = TRUE, row.names = 1) gmt_file <- system.file("extdata", "c2.cp.kegg.v7.4.symbols.gmt", package = "DOSE") gene_sets <- readGMT(gmt_file) ``` 然后,进行GSEA分析。可以选择使用差异表达基因列表作为输入,或者将基因表达数据与基因集数据一起传递。以下是基于基因表达数据进行GSEA分析的示例。 ```R gene_rank <- computeGeneRank(expression_data, method = "t.test") result <- enrichGSEA(gene_sets, gene_rank) ``` 最后,可以使用enrichplot包中的函数绘制GSEA结果的可视化,例如绘制富集和基因集热。 ```R dotplot(result, showCategory = 20) gene_heatmap(result, top = 10) ``` 通过这些代码,我们可以使用R语言实现GSEA生信分析,从而确定基因集在给定基因表达数据中的富集程度,并可视化展示分析结果。 ### 回答3: GSEA (基因集富集分析) 是一种用于分析生物学实验数据的生物信息学工具,它可以确定在给定条件下,特定基因集中的基因与实验结果相关性的显著性。下面是一个用R语言进行GSEA生信分析的代码示例: 1. 导入所需的R包。 ```R library(clusterProfiler) ``` 2. 导入基因表达数据。 ```R expression_data <- read.table("expression_data.txt", header = TRUE, sep = "\t") ``` 3. 根据实验分组信息创建一个分组向量。 ```R group <- c(rep("Group A", 3), rep("Group B", 3)) ``` 4. 根据基因的符号名称创建一个基因符号向量。 ```R gene_symbols <- c("Gene1", "Gene2", "Gene3", "Gene4", "Gene5", "Gene6") ``` 5. 创建一个基因集对象。 ```R gene_set <- list( GroupA_genes = c("Gene1", "Gene2", "Gene3"), GroupB_genes = c("Gene4", "Gene5", "Gene6") ) ``` 6. 运行GSEA分析。 ```R gsea_result <- gseGO(expression_data, geneSet = gene_set, nPerm = 1000, minGSSize = 3, maxGSSize = 500, pvalueCutoff = 0.05) ``` 7. 查看GSEA结果。 ```R print(gsea_result) ``` 这段代码中,首先导入了clusterProfiler包,它包含了进行GSEA分析所需的函数。然后,基因表达数据被读入到一个名为expression_data的数据框中。接下来创建了一个分组向量,它指定了每个样品所属的实验组。然后,基因符号向量被创建,其中包含了基因的符号名称。根据实验组信息和基因符号,一个基因集对象被创建。最后,调用gseGO函数运行GSEA分析,其中包括参数,如基因集、置换次数、最小/最大基因集大小和显著性阈值。最后,打印GSEA分析的结果。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

码农耕地人~go

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值