SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
het:hom distribution in case-control study nans_bn Bioinformatics 0 11-25-2013 06:40 AM
Download Free case study - Next Generation Sequencing ladin General 6 08-21-2012 10:31 AM
GATK error with a VCF including missing genotypes a11msp Bioinformatics 0 04-23-2012 06:33 AM
Case Study using SeqMan NGen at the Next-Gen Sequencing Conference DNASTAR Vendor Forum 0 09-23-2010 12:56 PM

Reply
 
Thread Tools
Old 05-28-2014, 06:24 AM   #1
c_ro87
Member
 
Location: Buenos Aires

Join Date: Feb 2012
Posts: 12
Default Missing genotypes, in case control study using Plink

Hello,

I have a sequencing in ION PGM.

We sequenced 96 barcodes (individuals) and 310 amplicons (chromosomal regions).

32 barcodes are controls and 64 are cases.

We did the variant calling and get 96 VCF files.
We combine them in a single VCF file using GATK. We have 90 different SNPs in the sample.

We convert the single VCF file to plink format (map and ped files) using vcftools.

Now we want to use plink to make a association test,

the ped file looks like this

1 1 0 0 0 2 0 0 0 0 0 0 0 0 0 0
2 2 0 0 0 1 0 0 0 0 0 0 0 0 T C
3 3 0 0 0 2 0 0 0 0 0 0 0 0 T C
4 4 0 0 0 2 0 0 C T 0 0 0 0 T C
5 5 0 0 0 1 0 0 0 0 0 0 0 0 0 0
6 6 0 0 0 2 0 0 0 0 0 0 0 0 0 0
7 7 0 0 0 2 0 0 0 0 0 0 0 0 0 0
8 8 0 0 0 1 0 0 0 0 0 0 0 0 T C
9 9 0 0 0 2 0 0 0 0 0 0 0 0 T C
10 10 0 0 0 2 0 0 0 0 0 0 0 0 T C
11 11 0 0 0 1 0 0 0 0 0 0 0 0 0 0
12 12 0 0 0 2 0 0 C T 0 0 0 0 T C


You can see that there are a lot of missing genotypes, i would like to know what's the standar in this case?

assume that the missing genotypes are references? because most of them probably are, and other could be missing data, but we can't know that, only checking the bam file i guess

If assume the Missing as reference, is there any command in plink to add them automatically?

thanks

Cristian
c_ro87 is offline   Reply With Quote
Reply

Tags
case, control, genotypes, missing, plink

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 07:46 PM.


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