[BioC] BUG? when setting mcols of a GRanges to a DataFrame one of whose columns is an ordered factor

Cook, Malcolm MEC at stowers.org
Tue Aug 13 22:29:13 CEST 2013


Great Herve,

Any chance of a back port to released version?

THx,

Malcolm

 >-----Original Message-----
 >From: Hervé Pagès [mailto:hpages at fhcrc.org]
 >Sent: Tuesday, August 13, 2013 2:33 PM
 >To: Cook, Malcolm
 >Cc: 'bioconductor at r-project.org'
 >Subject: Re: [BioC] BUG? when setting mcols of a GRanges to a DataFrame one of whose columns is an ordered factor
 >
 >Hi Malcolm,
 >
 >On 08/13/2013 11:39 AM, Cook, Malcolm wrote:
 >> Hiya,
 >>
 >> I find that setting mcol of a GRanges to a DataFrame one of whose columns is an ordered factor does "the wrong thing".
 >>
 >> Not sure how to approach this other than report it for now.
 >>
 >> Best,
 >>
 >> Malcolm
 >>
 >>
 >>
 >> viz:
 >>
 >>> library(GenomicRanges)
 >>> example(GRangesList)
 >>> gr1
 >> GRanges with 1 range and 2 metadata columns:
 >>        seqnames    ranges strand |     score        GC
 >>           <Rle> <IRanges>  <Rle> | <integer> <numeric>
 >>    [1]     chr2    [3, 6]      + |         5      0.45
 >>    ---
 >>    seqlengths:
 >>     chr2
 >>       NA
 >>> mc1<-mcols(gr1)
 >>> mc1$f<-factor('x')
 >>> mcols(gr1)<-mc1
 >>> gr1
 >> GRanges with 1 range and 3 metadata columns:
 >>        seqnames    ranges strand |     score        GC        f
 >>           <Rle> <IRanges>  <Rle> | <integer> <numeric> <factor>
 >>    [1]     chr2    [3, 6]      + |         5      0.45        x
 >>    ---
 >>    seqlengths:
 >>     chr2
 >>       NA
 >>>   mc1$fo<-ordered('x')
 >>> mcols(gr1)<-mc1
 >>> gr1
 >> GRanges with 1 range and 4 metadata columns:
 >>        seqnames      ranges  strand   |       score          GC          f          fo
 >>         "<Rle>" "<IRanges>" "<Rle>" "|" "<integer>" "<numeric>" "<factor>" Character,2
 >>    [1]   "chr2"    "[3, 6]"     "+" "|"         "5"      "0.45"        "x"         "x"
 >>    ---
 >>    seqlengths:
 >>     chr2
 >>       NA
 >
 >Thanks for the report. Note that this is a display bug only. Fixed in
 >IRanges 1.19.24:
 >
 >   > gr1
 >   GRanges with 1 range and 4 metadata columns:
 >         seqnames    ranges strand |     score        GC        f        fo
 >            <Rle> <IRanges>  <Rle> | <integer> <numeric> <factor> <ordered>
 >     [1]     chr2    [3, 6]      + |         5      0.45        x         x
 >     ---
 >     seqlengths:
 >      chr2
 >        NA
 >
 >Cheers,
 >H.
 >
 >
 >>> sessionInfo()
 >> R version 3.0.1 (2013-05-16)
 >> Platform: x86_64-unknown-linux-gnu (64-bit)
 >>
 >> locale:
 >>   [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C               LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8
 >LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8    LC_PAPER=C                 LC_NAME=C                  LC_ADDRESS=C
 >LC_TELEPHONE=C             LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
 >>
 >> attached base packages:
 >> [1] parallel  stats     graphics  grDevices datasets  utils     methods   base
 >>
 >> other attached packages:
 >> [1] GenomicRanges_1.12.4 IRanges_1.18.2       BiocGenerics_0.6.0
 >>
 >> loaded via a namespace (and not attached):
 >> [1] stats4_3.0.1 tools_3.0.1
 >>>
 >>
 >> _______________________________________________
 >> 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
 >>
 >
 >--
 >Hervé Pagès
 >
 >Program in Computational Biology
 >Division of Public Health Sciences
 >Fred Hutchinson Cancer Research Center
 >1100 Fairview Ave. N, M1-B514
 >P.O. Box 19024
 >Seattle, WA 98109-1024
 >
 >E-mail: hpages at fhcrc.org
 >Phone:  (206) 667-5791
 >Fax:    (206) 667-1319



More information about the Bioconductor mailing list