data(simChIP)
## Not run:
# ## 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")
# ## End(Not run)
Run the code above in your browser using DataLab