又现Science，这种图片的含金量还在上升~

文章预览

生信益站，一点就有益 ！ 祝友友们天天开心，月月发 CNS~ 设我为🌟星标，每天就能第一时间看到推送啦~ 这些图来自《Science》 DOI: 10.1126/science.adn8105 DOI: 10.1126/science.adn8105 仔细看了下，站长之前画过类似的图 DOI: 10.1126/science.adn8105 ❝ 也就是 抛弃火山图的第5天 （篇）—— 差异基因等级条形图 ：展示 FDR、log2FC、以及感兴趣基因的label 。 感兴趣的可以了解下~ 数据读取 所有基因差异表格： df "WT-vs-YY.DEseq2.tsv",sep= "\t" ,head= T ) 显著（FDR < 0.05 & |log2FC|>1）差异基因表格： df.sig "WT-vs-YY.DEseq2.Sig.tsv",sep= "\t" ,head= T ) 基础图形 library (ggplot2) library (dplyr) df %>%       mutate(FDR_Region = case_when(         FDR 0.001 ~  " < 0.001" ,         FDR 0.05 ~  " < 0.05" ,         FDR >=  0.05  ~  ">=0.05"     ))  %>%      ggplot(aes(x = reorder(GeneID,logFC),           ………………………………