Question: PLINK - Sample(s) present with no genotype data
0
gravatar for liorglic
3 months ago by
liorglic310
liorglic310 wrote:

Hello,
I am running PLINK v1.90b6.12 , trying to perform a PCA analysis on a vcf file.
I run the following command:

plink --vcf all.maf005.rehead.vcf --pca --allow-extra-chr

and getting the ouput:

PLINK v1.90b6.12 64-bit (28 Oct 2019)          www.cog-genomics.org/plink/1.9/
(C) 2005-2019 Shaun Purcell, Christopher Chang   GNU General Public License v3
Logging to plink.log.
Options in effect:
  --allow-extra-chr
  --pca
  --vcf all.maf005.rehead.vcf

128770 MB RAM detected; reserving 64385 MB for main workspace.
--vcf: plink-temporary.bed + plink-temporary.bim + plink-temporary.fam written.
17 variants loaded from .bim file.
3820 people (0 males, 0 females, 3820 ambiguous) loaded from .fam.
Ambiguous sex IDs written to plink.nosex .
Using up to 15 threads (change this with --threads).
Before main variant filters, 3820 founders and 0 nonfounders present.
Calculating allele frequencies... done.
Total genotyping rate is 0.198984.
17 variants and 3820 people pass filters and QC.
Note: No phenotypes present.
Relationship matrix calculation complete.
Error: Sample(s) present with no genotype data.  Use --mind to filter them out.

However, there is no missing data in my vcf - all genotypes are "./." or "1/1".
If I run:

plink --vcf all.maf005.rehead.vcf --pca --allow-extra-chr --mind

the ouput is:

PLINK v1.90b6.12 64-bit (28 Oct 2019)          www.cog-genomics.org/plink/1.9/
(C) 2005-2019 Shaun Purcell, Christopher Chang   GNU General Public License v3
Logging to plink.log.
Options in effect:
  --allow-extra-chr
  --mind
  --pca
  --vcf all.maf005.rehead.vcf

128770 MB RAM detected; reserving 64385 MB for main workspace.
--vcf: plink-temporary.bed + plink-temporary.bim + plink-temporary.fam written.
17 variants loaded from .bim file.
3820 people (0 males, 0 females, 3820 ambiguous) loaded from .fam.
Ambiguous sex IDs written to plink.nosex .
Error: All people removed due to missing genotype data (--mind).
IDs written to plink.irem .

Can anyone help figure out this issue?
Thanks!

plink • 217 views
ADD COMMENTlink modified 3 months ago by chrchang5237.1k • written 3 months ago by liorglic310
2
gravatar for chrchang523
3 months ago by
chrchang5237.1k
United States
chrchang5237.1k wrote:

./. encodes a missing genotype.

ADD COMMENTlink written 3 months ago by chrchang5237.1k

Stupid me - I thought ./. means ref. Now I see it should be 0/0. The vcf was created with paftools call - apparently it puts and unknown genotype whenever it doesn't find a variant. Thanks anyway.

ADD REPLYlink written 3 months ago by liorglic310
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: 1478 users visited in the last hour