2013-04-11 54 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

相關問題