Burden tests in plinkseq (v0.08) running exponentially slower with more variants?
0
0
Entering edit mode
9.6 years ago
vlaufer ▴ 290

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 • 2.3k views
ADD COMMENT

Login before adding your answer.

Traffic: 2732 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6