Trouble with PLINK's logistic regression analysis and covariatesTrouble with PLINK's logistic regression analysis and covariates

Hi,

I am new to PLINK and I am having trouble using one of the commands. I am trying to perform logistic regression analysis on several Ped and map files. Here is the command I put in:

plink --file TWBR10907-07_TWB2_QC --covar TWB2.0_2.6_2.txt --logistic --covar-name FID,IID,SEX,PARKISON --out logistic_results

TWB2.0_2.6_2.txt:
FID IID SEX PARKISON
TV2N_20203805G  TV2N_20203805G  1   0
TV23_20182701G  TV23_20182701G  1   0
TV29_20181602E  TV29_20181602E  1   0
TV23_20182809B  TV23_20182809B  2   0
TV23_20180107G  TV23_20180107G  2   0
TV2N_20203806G  TV2N_20203806G  2   0
TV23_20180108G  TV23_20180108G  1   0
TV29_20181603F  TV29_20181603F  1   0
TV21_20180301B  TV21_20180301B  1   0

I get back the following message and error:

283083 variants loaded from .bim file.
81236 people (25584 males, 55652 females) loaded from .fam.
81236 phenotype values loaded from .fam.
Using 1 thread (no multithreaded calculations invoked).
Error: Missing --covar-name token in --covar file header line.

I've searched online and I can't find anything on this error. Can anyone please help me out?

Exia

plink gwas snp

• 1,084 views • link

written 2.2 years ago by F110152169 •

Try: --covar-name SEX,PARKISON

FID and IID are Family ID and Individual ID, they are not covariates.

Side note, if PARKISON means "Parkinson's disease", then it is a typo.

• 0 views • link

written 2.2 years ago by zx8754

It's Parkinson's disease. How do I get the right order?

• 0 views • link

written 2.2 years ago by F110152169 •

0 answers

No answers yet.

More posts like this

Processing WES VCF for case control GWAS analysis
written 3.2 years ago by sonsunjirachote •

Dear community I'm currently trying to try a GWAS for case control study using PLINK using whole exome VCF files. I understand that I need …
plink logistic regression error
written 3.9 years ago by dec986

I'm running the following command with Plink 1.90b6.26: plink --bfile European.imputed.merged.cleanedx --logistic --pheno ../European.pheno.plink.tsv --pheno-name event --covar ../European.pheno.plink.tsv --covar-name age,Sex,PC1,PC2 –-ci 0.95 --out t2d_basic_model.txt where …
Adding additional covariates to plink logistic regression model on HLA returns p-vals etc as NA
written 4.3 years ago by 4galaxy77

I am running a logistic regression association test on the imputed HLA genotypes, using TB as a phenotype, using the HLA-tapas wrapper around plink v1.90. …
Plink Alternative Phenotype File Columns not being Read
written 4.6 years ago by dtlloyd •

Hi, I have a plink alternative phenotype file with the following format: ``` FID IID Phenotype 1 2 1 1 3 0 .. ``` As …
Plink1.07 generate ERRORs when reading covaries
written 5.4 years ago by heqianchange •

Hi, everyone. I am using Plink1.07 to do GWAS. When I go to logistic association analysis with covaries will generate errors. I make the covariate …
How to resolve "Error: All remaining samples for --glm phenotype 'pheno' are cases" in Plink?
written 6.2 years ago by anamaria

Hello, I was running GWAS via: plink2 --threads 24 --pfile temporary17 --recover-var-ids NVCFchr17.vcf.gz 'force' --pheno pheno_hes_nephropathy_hba1c.txt --pheno-name pheno --covar pheno_hes_nephropathy_hba1c.txt --covar-name sex,age,PC1,PC2,PC3,PC4,PC5,PC6,PC7,PC8,PC9,PC10,TD,array,HBA1C --glm genotypic cols=+beta,+err …
How to update Family ID in bed file
written 10.0 years ago by rajneeshpaliwal •

Hi, I am new user for plink. I have a PED file with Individual IDs for both Family IDs (FID) and individual ID. But now …
Trouble with PLINK's logistic regression analysis and covariates
written 10.5 years ago by Diego.Morales •

Hi, I am new to PLINK and I am having trouble using one of the commands. I am trying to perform logistic regression analysis on …
PLINK2 association study using covariate regression model
written 11.4 years ago by tracylee7001 •

Hello all, I'm using PLINK2 to do association study in a case-control trait. I have 4 covariates. Here is my command: ``` plink2 \ --vcf …
PLINK error: no --covar values loaded
written 11.5 years ago by kindlychung •

Here is the command that I used: ``` plink \ --allow-no-sex \ --bfile ./collr_task_Curliness-sex_age_s400_p1.00e-01/RS123_3069obs_shift_0000 \ --covar curliness_RS_withIID.tab \ --covar-name sex,age \ --linear hide-covar \ --out …