merging VCF files
1
0
Entering edit mode
@stephanie-m-gogarten-5121
Last seen 10 days ago
University of Washington
Hi all, Does VariantAnnotation currently have a method to merge VCF objects? I've been looking through the documentation and code and haven't found anything like this. If not, I think it would be a useful feature to add. My use case: I have two VCF files, with the same samples (but in different order in each file). The two files have non-overlapping variants. I would love to have an rbind(VCF, VCF) method; then I could do something like: vcf2 <- vcf2[,colnames(vcf1)] vcf <- rbind(vcf1, vcf2) cbind() would also be useful, for combining files with the same variants but different samples. thanks, Stephanie
VariantAnnotation VariantAnnotation VariantAnnotation VariantAnnotation • 2.5k views
ADD COMMENT
0
Entering edit mode
@valerie-obenchain-4275
Last seen 2.2 years ago
United States
Good suggestion. I'll implement something in the next couple of weeks. Valerie On 01/11/2013 02:22 PM, Stephanie M. Gogarten wrote: > Hi all, > > Does VariantAnnotation currently have a method to merge VCF objects? > I've been looking through the documentation and code and haven't found > anything like this. If not, I think it would be a useful feature to add. > > My use case: I have two VCF files, with the same samples (but in > different order in each file). The two files have non-overlapping > variants. I would love to have an rbind(VCF, VCF) method; then I could > do something like: > > vcf2 <- vcf2[,colnames(vcf1)] > vcf <- rbind(vcf1, vcf2) > > cbind() would also be useful, for combining files with the same variants > but different samples. > > thanks, > Stephanie > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor
ADD COMMENT

Login before adding your answer.

Traffic: 781 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6