Burden Test for whole exome sequencing data running into errors
Hi,
I am trying to run a burden test on vcf files using PLINK SEQ. When running the program, I run into the error when trying to get the gene based associations. My code is posted below, and the error is posted below. Any help would be appreciated.
I ran gatk on my vcf file to filter out low quality variants. Then fed that file into plink seq
pseq "vcf.file" assoc --phenotype phe1 --mask loc.group=refseq > burdentest.txt
Error: plinkseq warning: convergence problem in tqli() (repeated 3103 times)
• 933 views
•
link
updated
by
Ram
•
written
by
shraimrawan •
•
written
by
shraimrawan •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to perform burden test using plink /seq
written by aiswaryabioinfo •
-
Using admixture to estimate ancestry from a single VCF with no family information
written by dec986
-
Fam file information lost on conversation to VCF file
written by nick.gold •
-
how to merge genetic map of ChrX nonPAR with PAR1 and PAR2 for imputation?
written by nchuang
-
Writing PLINK PED and MAP files... Error: Could not open temporary file.
written by angelaparody •
-
Recoding lgen into Ped using plink
written by njbernstein •
-
Convert plink files to VCF: Reference allele file
written by Ishtiaq Ahmad Khan
-
my VCF file was said to contain no header/meta data when using vcftools
written by CrazyB
-
i-stats for multiple samples
written by potassiumiodide0990
-
Error In Loading Vcfs To Plink Seq
written by Sheila
Please use the formatting bar (especially the
codeoption) to present your post better. You can use backticks for inline code (`text` becomestext), or select a chunk of text and use the highlighted button to format it as a code block. I've done it for you this time.