Hi everyone,
I'm trying to add another 'y axis' to a heatmap that I generate
thusly:
>
heatmap(t(pemMatrix),Rowv=NA,Colv=NA,col=heat.colors(3),scale="row",main="PEM
scores by Tissue Type", xlab = "Tissues", ylab =
"Contigs",margins = c(7,8),
na.rm = TRUE)
The first y axis (on side 4, in axis() terms) indexes the columns in the
transpose of pemMatrix (numeric data values). I want to add another y axis
on the other side (side 2 in axis terms) to display the results of another
test. I've checked the wiki and ?image, ?heatmap, and have been able to add
something like what I want through the following:
>
heatmap(t(pemMatrix),Rowv=NA,Colv=NA,col=heat.colors(3),scale="row",main="PEM
scores by Tissue Type", xlab = "Tissues", ylab =
"Contigs",margins = c(7,8),
na.rm = TRUE)> par(new = TRUE, las=1)
> axis(2,at = axTicks(4), labels = as.character(kwResults[,2]))
but the axis on side 2 is never scaled as the first axis on side 4. I've
also checked par(), and the ylim for the graph is the same as what I set
with at in the call to axis(). Can anyone shed some light on how heatmap is
generating the y-axis on side 4?
Thanks,
Lee Zamparo
[[alternative HTML version deleted]]