Question: Get SNP information with association test with PLINK
0
gravatar for landscape95
4 weeks ago by
landscape95160
landscape95160 wrote:

Hi everyone,

Does anyone know how to get the SNP information containing CHR ID A1 REF etc with association test on quantitative trait from plink (output file: .qassoc). I done the association test and I want to obtain these information also. In --linear with plink2, I have it, but with plink1.9 association with quantitative trait, I don't have that information (--assoc in plink1.9 runs a lot faster than --linear in plink2 without covariates while --linear with in plink2 gives me memory error:

Error: Out of memory. The --memory flag may be helpful).

I understand that I can use the .bim file, but having integrated information would be very helpful.

Your help is really appreciated!

plink gwas • 149 views
ADD COMMENTlink modified 4 weeks ago by zx87546.8k • written 4 weeks ago by landscape95160
2
gravatar for chrchang523
4 weeks ago by
chrchang5234.8k
United States
chrchang5234.8k wrote:

Can you rerun the analysis with covariates with the latest plink2 build, and post the .log file if it still fails?

ADD COMMENTlink written 4 weeks ago by chrchang5234.8k

Thanks for your prompt reply,

Here is the code and output of plink1.9:

plink --bfile test --assoc --adjust --out output --allow-no-sex --threads 100

PLINK v1.90b4.1 64-bit (30 Mar 2017)
Options in effect:
  --adjust
  --allow-no-sex
  --assoc
  --bfile test
  --out output
  --threads 100

Start time: Wed Feb 20 11:27:04 2019

Random number seed: 1550622424
257682 MB RAM detected; reserving 128841 MB for main workspace.
1133273 variants loaded from .bim file.
230337 people (0 males, 0 females, 230337 ambiguous) loaded from .fam.
Ambiguous sex IDs written to output.nosex .
230337 phenotype values loaded from .fam.
Using 1 thread (no multithreaded calculations invoked).
Before main variant filters, 230337 founders and 0 nonfounders present.
Calculating allele frequencies... done.
Total genotyping rate is 0.992698.
1133273 variants and 230337 people pass filters and QC.
Phenotype data is quantitative.
Writing QT --assoc report to output.qassoc
... done.
--adjust: Genomic inflation est. lambda (based on median chisq) = 1.79545.
--adjust values (1133273 variants) written to
output.qassoc.adjusted .

End time: Wed Feb 20 11:46:25 2019

Here is the code and output with plink2:

plink2 --bfile test --linear --adjust --out output --threads 100

PLINK v2.00a2LM 64-bit Intel (2 Feb 2019)
Options in effect:
  --adjust
  --bfile test
  --glm
  --out output
  --threads 100

Start time: Wed Feb 20 11:35:19 2019

Random number seed: 1550622919
257682 MiB RAM detected; reserving 128841 MiB for main workspace.
Using up to 100 threads (change this with --threads).
230337 samples (0 females, 0 males, 230337 ambiguous; 230337 founders) loaded
from test.fam.
1133273 variants loaded from test.bim.
1 quantitative phenotype loaded (230337 values).
Calculating allele frequencies... done.

Error: Out of memory.  The --memory flag may be helpful.
ADD REPLYlink written 4 weeks ago by landscape95160
2

Hmm, that's a bug, it isn't actually running out of memory, instead some other error is occurring that's getting misreported as "out of memory". Trying to replicate this now.

Update: Replicated the crash, will post a bugfix tonight.

ADD REPLYlink modified 4 weeks ago • written 4 weeks ago by chrchang5234.8k

Thank you so much, Chang.

ADD REPLYlink written 4 weeks ago by landscape95160
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1086 users visited in the last hour