PLINK gxe not reading covariate file
I am trying to conduct a GxE in PLINK using a covariate file but PLINK doesn't seem to be recognizing the file. Can someone tell me what I'm doing wrong?
plink --bfile omni --gxe --covar omni.txt
The covariate file omni.txt has the same # of observations as the .fam file and contains FID, IID, Covar1 coded like so...
FID IID Covar1
## ## 1, 2, or -9
The log says that it reads the .bim and .fam files but doesn't even mention the covariate file.
• 4,344 views
•
link
updated
by
Ram
•
written
by
jennavive00 •
•
written
by
jennavive00 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Trouble with PLINK's logistic regression analysis and covariatesTrouble with PLINK's logistic regre…
written by F110152169 •
-
Processing WES VCF for case control GWAS analysis
written by sonsunjirachote •
-
Plink Alternative Phenotype File Columns not being Read
written by dtlloyd •
-
Covariates on non-phenotype fam file?
written by kstafford32 •
-
Plink1.07 generate ERRORs when reading covaries
written by heqianchange •
-
Plink: Losing FAM file data when merging multiple BED files
written by acererak •
-
PLINK unable to read covariate file
written by zichunxu0824 •
-
plink file format - covariate file
written by aj291995 •
-
PLINK error: no --covar values loaded
written by kindlychung •
-
plink analysis including covariates and gxe
written by ddzhangzz
Hi,
Can you post or send me the .log file?
PLINK v1.90b3j 64-bit (2 Apr 2015)
6 arguments: --bfile OMNI --covar omni.txt --gxe
Hostname: BR1
Working directory: W:\AIDGWAS
Start time: Thu Jun 18 11:18:48 2015
Random number seed: 1434651528
12245 MB RAM detected; reserving 6122 MB for main workspace.
608811 variants loaded from .bim file.
5582 people (2927 males, 2655 females) loaded from .fam.
5582 phenotype values loaded from .fam.
End time: Thu Jun 18 11:18:49 2015
This looks like a bug in the command-line parser/dispatcher; I will try to post a fix later today.
Thank you! :)
Okay, I've investigated this; the bug is a missing error message. PLINK should have printed
Error: --gxe requires a scalar phenotype.before exiting. (Today's development build will do that.)This brings us back to your dataset:
--gxedoes not handle case/control phenotypes, so you will need to use another analysis method.I'm a little confused. If it doesn't do case/control studies then what is it used for? It says here https://www.cog-genomics.org/plink2/order that
--gxeis used for quantitative phenotype + case/control covariate association testThe covariate is case/control, but the main phenotype has to be continuous.