[Bioc-devel] Reduce Gap Width

Martin Morgan mtmorgan at fhcrc.org
Tue Aug 9 16:13:51 CEST 2011


Hi Dario --

On 08/09/2011 12:00 AM, Dario Strbenac wrote:
> Could it be made possible to neatly reduce all GRanges on a chromosome ?
>
> reduce(GRanges("chr1", IRanges(c(1000, 9000), width = 1000)), min.gapwidth = Inf)

try Machine$integer.max instead of Inf; but maybe you're looking for 
range(GRanges(<...>))

Martin

>
> rather than
>
> reduce(GRanges("chr1", IRanges(c(1000, 9000), width = 1000)), min.gapwidth = 999999999)
>
> --------------------------------------
> Dario Strbenac
> Research Assistant
> Cancer Epigenetics
> Garvan Institute of Medical Research
> Darlinghurst NSW 2010
> Australia
>
> _______________________________________________
> Bioc-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/bioc-devel


-- 
Computational Biology
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N. PO Box 19024 Seattle, WA 98109

Location: M1-B861
Telephone: 206 667-2793



More information about the Bioc-devel mailing list