DReichLab / DReichLab/AdmixTools

questions about qpDstat

Open
#123 2 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Dominant language
C
Stars
235
Forks
76
PR merge metrics
No merged PRs in 30d

Description

Dear Professor,
I encountered an issue while running qpDstat.

  • My command line is as follows:

/data/00/user/user187/00.apps/AdmixTools/src/qpDstat -p file.par > allGenomes_chr.D.out

My input files are as follows (showing the first 20 lines of each):

  • qpData/allGenomes_chr.geno:

2002200222222022000
2002222222222222222
2222222222222222202
2222222222222202222
2222222220222222222
2222222222222222022
2222222222222000222
2222222220222222222
2220022000000000022
2222222002022202222
2000022000000000000
2222222002222020222
2002222222222222222
2222222222222222200
2222222222222202222
2222222222222222022
2222202220222222222
2222220002000222222
2000020000000000000
2222222220222222222

  • qpData/allGenomes_chr.snp:

1:15229263 1 0.0 15229263 C T
1:15229267 1 0.0 15229267 T A
1:15229268 1 0.0 15229268 A G
1:15229270 1 0.0 15229270 C A
1:15229271 1 0.0 15229271 G A
1:15229272 1 0.0 15229272 A G
1:15229273 1 0.0 15229273 T A
1:15229275 1 0.0 15229275 T A
1:15229279 1 0.0 15229279 T C
1:15229280 1 0.0 15229280 A G
1:15229281 1 0.0 15229281 T C
1:15229284 1 0.0 15229284 C T
1:19870759 1 0.0 19870759 G A
1:19870760 1 0.0 19870760 G A
1:19870778 1 0.0 19870778 A C
1:19870781 1 0.0 19870781 G A
1:19870782 1 0.0 19870782 G A
1:19870785 1 0.0 19870785 T C
1:19870786 1 0.0 19870786 C T
1:19870787 1 0.0 19870787 C T

  • qpData/Dstat_all.pop:

Papt Con Des Sor
Teu Con Plp Sor
Eur Con Slp Sor
Urot Neu Sygs Sor
Syhn Gal Spn Sor
Sctn Syhn Slp Sor
Eur Mog Des Sor
Sygs Toc Con Sor
Sygs Teu Spn Sor
Urot Syhn Sygs Sor
Eur Spn Urop Sor
Sctn Syhn Eur Sor
Dym Teu Con Sor
Syhn Eur Urop Sor
Urot Mog Con Sor
Dym Spn Urop Sor
Syhn Mog Des Sor
Dym Sygs Urop Sor
Sctn Mog Plp Sor
Syhn Sygs Spn Sor

  • qpData/allGenomes_chr.ind:

Sctn U Sctn
Dym U Dym
Urot U Urot
Toc U Toc
Teu U Teu
Con U Con
Urop U Urop
Syhn U Syhn
Sygs U Sygs
Sor U Sor
Papt U Papt
Eur U Eur
Mog U Mog
Spn U Spn
Neu U Neu
Slp U Slp
Plp U Plp
Gal U Gal
Des U Des

  • The output log and error message are:

cat allGenomes_chr.D.out/data/00/user/user187/00.apps/AdmixTools/src/qpDstat: parameter file: file.par

THE INPUT PARAMETERS##PARAMETER NAME: VALUE

DIR: qpData
SSS: allGenomes_chr
indivname: qpData/allGenomes_chr.ind
snpname: qpData/allGenomes_chr.snp
genotypename: qpData/allGenomes_chr.geno
poplistname: qpData/Dstat_all.pop
printsd: YES
f4mode: NO

qpDstat version: 1152inbreed set NO

qpData/allGenomes_chr.snp: genetic distance set from physical distance

  • And the error log:

cat /data/01/user187/0.data_all/01.Tal/11.diff_tree/01.2.4d_windows/09.D-statistic/qpDstat/qpDstat_all_in_one_pop_file/run.sh.job20251118150335/run.sh.split.sh.1.sl.err
malloc(): invalid size (unsorted)
/slurmState/slurmSpool/slurmd/job1138686/slurm_script: line 10: 2649223 Aborted (core dumped) /data/00/user/user187/00.apps/AdmixTools/src/qpDstat -p file.par > allGenomes_chr.D.out

Could you help me understand the issue I am facing? How can I solve it?
Looking forward to your reply.

Contributor guide

No contributing guide indexed for this repository

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Start by reproducing the crash with file.par and the supplied qpData/allGenomes_chr.* inputs, then inspect the qpDstat input parsing and allocation paths implicated by the malloc error. Done means identifying the malformed input or defect, preventing the abort, and confirming the same command completes successfully.

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

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.