... ave() will give the sequential numbers, so something like
with(mydf, paste0(ID, ave( ID, ID, FUN=seq_along )))
should do it.
I'd like to be able to specify the starting 'centers' for dmn().
IIUC DirichletMultinomial::dmn(count, k) will initialize the EM
using a kmeans heuristic for selecting the starting point. Replicate
the same data can yield stark differences in the result.
I have a datas ...
... Michael Lawrence writes:
Thanks for your reply and for GenomicRanges.
Please see inline comments below.
> On Tue, Jan 7, 2014 at 9:30 AM, Charles Berry
> > Vince S. Buffalo ...> writes:
> > >
> > > Hi Bioconductor folks,
> > ...
... Vince S. Buffalo writes:
> Hi Bioconductor folks,
> I'm trying to create some GRanges summaries, but I think I may be
> an obvious solution. I have fixed-width windows as a GRanges object,
> for each window/row I'd like to add a metadata column that indicates
... Peter Hickey writes:
> Hi all,
> Suppose I've created a Hits object by
> all_hits <- findOverlaps(query, subject, select = 'all', type =
> and now I want to get the first (or last) hit for each query. Of
gr1 <- GRanges(rep(letters,each=10),
... Ludo Pagie writes:
> Hi all,
> I don't understand what's going on here.I'm trying to compute the
> coverage of a few ranges. I specify a weight for these ranges and
> some positions I get a coverage even though I know these positions
> only covered by ranges which ...
... Valerie Obenchain writes:
> Hi Andrew,
> I can't reproduce this problem. Can you successfully save the files
> different order, i.e., is it just the addition of 'introns' that is
> giving you trouble?
> I assume these are custom annotation GRanges that did not co ...
... Dario Strbenac writes:
> Is there a way to move BSgenome.Hsapiens.UCSC.hg19 from an R 2.15
to an R 3.0 library without
> downloading it all again ? I installed it as a Windows binary.
If so, I w ...