用于微阵列数据的稳健多阵列平均值的python脚本

时间:2010-02-18 23:55:40

标签: python arrays bioinformatics dna-sequence

我试过谷歌没有运气。我已经看到一些对使用python但没有代码完成的稳健多阵列平均的弱引用。我对重新发明轮子不是那么感兴趣。关于python模块的任何建议,脚本......

如果我能找到一个很好的解释或算法的例子,我会写一个python实现来分享。

如果你不确定我在说什么,你可以看看这个,虽然这不是定义。 http://www.mathworks.com/access/helpdesk/help/toolbox/bioinfo/ref/gcrma.html

3 个答案:

答案 0 :(得分:3)

尽管这个问题是在六年多前提出的,当我最近问自己同样的问题时,我仍然找不到一个很好的Python实现RMA。因此,我对该方法进行了大量研究,并自己写了一篇:

我还写了一篇关于它的文章,其中涉及一些技术细节:

pyAffy: An efficient Python/Cython implementation of the RMA method for processing raw data from Affymetrix expression microarrays

答案 1 :(得分:2)

你可能想考虑为R使用python接口,它有gcrma包。 RPy是一个模块,允许您使用系统上安装的所有R模块。 Bioconductor有一个gcrma module用于R.我找不到任何用于Python的模块。

答案 2 :(得分:0)

我是一名从事基因组学工作的生物学家,为5-6个微阵列项目做出了贡献。我也非常精通Python,可以在R中做得很好。

Python是一种很棒的语言,但截至目前,我发现没有任何模块可以轻松满足您的要求。另一方面,R为微阵列分析提供了大量的包装。

鉴于微阵列数据分析是如何通过医疗使用标准微阵列芯片/技术驱动的,生物学研究必须使用定制开发(主要是R)包进行分析。不幸的是,我认为这种情况在更成熟的编程语言中为微阵列数据分析留下了一个不那么有趣的位置,而这些语言并没有被研究生群体所认识。

希望你找到一个python替代品,在这种情况下,请告诉我们!

干杯