[BioC] plotting probes along a chromosome

J.Oosting at lumc.nl J.Oosting at lumc.nl
Sat Jul 16 19:26:41 CEST 2011


A similar question with a number of possible solutions can be found on BioStar
http://biostar.stackexchange.com/questions/378/drawing-chromosome-ideogams-with-data/

Jan

________________________________________
Hello List,

I am trying to graph the position of probes ~50 nucleotides in length along
a chromosome and am looking for the best way to accomplish this. So far I
have tried the genomeGraphs package and can plot a chromosome and its bands
just fine. Next I attempt to plot assign the data coordinates using
>expres = new("GenericArray", intensity = intensity, probeStart =
exonProbePos, dp = DisplayPars(color="darkred", type="point"))
function. This accepts the values just fine. However, when I attempt to plot
it using ;

>
gdPlot(list(ideog,expres), minBase = minbase, maxBase =maxbase), I get the
following error:

Error in unit(x, default.units) : 'x' and 'units' must have length > 0

I am confused by this error because the input files that I give to "expres"

seem to match the ones given by the data provided by the package

I import the data:

"chr" "location" "strand" "expr"
"1" 4 180491934 -1 1.0
"2" 4 180491939 1 1.0
"3" 4 180491945 -1 1.0
"4" 4 180491963 1 1.0
"5" 4 180491993 -1 1.0
"6" 4 180491996 1 1.0
"7" 4 180491999 -1 1.0

using
>z<-read.table('p', header=TRUE)
the expres variable seems to accept a matrix of type "double"
which I build using
>expres_matrix<-as.matrix(z$expr)
and a list of "integers" wihich I give as z$location

when I do a class and typeof on these variables, they match the class/type
of the ones provided in the package but
when I run the plotting program, I get the error

I was initially trying this as a quick and dirty technique but am now stuck,
If anyone can suggest alternate/better methods, I am all ears.

Thanks

        [[alternative HTML version deleted]]

_______________________________________________
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