SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
VCF file missing GT field - How to fix? user 31888 Bioinformatics 0 04-13-2016 02:46 PM
Samtools VCF file: NS field value missing/not printing probean Bioinformatics 0 02-07-2014 10:04 AM
rare variant missing in the vcf file crang Bioinformatics 10 04-04-2013 08:27 AM
VCF file, missing AA, AC and AN tags phatjoe Bioinformatics 0 01-15-2012 04:17 PM
How to calculate proportion of reads with each base at every reference position lindseyjane Bioinformatics 7 01-19-2010 07:43 AM

Reply
 
Thread Tools
Old 08-18-2016, 03:19 AM   #1
LeonDK
Member
 
Location: Denmark

Join Date: Sep 2014
Posts: 69
Default Calculate the proportion of missing data per sample from VCF file

I have a VCFv4.2 (Sanger-imputation-service genotype data) file and I need to calculate the proportion of missing data per sample.

Each sample is encoded like so:
Code:
GT:ADS:DS:GP
0|0:0.25,0.15:0.4:0.6375,0.325,0.0375
I am uncertain as to how 'missing' is identified - Any suggestions?
LeonDK is offline   Reply With Quote
Old 08-18-2016, 04:00 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Missing data will have ".|." or "./." as the genotype. However, since this is imputed I'm not sure how much there will be in the way of missing values (they should have been largely imputed).
dpryan is offline   Reply With Quote
Old 08-18-2016, 05:11 AM   #3
LeonDK
Member
 
Location: Denmark

Join Date: Sep 2014
Posts: 69
Default

Quote:
Originally Posted by dpryan View Post
Missing data will have ".|." or "./." as the genotype. However, since this is imputed I'm not sure how much there will be in the way of missing values (they should have been largely imputed).
Running the following on the vcf files
Code:
for f in `ls -1`; do gunzip -c $f | grep -c '\.|\.'; done
Yields nothing but zeros, so it looks like you are correct.
LeonDK is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:10 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO