Hi Katsuhito,
I've cc'd your message to the Bioconductor mailing list. We have two
lists, one for general questions and the other for bug reports/feature
requests. Please post future questions to one of these lists instead of
sending them to a single person. The lists reach a wider audience and
others can chime in with their responses/experience. You can find info
about the mailing lists here,
http://www.bioconductor.org/help/mailing-list/
This bug was fixed in the devel version a few weeks ago. I've pushed the
fix to release, version 1.8.6, which will be available via biocLite()
tomorrow ~10am PST.
Valerie
On 11/20/2013 01:05 PM, Yasuno, Katsuhito wrote:
Dear Valerie,
When I read VCF files, I realized that elements of PL field (Number = G)
always have
extra slots with NA. In fact, when writing the corresponding vcf object
to a VCF file, it causes
NA's to be printed. Can your package handle Number = G fields? Thanks
ID10-1 ID100-1 ID100-35 ID100-36 ID101-1
21:9908330_A/G Integer,4 Integer,4 Integer,4 Integer,4 Integer,4 <- 1
alternate allele
21:9908348_C/T Integer,4 Integer,4 Integer,4 Integer,4 Integer,4
21:9908370_G/A Integer,4 Integer,4 Integer,4 Integer,4 Integer,4
21:9909124_C/T Integer,4 Integer,4 Integer,4 Integer,4 Integer,4
21:10862573_T/C Integer,4 Integer,4 Integer,4 Integer,4 Integer,4
21:10862598_C/A Integer,9 Integer,9 Integer,9 Integer,9 Integer,9 <- 2
alternate alleles
> sessionInfo()
R version 3.0.2 (2013-09-25)
Platform: x86_64-unknown-linux-gnu (64-bit)
locale:
[1] C
attached base packages:
[1] parallel stats graphics grDevices utils datasets
methods base
other attached packages:
[1] VariantAnnotation_1.8.5 Rsamtools_1.14.1
Biostrings_2.30.1 GenomicRanges_1.14.3 XVector_0.2.0
[6] IRanges_1.20.5 BiocGenerics_0.8.0
loaded via a namespace (and not attached):
[1] AnnotationDbi_1.24.0 BSgenome_1.30.0
Biobase_2.22.0 DBI_0.2-7 GenomicFeatures_1.14.2
[6] RCurl_1.95-4.1 RSQLite_0.11.4
XML_3.98-1.1 biomaRt_2.18.0 bitops_1.0-6
[11] rtracklayer_1.22.0 stats4_3.0.2
tools_3.0.2 zlibbioc_1.8.0
Best,
Katsuhito
############################
Katsuhito Yasuno, PhD
Department of Neurosurgery
Yale University School of Medicine
1 Gilbert Street, TAC S327
New Haven, CT 06510
katsuhito.yas...@yale.edu
(203) 785 3056
############################
_______________________________________________
Bioc-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/bioc-devel