data(simChIP)
## scatter plots
library(geneplotter)
simChIP[ , 3:14] <- logb(simChIP[ , 3:14], 2)
par(mfrow=c(2,2))
smoothScatter(simChIP[ , 7], simChIP[ , 11], xlab = "H3 sample 1",
ylab = "H3K27me3 sample 1")
smoothScatter(simChIP[ , 8], simChIP[ , 12], xlab = "H3 sample 2",
ylab = "H3K27me3 sample 2")
smoothScatter(simChIP[ , 9], simChIP[ , 13], xlab = "H3 sample 3",
ylab = "H3K27me3 sample 3")
smoothScatter(simChIP[ , 10], simChIP[ , 14], xlab = "H3 sample 4",
ylab = "H3K27me3 sample 4")Run the code above in your browser using DataLab