草庐IT

统计绘图 | 这个能不能作为【火山图】?

1.火山图?之前瞥到了小明的跟着Cell学作图:R语言ggplot2作图展示差异表达的基因,但没有细看,根据印象做了下面这个图。差异表达结果合并awk'NR==1{print}FNR>1{print}'*vs0.txt>allvs0.txt在R中读取library(tidyverse)df_deg作图,其中不显示padj大于0.5的基因;df_deg的contrast列:1vs0|2vs0|3vs0|4vs0|5vs0,不同的组的差异表达结果;df_deg的change列:Up|Down|Stable,可根据log2FoldChange和padj定义;my_gene是我想突出的基因构成的向量。

统计绘图 | 这个能不能作为【火山图】?

1.火山图?之前瞥到了小明的跟着Cell学作图:R语言ggplot2作图展示差异表达的基因,但没有细看,根据印象做了下面这个图。差异表达结果合并awk'NR==1{print}FNR>1{print}'*vs0.txt>allvs0.txt在R中读取library(tidyverse)df_deg作图,其中不显示padj大于0.5的基因;df_deg的contrast列:1vs0|2vs0|3vs0|4vs0|5vs0,不同的组的差异表达结果;df_deg的change列:Up|Down|Stable,可根据log2FoldChange和padj定义;my_gene是我想突出的基因构成的向量。