方法与软件
方法:各步骤的方法写在对应页面的“?”里,这里集中列出阈值和版本。
结果:阈值全部来自配置文件。
阈值
| 项目 | 标准 |
|---|---|
| 差异基因 | 校正后 P<0.05,|log2FC|≥1(服务商判定) |
| 差异代谢物 | P<0.05,FC≥2,VIP≥1(Welch t 检验,P 未校正) |
| 通路富集 | 超几何检验 + Fisher 合并 + BH,FDR<0.05 |
| 基因-代谢物相关 | Pearson,|r|≥0.8,BH FDR<0.05 |
| 基因纳入 | 至少 3 个样本 FPKM≥1,或任一比较差异 |
| 通路范围 | KEGG ko 层级,去总览图及人类疾病、器官系统、药物开发类;只保留 ath、osa、aof、vvi、gmx、sly、zma 至少一种具有的通路;排除甜菜红素 ko00965 |
软件版本
| software | version |
|---|---|
| R | 4.4.3 |
| data.table | 1.18.6.1 |
| ggplot2 | 4.0.3 |
| WGCNA | 1.74 |
| Mfuzz | 2.66.0 |
| mixOmics | 6.30.0 |
| OmicsPLS | 2.1.0 |
| fgsea | 1.32.2 |
| vegan | 2.7.5 |
| ade4 | 1.7.24 |
| ropls | 1.38.0 |
| MOFA2 | 1.16.0 |
| GENIE3 | 1.28.0 |
| igraph | 2.3.4 |
| ggraph | 2.2.2 |
| circlize | 0.4.18 |
| pheatmap | 1.0.13 |
| limma | 3.62.1 |
| Python | 3.11.18 |
| pandas | 3.1.1 |
| Pillow | 13.0.0 |
| KEGG (REST 快照) | pathway 587 2026/11/05;ko 28,565 2026/11/06 |