Re: [Bioc-devel] Changes in AnnotationDbi

2015-06-09 Thread Simon Anders
Hi Martin On 09/06/15 15:35, Martin Morgan wrote: In case you missed it in Marc's reply, and acknowledging that this is different from your suggestion, there is mapIds() for doing this on a single column basis, which is the common use case where one doesn't care too much about multiple mapping

[Bioc-devel] strandless introns with VariantAnnotation::locateVariants()

2015-06-09 Thread Robert Castelo
hi, currently, the annotation of variants in intronic regions by VariantAnnotation and the locateVariants() function does not assign strand to annotations in introns: library(VariantAnnotation) example(locateVariants) loc_all[loc_all$LOCATION == "intron"] GRanges object with 2 ranges and 9 me

Re: [Bioc-devel] Changes in AnnotationDbi

2015-06-09 Thread James W. MacDonald
As select() works currently, the returned keys are in identical order as the input, with extra rows inserted as needed. And any one-to-nothing mapping results in a NA returned. So by definition my (admittedly naive) method is guaranteed to work. But your point is well taken - match() is safer. But

Re: [Bioc-devel] Gene annotation: TxDb vs ENSEMBL/NCBI inconsistency

2015-06-09 Thread Martin Morgan
On 06/08/2015 11:43 PM, Rainer Johannes wrote: dear Robert and Ludwig, the EnsDb packages provide all the gene/transcript etc annotations for all genes defined in the Ensembl database (for a given species and Ensembl release). Except the column/attribute "entrezid" that is stored in the internal

Re: [Bioc-devel] Changes in AnnotationDbi

2015-06-09 Thread Martin Morgan
On 06/09/2015 02:52 AM, Simon Anders wrote: Hi My two cents: On 04/06/15 19:50, James W. MacDonald wrote: In other words, for me it is a common practice to do something like this: fit <- lmFit(eset, design) fit2 <- eBayes(fit) gns <- select(, featureNames(eset), c("ENTREZID","SYMBOL")) gns <-

[Bioc-devel] BioC 2015 Scholarship Application Deadline June 10

2015-06-09 Thread Martin Morgan
See https://bioconductor.org/BioC2015/ for details. -- Computational Biology / Fred Hutchinson Cancer Research Center 1100 Fairview Ave. N. PO Box 19024 Seattle, WA 98109 Location: Arnold Building M1 B861 Phone: (206) 667-2793 ___ Bioc-devel@r-project.

Re: [Bioc-devel] Changes in AnnotationDbi

2015-06-09 Thread Simon Anders
Hi My two cents: On 04/06/15 19:50, James W. MacDonald wrote: In other words, for me it is a common practice to do something like this: fit <- lmFit(eset, design) fit2 <- eBayes(fit) gns <- select(, featureNames(eset), c("ENTREZID","SYMBOL")) gns <- gns[!duplicated(gns[,1]),] fit2$genes <- gns

Re: [Bioc-devel] Gene annotation: TxDb vs ENSEMBL/NCBI inconsistency

2015-06-09 Thread Rainer Johannes
dear Ludwig, On 09 Jun 2015, at 10:46, Ludwig Geistlinger mailto:ludwig.geistlin...@bio.ifi.lmu.de>> wrote: Dear Johannes, Thx for providing the great EnsDb packages! One question: As of now, I am able to choose between TxDb and EnsDb for genomic coordinates of genomic features such as genes

Re: [Bioc-devel] Gene annotation: TxDb vs ENSEMBL/NCBI inconsistency

2015-06-09 Thread Ludwig Geistlinger
Dear Johannes, Thx for providing the great EnsDb packages! One question: As of now, I am able to choose between TxDb and EnsDb for genomic coordinates of genomic features such as genes, transcripts, and exons. For the sequences themselves I need the corresponding BSgenome package. While it is e