DReichLab / DReichLab/AdmixTools

Question about “duplicate in list” error and memory usage when running qpDstat

Open
#121 1 comment 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 Author,

We are currently running qpDstat from the AdmixTools package using the following command:

/data/01/user157/software/AdmixTools/bin/qpDstat -p parDstat_alltest > test.out

Our parameter file is as follows:

indivname: snp.eigenstrat.ind
snpname: snp.eigenstrat.lower.chr.snp
genotypename: snp.eigenstrat.chr.geno
poplistname: test.popfilename.txt
printsd: YES
f4mode: NO
blgSize: .01

The .ind, .geno, and .snp files are as follows:

Image Image Image

We have 14 species within a single genus, and for each triplet of ingroup species, we used a species from another genus (M_un) as the outgroup.
In the popfilename.txt, we therefore listed all possible three-species combinations (each with the same outgroup) and included all three permutations for each triplet, as shown in the attached example (see screenshot).

Image

However, during execution we encountered the following error:

fatalx:
(loadlist) duplicate in list: A_ca
/slurmState/slurmSpool/slurmd/job1133823/slurm_script: line 10: 3990729 Aborted (core dumped)

Could you please let us know what might be causing this “duplicate in list” error?
We checked our population list file, but we are sure that no duplicate population names is within a single line.

In addition, our dataset includes 278,518,514 SNPs, which seems to require very large memory.
We currently allocate 2 CPUs and 250 GB of memory, but in some runs qpDstat still exceeds memory limits.
Do you have any recommendations on the typical memory requirements or efficient settings for datasets of this size (e.g., adjusting blgSize or using a subset of SNPs)?

Thank you very much for your time and for maintaining such an important tool for population genetic analysis.

Best regards,
Na Wan

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 with the qpDstat entry point and the parDstat_alltest and popfilename.txt inputs, then inspect how duplicate population names and large SNP datasets are handled. Compare the reported duplicate error and memory failures with the tool's documented parameters, including blgSize; done means identifying the cause and documenting tested memory-efficient settings.

Written by the indexing model from the issue text.

Assessment

Tech stack
c
Domain
bioinformatics, data
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.