2013-04-11 80 views
0

当我使用这个命令的Plink:为什么Hardy-Weinberg测试在GenABEL和Plink中有所不同?

Plink --file filename --hardy 

将创建一个与每个标记不同的p值的新文件filename.hwe,而相形之下,在GenABEL

library(GenABEL) 
... 
snpSummary<-summary(gtdata(filename)) 
hwe_test<-snpSummary$Pexact 

为什么我收到其他p值为相同的标记?

+2

如果您向我们提供一些样本数据,那么它将非常有用,然后准确显示两者之间结果的差异。还要在'R'代码的顶部放置一个'library()'语句,以便我们知道您使用的是哪些软件包,我们可以将其剪切并粘贴到我们的浏览器中。 – nograpes 2013-04-11 17:28:01

回答

2

有不同的HWE测试。例如,卡方合体吻合度与Fisher精确测试vs杂合体过量/缺陷。

我不确定,但您应该检查每个软件的手册,以确认他们首先使用相同的测试。

+0

谢谢KKeenan,该软件不提供任何有关您说的信息!只有命令才能执行。 – 2013-04-12 16:22:11

相关问题