DReichLab / DReichLab/AdmixTools
D-statistics (qpDstat)
Nobody has claimed this yet.
- Dominant language
- C
- Stars
- 235
- Forks
- 76
- PR merge metrics
- No merged PRs in 30d
Description
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....
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start at the qpDstat output path and compare the reported columns with the expected result format in the issue. Determine why version 450 emits only D stat and Z, and document or correct the criteria behind the “best” notation; done means the output behavior and missing fields are explained or fixed.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- c
- Domain
- bioinformatics
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 25/100