用途:为差异表达(DE)分析生成标准图件面板,分为诊断图与结果图两类,把图形当作底层统计模型的诊断依据。主要能力包括:离散度图 plotDispEsts(dds)(DESeq2)与 plotBCV(y)(edgeR);p 值直方图(用 ggplot2 绘制);基于 VST/rlog 的 PCA 与样本距离热图;MA 图 plotMA(res) 与 plotMD(qlf);火山图(ggplot2 包装或 EnhancedVolcano);Top 差异基因热图;单基因 plotCounts。
必要用法要点:先看离散度图与 p 值直方图,再信任结果;PCA 必须先做方差稳定(vst 或 rlog),不可用原始 counts,QC 用 blind=TRUE,结果图用 blind=FALSE;n<30 且文库大小差异超 4 倍时用 rlog,n>30 用 vst。火山图横轴用收缩后 LFC(lfcShrink 的 apeglm),纵轴为未收缩 Wald p 值,坐标轴需注明;标签超过 10 个基因时设 max.overlaps=Inf 或 options(ggrepel.max.overlaps = Inf),否则默认 10 会静默丢标签。热图看模式用 scale='row'(按基因 z 分数),QC 热图看绝对水平用 scale='none'。选高变基因建议用 matrixStats::rowMads 而非 rowVars。超过 3 个基因集比较用 UpSet 而非 Venn。
限制:资料仅给出参考版本范围(DESeq2 1.42+、edgeR 4.0+、limma 3.58+、ggplot2 3.5+、pheatmap 1.0+、RColorBrewer 1.1+、ggrepel 0.9+、EnhancedVolcano 1.20+、matrixStats 1.2+),使用前需用 packageVersion 核对版本并按实际 API 调整示例;EnhancedVolcano 1.10-1.11 参数名不同。n=3 时结果应视为探索性。其余运行环境、收费与完整安装方式暂未提供。