plink2 fails to calculate allele frequency
Hello everyone. I am trying to calculate the equal frequency in groups. code show as below:
plink2 --vcf <VCF path> --freq --pheno <population-file path> --loop-cats population
But I now get an error:
r + plink2-temporary.psam written. 99 samples (0 females, 0 males, 99 ambiguous; 99 founders) loaded from plink2-temporary.psam. Error: plink2-temporary.pvar has a split chromosome. Use --make-pgen + --sort-vars to remedy this.
I tried adding at the end of my code:
plink2 --vcf <VCF path> --freq --pheno <population-file path> --make-pgen --sort-vars
But still the same error occurs. I don't know how to fix this, hope someone can help me take a look, thanks a lot
1 answer
Remove —freq from your command. Wait until —make-pgen + —sort-vars has completed, and then run —freq on the resulting sorted .pgen fileset rather than the unsorted VCF.
• 0 views
•
link
written
by
chrchang523
Log in to answer this question.
More posts like this
-
plink2 giving multiple errors
written by puneet •
-
Can't figure out plink --sample-diff
written by curious
-
VCF to Plink files
written by hi.there •
-
plink2 cannot make bed file
written by dec986
-
Plink Alternative Phenotype File Columns not being Read
written by dtlloyd •
-
Removing extra chromosomes from data
written by nanodano •
-
Error: Failed to open files
written by leslie.ecker •
-
Plink2 error: Skipping --glm regression on phenotype 'PHENO1' since variance inflation factor for c…
written by landscape95
-
PLINK2 association study using covariate regression model
written by tracylee7001 •
-
PLINK error: no --covar values loaded
written by kindlychung •
