Would it make sense for an OrganismDb object to have an associated
BSgenome? This association is implicit anyway. Or do we need an even
broader umbrella object? This would make some things easier. For example,
predictCoding(variants, Homo.sapiens). Just nice to have everything in one
bundle. Also,
Let's say I have a vector of gene IDs where some are NA, and are some are
repeated, and I want to map them to gene symbols, where I get NAs for the
NA IDs or IDs without a symbol. What is the best way to do this?
I tried select() but it gave me a table with unique entries; not very
convenient. It
hi Michael,
this souns like if you had a list of variants where you have annotated
their Entrez Gene IDs, which sometimes are NA because those variance do
not overlap any gene and sometimes are repeated Entrez Gene IDs when two
or more of those variants overlap the same gene :)
at least is t
Wow, VariantFiltering is awesome. It's impressive how it integrates all of
the annotation resources so seamlessly. And the shiny app looks very
useful. This should be a model package for where we want to bring
Bioconductor.
You should look into using VRanges instead of GRanges for the return value
hi, thanks of the compliments to the package, i'm happy to hear you
liked it! i must acknowledge that part of the design of the package is
the result of conversations i had with Martin, Marc and specially
Valerie during the review process.
i only got to know about VRanges once the package was
>>I can't install this DESEq2 on the RHEL 5.5 server:
>>
>>* installing *source* package ‘RcppArmadillo’ ...
>>** package ‘RcppArmadillo’ successfully unpacked and MD5 sums
checked
>>* checking LAPACK_LIBS: divide-and-conquer complex SVD unavailable via
R-
>>sup
hi Jarod,
What version of Rcpp and RcppArmadillo is installed?
Mike
On Wed, Jun 18, 2014 at 11:35 AM, jarod...@libero.it wrote:
>
>>>I can't install this DESEq2 on the RHEL 5.5 server:
>>>
>>>* installing *source* package ‘RcppArmadillo’ ...
>>>** package ‘RcppArmadilloâ€
Hi Kristof,
This has been fixed in both the release and devel versions of affy.
Ideally we would change the API so the user has to tell plotLocation()
what array they are using, so the number of rows can be extracted from a
reputable source.
However, the current implementation is pretty frag
Hi Michael,
On 06/18/2014 06:03 AM, Michael Lawrence wrote:
Let's say I have a vector of gene IDs where some are NA, and are some are
repeated, and I want to map them to gene symbols, where I get NAs for the
NA IDs or IDs without a symbol. What is the best way to do this?
I tried select() but i
Thanks for your help
I have the 2.0 when i try to install 3.8 I have the error on version of gcc
on Rhel5.5
>Messaggio originale
>Da: michaelisaiahl...@gmail.com
>Data: 18/06/2014 17.44
>A: "jarod...@libero.it"
>Cc: "bioc-devel@r-project.org"
>Ogg: Re: [Bioc-devel] Problem on Installat
That is a good start. But for convenience, I would favor something that
just returns the vector corresponding to "column" rather than a data.frame.
Thanks,
Michael
On Wed, Jun 18, 2014 at 10:11 AM, Hervé Pagès wrote:
> Hi Michael,
>
>
> On 06/18/2014 06:03 AM, Michael Lawrence wrote:
>
>> L
Seconded, this would be so useful. I still use mget() for heavens sake.
Meanwhile I'm going to try VariantFiltering. Thanks for starting this
conversation Michael.
--t
> On Jun 18, 2014, at 11:18 AM, Michael Lawrence
> wrote:
>
> That is a good start. But for convenience, I would favor so
Hi Michael,
The fact that duplicate keys are not being allowed in OrganismDb objects
is an inconsistency that I will be looking into. That is, if you do the
same thing with orgDb or TranscriptDb objects, duplicates are allowed
and passed along. So I plan to fix this. I have checked in a
pre
are you installing with:
biocLite("DESeq2")
hmm these don't look like Rcpp* version numbers. What do you get with
library(Rcpp)
library(RcppArmadillo)
sessionInfo()
On Wed, Jun 18, 2014 at 1:13 PM, jarod...@libero.it wrote:
> Thanks for your help
> I have the 2.0 when i try to install 3.8 I
We'll try a single arg to REDUCER and see how it goes.
BTW I'm also going to swap out DataFrame for Vector in the rowData.
DataFrame has been more difficult than anticipated (storing names,
subsetting to get ranges out) and doesn't give any clear advantage over
Vector.
Val
On 06/17/2014 0
Hi Peter,
Just added support for "dbSNP" seqlevels style for Human (in
GenomeInfoDb 1.1.9, will become available tomorrow):
library(SNPlocs.Hsapiens.dbSNP.20120608)
myrsids <- c("rs2639606", "rs75264089", "rs73396229", "rs55871206",
"rs10932221", "rs56219727", "rs73709730", "r
On 06/18/2014 05:26 AM, Michael Lawrence wrote:
Would it make sense for an OrganismDb object to have an associated
BSgenome? This association is implicit anyway. Or do we need an even
broader umbrella object? This would make some things easier. For example,
predictCoding(variants, Homo.sapiens).
Sounds nice but we have to be careful. We may need a parameter:
Homo.sapiens("hg19") ... because
broad adoption of new builds may take considerable time.
On Wed, Jun 18, 2014 at 6:13 PM, Hervé Pagès wrote:
> On 06/18/2014 05:26 AM, Michael Lawrence wrote:
>
>> Would it make sense for an Orga
Multiple inheritance might work. May be a good time to think about an
abstract "sequence source", i.e., anything that supports getSeq(). Right
now, we use ANY for that, which is not ideal. Same for "transcript source".
We've already implemented the TxDb API directly on top of biomaRt.
On Wed, Ju
19 matches
Mail list logo