Illumina甲基化阵列分析软件

methylize的Python项目详细描述


methylize是一个python包,用于分析Illumina甲基化数组的输出。它补充了methylprep和{}。查看ReadTheDocs.

ReadthedocsimageCircleCICodacy BadgeCoverage Status

  1. Overview
  2. Demonstrating differentially methylated probe (DMP) detection (volcano plot) and mapping to chrosomes (manhattan plot)
  3. About BumpHunter

甲基化包装

methylize包包含用于处理来自本地文件的数据的高级api和允许您在运行methylprepmethylcheck之后分析数据的低级功能。为了获得最大的可用性,将methylize与处理过的示例数据(beta值或m值的数据框和包含有关示例的元数据的单独数据框)一起导入Jupyer笔记本。在

Methylize允许您对所有探针进行线性或逻辑回归,并确定DNA发生差异修饰的甲基组中的兴趣点。然后,您可以使用这些回归结果来创建volcano plots和{em1}$manhattan plots。在

曼哈顿地块示例

Manhattan Plot

Manhattan Plot (alternate coloring)

火山图样本

Volcano Plot

可定制:打印尺寸、调色板和截止p值线可由用户设置。 导出:可以导出所有探测统计信息,也可以只将重要探测导出为CSV或python pickled DataFrame。在

安装

pipinstallmethylize

差异甲基化位置/探针(DMP)检测

diff_meth_pos()函数搜索单个差异甲基化位置/探针 (DMPs)通过回归给定条件下每个样品的甲基化M值 根据基因组定位数据。在

表型可以作为

  • 基于字符串的列表
  • 整数二进制数据
  • 数字连续数据
  • (TODO:使用methyprep生成的元数据数据帧作为输入)

当运行logistic回归时,该函数将表型的字符串标签合并为0s和1s。 logistic回归只允许2个表型。线性回归可以有两个以上的表型。在

输入和参数


^{pr2}$

返回

A pandas dataframe of regression statistics with a row for each probe analyzed
and columns listing the individual probe's regression statistics of:
    - regression coefficient
    - lower limit of the coefficient's 95% confidence interval
    - upper limit of the coefficient's 95% confidence interval
    - standard error
    - p-value (phenotype group A vs B - likelihood that the difference is significant for this probe/location)
    - q-value (p-values corrected for multiple testing using the Benjamini-Hochberg FDR method)
    - FDR_QValue: p value, adjusted for multiple comparisons

The rows are sorted by q-value in ascending order to list the most significant
probes first. If q_cutoff is specified, only probes with significant q-values
less than the cutoff will be returned in the dataframe.

如果进度条缺失: 如果在jupyterlab笔记本中看不到进度条,请尝试以下操作: -conda安装-c conda forge nodejs -jupyter labextension安装@jupyter widgets/jupyterlab管理器

欢迎加入QQ群-->: 979659372 Python中文网_新手群

推荐PyPI第三方库


热门话题
在数据库中存储密码的java加密方法   java正则表达式否定整个正则表达式   java为什么要得到这个Hashmap输出?   看不到玻璃鱼爪哇SE 6   类对象包装器中的Java基本数据字段   java从文本文件中读取整数并存储到单独的变量中?(扫描仪)   优化大型Java数据阵列的处理和管理   如何使用Java XML包装类创建对象   java为ExecutorService invokeAll()创建包装器   java如何在Android Studio 1.0.0中设置Facebook SDK?获取SDK位置未找到错误   java在尝试从线程启动动画时调用了FromErrorThreadException   java根据哈希确认文件内容   通过java在neo4j中获取索引值相同的所有节点?   java为什么我的Validare邮政编码(布尔)程序返回false?   java会话自动从servlet/jsp生成,尽管存在以下条件:<%@page session=“false”%>   创建新LANsocket时拒绝java连接   java如何多线程更新由sql代码更新的数据库?   安卓 Java使用类作为集合来添加项   安卓为什么我的清单文件不声明java包?