D-statistics (qpDstat)

[kafka2184]

Hello In my log result I only have the value "Dstat" and "Z" but no BABA et ABBA ... and "best" notation are on which criteria ? can help me please ? thanks in advance.

Should have : result: Pop1 (W) Pop2 (X) : Pop3 (Y) Pop4 (Z) D-stat Z BABA ABBA #_of_SNPs_that_all_populations_have_data

my out file : qpDstat version: 450 0 "N1" 1 1 "N9" 1 2 "N5" 1 3 "T1" 1 4 "D1" 1 5 "D22" 1 jackknife block size: 0.050 snps: 54984 indivs: 6 number of blocks for jackknife: 68 D stat Z result: "N1" "N9" : "N5" "T1" -0.0998 -1.083 best result: "N1" "N5" : "N9" "T1" -0.3463 -5.492 result: "N1" "T1" : "N9" "N5" -0.2576 -3.761 result: "N1" "N9" : "N5" "D1" -0.1446 -1.591 best result: "N1" "N5" : "N9" "D1" -0.3724 -6.114 result: "N1" "D1" : "N9" "N5" -0.2441 -3.408 etc....

[bumblenick]

You only get the xtra columns with popfilename: You are running with poplistname:

No, there is no logic here...

Nick

On Wed, Nov 8, 2017 at 8:46 AM, kafka2184 notifications@github.com wrote:

Hello In my log result I only have the value "Dstat" and "Z" but no BABA et ABBA ... and "best" notation are on which criteria ? can help me please ? thanks in advance.

Should have : result: Pop1 (W) Pop2 (X) : Pop3 (Y) Pop4 (Z) D-stat Z BABA ABBA

_of_SNPs_that_all_populations_have_data

my out file : qpDstat version: 450 0 "N1" 1 1 "N9" 1 2 "N5" 1 3 "T1" 1 4 "D1" 1 5 "D22" 1 jackknife block size: 0.050 snps: 54984 indivs: 6 number of blocks for jackknife: 68 D stat Z result: "N1" "N9" : "N5" "T1" -0.0998 -1.083 best result: "N1" "N5" : "N9" "T1" -0.3463 -5.492 result: "N1" "T1" : "N9" "N5" -0.2576 -3.761 result: "N1" "N9" : "N5" "D1" -0.1446 -1.591 best result: "N1" "N5" : "N9" "D1" -0.3724 -6.114 result: "N1" "D1" : "N9" "N5" -0.2441 -3.408 etc....

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/DReichLab/AdmixTools/issues/34, or mute the thread https://github.com/notifications/unsubscribe-auth/AQn_h9Itoj2QoAeXGXRor2Wk2Onb-67_ks5s0bC9gaJpZM4QWZu7 .

[kafka2184]

Hello thanks for answer, no I tried poplistname as well and no xtra column

THE INPUT PARAMETERS

PARAMETER NAME: VALUE

DIR: /home/mferrand/ SSS: ALL indivname: /home/mferrand/ALL.ind snpname: /home/mferrand/ALL.snp genotypename: /home/mferrand/ALL.geno popfilename: /home/mferrand/list_D_aN9.txt

qpDstat version: 450

number of quadruples 12 0 "N1" 1 1 "N5" 1 2 "D1" 1 3 "D22" 1 4 "T1" 1 5 "N9" 1 jackknife block size: 0.050 snps: 54984 indivs: 6 number of blocks for jackknife: 68 D stat Z result: "N1" "T1" : "N5" "N9" 0.2576 3.761 result: "N1" "D1" : "N5" "N9" 0.2441 3.408 result: "N1" "D22" : "N5" "N9" -0.0047 -0.675 result: "N5" "T1" : "N1" "N9" -0.0998 -1.083 result: "N5" "D1" : "N1" "N9" -0.1446 -1.591 result: "N5" "D22" : "N1" "N9" -0.3993 -6.860 result: "D1" "D22" : "N1" "N9" -0.2902 -4.143 result: "D1" "D22" : "N5" "N9" -0.2854 -3.621 result: "D1" "T1" : "N1" "N9" 0.0558 0.970 result: "D1" "T1" : "N5" "N9" 0.0265 1.316 result: "D22" "T1" : "N1" "N9" 0.3016 3.991 result: "D22" "T1" : "N5" "N9" 0.2684 3.879

end of qpDstat run