[BioC] Geneplotter, alongChrom problem solved

Giulio Di Giovanni perimessaggini at hotmail.com
Wed Mar 23 10:42:44 CET 2005


Now it's ok. So, alongChrom require the multi-colors selection...

Thanks a lot !

Giulio

>From: Wolfgang Huber <huber at ebi.ac.uk>
>To: Giulio Di Giovanni <perimessaggini at hotmail.com>
>CC: bioconductor at stat.math.ethz.ch
>Subject: Re: [BioC] Geneplotter, alongChrom problem
>Date: Tue, 22 Mar 2005 22:47:36 +0100
>
>Hi Giulio,
>
>try
>
>alongChrom(eset, "Y", newChrom, xloc = "equispaced",
>   plotFormat = "image", scale = "zscale", colors = cm.colors(256))
>
>It seems that the default value of the 'colors' parameter (namely, just the 
>single value "red"), is not good for the image plot format.
>
>   Cheers,
>    Wolfgang
>
>Giulio Di Giovanni wrote:
>>Hi to all,
>>
>>I have a problem using alongChrom in Geneplotter pkg, with the "image" 
>>plot.
>>With this type of plot, all that I obtain is a single color filled image, 
>>instead the multicolor image. I didn't receive any error messages.
>>
>>I tried with different expression set and for different chromosomes, but 
>>with no differences. For example, the run from commonly known  "An 
>>introduction to some graphics in Bioconductor" pdf document is:
>>
>>
>>>library(annotate)
>>>library(geneplotter)
>>>library(hgu95av2)
>>>data(eset)
>>>newChrom <- buildChromLocation("hgu95av2")
>>>newChrom
>>
>>Instance of a chromLocation class with the following fields:
>>        Organism:  Homo sapiens
>>        Data source:  hgu95av2
>>        Number of chromosomes for this organism:  25
>>        Chromosomes of this organism and their lengths in base pairs:
>>                 1 : 246127941
>>                 2 : 243615958
>>                 3 : 199344050
>>                 4 : 191731959
>>                 5 : 181034922
>>                 6 : 170914576
>>                 7 : 158545518
>>                 8 : 146308819
>>                 9 : 136372045
>>                 10 : 135037215
>>                 11 : 134482954
>>                 12 : 132078379
>>                 13 : 113042980
>>                 14 : 105311216
>>                 15 : 100256656
>>                 16 : 90041932
>>                 17 : 81860266
>>                 18 : 76115139
>>                 19 : 63811651
>>                 20 : 63741868
>>                 21 : 46976097
>>                 22 : 49396972
>>                 X : 153692391
>>                 Y : 50286555
>>                 M : 16571
>>
>>>
>>>par(mfrow = c(2, 1))
>>>for (sc in c("max", "relative")) cPlot(newChrom, fg = "blue",
>>
>>+ scale = sc)
>>
>>>
>>>cols <- c("red", "green", "blue")
>>>cols <- cols[eset$cov3]    # needs eset - expression set
>>>par(mfrow = c(3, 2))
>>>alongChrom(eset, "1", newChrom, xloc = "equispaced", plotFormat = 
>>>"cumulative",
>>
>>+ col = cols, lwd = 2)
>><environment: 027EDB48>
>>
>>>alongChrom(eset, "1", newChrom, xloc = "physical", col = cols,
>>
>>+ lwd = 2)
>><environment: 021B72B4>
>>
>>>alongChrom(eset, "1", newChrom, xloc = "equispaced", plotFormat = 
>>>"local",
>>
>>+ col = cols, lwd = 2)
>>
>>>alongChrom(eset, "1", newChrom, xloc = "equispaced", plotFormat = 
>>>"local",
>>
>>+ col = cols, type = "p", pch = 16)
>>
>>>alongChrom(eset, "1", newChrom, xlim = c(87511280, 127717880),
>>
>>+ xloc = "equispaced", plotFormat = "local", col = cols, type = "p",
>>+ pch = 16)
>>
>>>alongChrom(eset, "Y", newChrom, xloc = "equispaced", plotFormat = 
>>>"image", scale="zscale")
>>>
>>
>>All the plots are OK, but the last one as I said is a completly red 
>>image...
>>I have the same result for all the datasets that I tried...
>>
>>Please, someone could suggest me what I'm doing wrong ?
>>
>>Thank you to all in advance,
>>
>>Giulio
>>
>>_______________________________________________
>>Bioconductor mailing list
>>Bioconductor at stat.math.ethz.ch
>>https://stat.ethz.ch/mailman/listinfo/bioconductor
>
>
>--
>Best regards
>   Wolfgang
>
>-------------------------------------
>Wolfgang Huber
>European Bioinformatics Institute
>European Molecular Biology Laboratory
>Cambridge CB10 1SD
>England
>Phone: +44 1223 494642
>Fax:   +44 1223 494486
>Http:  www.ebi.ac.uk/huber
>-------------------------------------



More information about the Bioconductor mailing list