用途：从拟合好的 DESeq2 或 edgeR 模型提取差异表达估计与 p 值，处理缺失的 padj，按设计选择合适的 FDR 控制，生成下游工具所需的表格或排序列表。
主要能力：用 results() 或 topTags() 显式指定对比并转成数据框；区分 padj=NA 的三种成因（独立过滤、Cook's 距离离群、组内全零）并给出对应处理（independentFiltering=FALSE、cooksCutoff=FALSE、预处理过滤）；支持 BH、Storey q-value、IHW、ashr lfsr/svalue、BY、Holm/Bonferroni 等校正方法；用 TREAT 或 lfcThreshold= 控制“倍数变化”假设的 FDR；绘制 p 值直方图诊断模型问题；用 filter/arrange 筛选与排序；用 AnnotationDbi+org.db、biomaRt 或 mygene 注释，原核用 GFF；生成 GSEA preranked 排序列表；为 ORA 提供正确 universe。
用法要点：resultsNames(dds) 后 results(dds, name=..., alpha=0.05)，lfcShrink(type='apeglm')；edgeR 用 topTags(sort.by='none') 保持行序。注意列名差异（DESeq2 log2FoldChange/padj、edgeR logFC/FDR、limma logFC/adj.P.Val）。GSEA 用 stat 或收缩 LFC，不可用未收缩 LFC；ORA 的 universe 必须是测试过的基因集。
限制与风险：na.omit 会丢弃真实信号；事后按 padj<0.05 且 |LFC|>tau 过滤不能控制“变化幅度”的 FDR；IHW 在测试数少或协变量与处理相关时无益；lfsr 与 padj 不可互换；旧 HGNC 符号会静默丢基因，应连接稳定的 Ensembl/Entrez ID；n<=3/组 应视为探索性。资料未提供具体性能、收费或平台信息。