problem tick marker and text
Duijvesteijn, Naomi wrote:
Hi R-ers, I am struggling with my x-axis in a association plot. What I would like is to place the labels of the x-axis between the tick markers and normally the labels are printed at the place where the tick marker is placed. I don???t want to move the tick marker (it gives the switch between one chromosome and the next) but I just want to put the chromosome number in between the two tick markers. This example from mhtplot has the chromosome number at the place of the tickmarker. How can I change this to ???in between???? Thanks!, Naomi
Kudos for including reproducible code, BUT
1. do you really consider a dataframe of
dimension ~500K by 3 to be "minimal"?
2. there are over 2000 R-packages on CRAN;
wouldn't it be polite to say which one contains mhtplot()?
# fake example with Affy500k data
affy <-c(40220, 41400, 33801, 32334, 32056, 31470, 25835, 27457, 22864, 28501,
26273,
24954, 19188, 15721, 14356, 15309, 11281, 14881, 6399, 12400, 7125, 62
07)
CM <- cumsum(affy)
is CM used anywhere?
n.markers <- sum(affy)
n.chr <- length(affy)
test <- data.frame(chr=rep(1:n.chr,affy),pos=1:n.markers,p=runif(n.markers))
oldpar <- par()
par(las="1",cex=0.6)
colors <- rep(c("blue","green"),11)
par(cex.axis=1.3)
mhtplot(test,usepos=TRUE,colors=colors,gap=10000,pch=19,bg=colors)
box()
par(oldpar)
Here is one way, using about half the data. I use mtext() to place the labels, since axis= seems to be hard-coded in mhtplot. # make a vector of labels to place between ticks textvec <- letters[1:5] # plot mhtplot(test[1:2e5,],usepos=TRUE,colors=colors,gap=10000,pch=19,bg=colors) # now figure out where to put the labels. This would be # easy with seq() if the scale were uniform but it's not. # So we use the locator() function. xat <- locator()$x xat <- round(xat) xat #[1] 34316 85096 132458 173961 216929 # place the labels; adjust line= if you wish mtext(textvec, side=1, line=0, at=xat) If I had to do this frequently, I would hack the mhtplot() function. -Peter Ehlers
Disclaimer: De informatie opgenomen in dit bericht (en bijlagen) kan vertrouwelijk zijn en is uitsluitend bestemd voor de geadresseerde(n). Indien u dit bericht ten onrechte ontvangt, wordt u geacht de inhoud niet te gebruiken, de afzender direct te informeren en het bericht te vernietigen. Aan dit bericht kunnen geen rechten of plichten worden ontleend. ---------------------------------------------------------------------------- ---------------------------- Disclaimer: The information contained in this message may be confidential and is intended to be exclusively for the addressee. Should you receive this message unintentionally, you are expected not to use the contents herein, to notify the sender immediately and to destroy the message. No rights can be derived from this message. P Please consider the environment before printing this email ------------------------------------------------------------------------
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.