questions about qpDstat #123
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.