This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Burden tests in plinkseq (v0.08) running exponentially slower with more variants?

Hello, I am implementing tests like the following in plinkseq (v0.08)

pseq LOCUS assoc --tests calpha --perm 10000 --phenotype RA --mask loc.group=refseq maf=0.00:0.01

These run fine, and fairly quickly. However, if I change this latter test to:

pseq LOCUS assoc --tests calpha --perm 10000 --phenotype RA --mask loc.group=refseq maf=0.00:0.05

The program runs very slowly. I am not sure if it is hanging or not; a value of 0.03 runs slowly but does eventually complete.

Does anyone know the Big-O notation for this? Or, more empirically, how many variants could I run before I can expect a slowdown. Or am I likely running out of RAM, and if I increased this it would run? Any feedback is welcome.

Thanks.

plinkseq genetics association statistics

0 answers

No answers yet.

Log in to answer this question.