[BioC] correlation image and tickmarks
Arne.Muller at aventis.com
Arne.Muller at aventis.com
Tue Aug 3 10:21:32 CEST 2004
Hello,
I've a problem aligning tickmarks to an image. I've created a correlation matrix for 84 affy mouse chips. I'm visualizing the matrix as an image with colour coding according to the correlation coefficient.
The 84 chips are distributed over three factors, but the desgin is unbalanced, so that the tickarks and the lables for the axis are must not evenly distributed. A regular grid via the 'grid' function aligns perfectly with the image cells, but the tickmarks seem to be shifted, and not aligned perfectly with the image cells. This offset is even stonger for the y-axis. The thing is that I don't want 84 lables at the axis, it's enough to have one lable for all the different factor level combinations, which results in 28 labels.
Maybe you have an idea how to setup the command to align the tick marks.
kind regards,
Arne
Here are my commands:
> library(marrayPlots) # for the colors
> col <- maPalette(low='white', high='darkred', k=50)
> par(ps=8, cex=1, mar=c(1,5,5,1)) # space needed for lables @ axis 1 and 3
# a x and y tange from 1 to 84, x is the correlation matrix (dim = 84x84)
> image(1:84, 1:84, x, col=col, xaxt='n', yaxt='n', xlab='', ylab='')
# set up the axis
> axis(3, i, labels=names(l), las=2, tick=T)
> axis(2, i, labels=names(l), las=2, tick=T)
> grid(84, col='black', lty='solid')
# this is where i get the lables from
> l
NEW:4:0 NEW:4:100 NEW:4:250 NEW:4:500 NEW:4:1000 NEW:24:0
3 3 3 3 3 3
NEW:24:100 NEW:24:250 NEW:24:500 NEW:24:1000 OLD:4:0 OLD:4:100
3 3 3 3 4 3
OLD:4:250 OLD:4:500 OLD:4:1000 OLD:24:0 OLD:24:100 OLD:24:250
2 3 3 4 3 2
OLD:24:500 OLD:24:1000 PRG:4:0 PRG:4:100 PRG:4:250 PRG:4:1000
3 3 3 3 3 3
PRG:24:0 PRG:24:100 PRG:24:250 PRG:24:1000
3 3 3 3
# these are the positions along the axis for the tick marks, technical
# replicates from 1 to 3 (replicates of one factor level combination), 4 to 6
# ...
> i
[1] 3 6 9 12 15 18 21 24 27 30 34 37 39 42 45 49 52 54 57 60 63 66 69 72 75
[26] 78 81 84
More information about the Bioconductor
mailing list