我正在分析基因组数据,寻找一种快速从噪声中解析出重要峰值的方法,以获得一系列跨染色体的统计指标(例如Tajima's D)。
有没有人知道实现Garvilov和Adler (2011)中描述的峰值检测协议的脚本,多次测试本地最大值以检测1D中的峰值。安·斯塔特39(6) 3290-3319,doi: 10.1214/1111-AOS943
发布于 2013-11-26 08:18:19
我自己也遇到了类似的问题。我不熟悉你引用的论文,但我认为GWAS文献传统上定义了关于中性期望下的零分布的“显著性”,这是通过排列、自举和重新计算数据生成的。我很抱歉“要求澄清”,但是你能告诉我更多关于你的方法和理由吗?
https://stackoverflow.com/questions/19326976
复制相似问题