[BioC] distanceToNearest for GRanges?

Cook, Malcolm MEC at stowers.org
Wed Apr 11 18:46:07 CEST 2012


thanks Valerie,

Did you happen to vet the logic of my implementation, or did you just code it lower level from first principals?

Cheers,

~Malcolm

> -----Original Message-----
> From: Valerie Obenchain [mailto:vobencha at fhcrc.org]
> Sent: Wednesday, April 11, 2012 11:13 AM
> To: Janet Young
> Cc: Cook, Malcolm; bioconductor at r-project.org
> Subject: Re: [BioC] distanceToNearest for GRanges?
> 
> This has been added to GenomicRanges 1.9.4. You will also need IRanges
> 1.15.2.
> 
> Valerie
> 
> 
> On 03/27/2012 10:36 AM, Janet Young wrote:
> > Hi,
> >
> > thank you both - very helpful.
> >
> > Janet
> >
> >
> >
> > On Mar 27, 2012, at 6:35 AM, Valerie Obenchain wrote:
> >
> >> Hi Janet, Malcom,
> >>
> >> Yes, we can add distanceToNearest for GRanges. This will probably
> happen after the release (i.e., next week). I'll post back here when it's
> complete.
> >>
> >> Valerie
> >>
> >>
> >> On 03/26/12 16:26, Cook, Malcolm wrote:
> >>> Janet,
> >>>
> >>> In the mean time....
> >>>
> >>> distanceToNearest.GRanges<- function
> (query,subject=query,...,.AS=DataFrame) {
> >>>   ## PURPOSE: an implementation of distanceToNearest for GRanges in
> >>>   ## terms of 'nearest' and 'distance' whose return value is shaped
> >>>   ## like that of IRange's 'nearest', by default, but you can ask for
> >>>   ## the result .AS=data.table or .AS=data.frame if you prefer.
> >>>    if(missing(subject)) {
> >>>      nearestSubjectIndex<-nearest(query,...)
> >>>    } else {
> >>>      nearestSubjectIndex<-nearest(query,subject,...)
> >>>    }
> >>>    queryHasNearest<-! is.na(nearestSubjectIndex)
> >>>    queryIndex<-which(queryHasNearest)
> >>>    queryHavingNearest<-query[queryIndex]
> >>>    nearestSubjectIndex<-nearestSubjectIndex[queryIndex]
> >>>    nearestSubject<-subject[nearestSubjectIndex]
> >>>    d<-distance(queryHavingNearest,nearestSubject)
> >>>    .AS(query=queryIndex
> >>>               ,subject=nearestSubjectIndex
> >>>               ,distance=d)
> >>> }
> >>>
> >>> You've asked so many questions in this forum whose answer I made use
> of that I must owe you this tidbit....
> >>>
> >>> Lurkers, and fellow travelers, I welcome corrections and stylistic
> suggestions.
> >>>
> >>> Cheers,
> >>>
> >>> Malcolm Cook
> >>>
> >>>
> >>> ________________________________________
> >>> From: bioconductor-bounces at r-project.org [bioconductor-bounces at r-
> project.org] On Behalf Of Janet Young [jayoung at fhcrc.org]
> >>> Sent: Monday, March 26, 2012 4:42 PM
> >>> To: bioconductor at r-project.org
> >>> Subject: [BioC] distanceToNearest for GRanges?
> >>>
> >>> Hi there,
> >>>
> >>> Any plans to implement distanceToNearest for GRanges objects?   That
> would be really useful for me.
> >>>
> >>> thanks,
> >>>
> >>> Janet
> >>>
> >>>
> >>> -------------------------------------------------------------------
> >>>
> >>> Dr. Janet Young
> >>>
> >>> Tapscott and Malik labs
> >>>
> >>> Fred Hutchinson Cancer Research Center
> >>> 1100 Fairview Avenue N., C3-168,
> >>> P.O. Box 19024, Seattle, WA 98109-1024, USA.
> >>>
> >>> tel: (206) 667 1471 fax: (206) 667 6524
> >>> email: jayoung  ...at...  fhcrc.org
> >>>
> >>>
> >>> -------------------------------------------------------------------
> >>>
> >>> _______________________________________________
> >>> 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
> >>> _______________________________________________
> >>> 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



More information about the Bioconductor mailing list