SuperAPI AI 工具集

volcano-plot-script

编程与开发

查看项目

更新于 10-15 16:45

技能介绍

用途:从差异表达基因(DEG)分析结果生成火山图(volcano plot)的 R/Python 代码,用于展示基因表达数据中统计显著性(p 值)与变化幅度(倍数变化)的关系,适合生物信息学分析与论文插图。触发场景包括 RNA-seq DEG 结果可视化、p 值与 fold-change 散点图、寻找显著上调/下调基因、标注目标基因、多条件比较等。

主要能力:提供打包的可执行入口 scripts/main.py,可按参数生成可定制的 R 或 Python 脚本及发表级图表(PNG/PDF/SVG);支持自定义 log2FC 与 p 值阈值、基因注释与高亮、按颜色区分上调/下调/不显著基因、输出带标签的显著基因列表。参考材料在 references/,可复用资产包括 assets/example_volcano.R。

使用方式:输入需包含基因标识(symbol 或 ENSEMBL ID)、log2 倍数变化值、校正或原始 p 值,可选基因注释与通路信息。常用参数:--input(DEG 结果 CSV/TSV,必填)、--output(默认 volcano_plot.png)、--log2fc-col(默认 log2FoldChange)、--pvalue-col(默认 padj)、--gene-col(默认 gene)、--log2fc-thresh(默认 1.0)、--pvalue-thresh(默认 0.05)、--label-genes、--top-n(默认 10)及上/下/非显著颜色参数。示例:python scripts/main.py --input deg_results.csv --output volcano_plot.png。运行前可用 python -m py_compile scripts/main.py 校验脚本,python scripts/main.py --help 查看帮助。Python 依赖 pandas、matplotlib、seaborn、numpy;R 依赖 ggplot2、dplyr、ggrepel。

适用场景:RNA-seq DEG 可视化、显著基因识别、通路或标记基因标注、论文级图表制作、多实验条件比较。

限制与风险:技术难度中等,需理解 fold change、p 值、FDR 及绘图库;脚本本地执行,涉及代码执行与文件读写风险(中),无外部 API 调用;要求输入路径与输出目录受工作空间限制,输入需校验;当前生命周期为 Draft,已知问题为无,计划增加 Plotly 交互图、多比较组支持与通路富集集成。若输入缺失或任务超出文档范围,应明确指出缺失字段并停止或请求最小补充信息,不臆造文件、引用或执行结果。价格、收费、具体平台与效果数据暂未提供。