SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
filtering variants to just those in CDS kjmclean Bioinformatics 1 01-30-2016 06:01 AM
question on filtering variants kjaja Bioinformatics 1 08-07-2013 02:49 AM
Filtering out common variants frz Bioinformatics 1 08-07-2013 02:43 AM
filtering the variants kjaja Bioinformatics 1 08-07-2013 02:38 AM

Reply
 
Thread Tools
Old 11-28-2019, 05:01 AM   #1
an_na
Junior Member
 
Location: Poland

Join Date: Mar 2019
Posts: 1
Default filtering variants with R

Hi,

I'm learning data analysis with R, using RStudio.
I have my data in one file. I have variants from 3 tissue types, from 130 individuals. I would like to filter the variants within each individual, so that I can get: 1. unique variants for all 3 tissues, 2. variants common to all tissues, 3. variants common between 2 tissue types.
The file that I'm working on has 118 columns (vcf annotated file), I would like to keep all columns.
I've tried VennDiagram package, but it filters out variants that are present in any other individual, so in the overlapping files I loose some variants.
I understand that I would need to use loop filtering, so that 1. I pick individuals by "ID" column, then sort variants by "variant" column and compare them between "tissue" column to output into all filtered types of variants (unique and common).

Does anyone know how to do that???

Cheers,
A.
an_na is offline   Reply With Quote
Old Yesterday, 05:50 AM   #2
Melissa
Senior Member
 
Location: Switzerland

Join Date: Aug 2008
Posts: 121
Default

To find overlapping variants between two files, you can use bcftools. If you want to generate a simple venn diagram, you can provide lists to http://bioinformatics.psb.ugent.be/webtools/Venn
Melissa is offline   Reply With Quote
Reply

Tags
data analysis help, variant filtering

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 06:27 AM.


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