r - Label on X-axis for each group of points -

-

i'm creating script allows me data database, , visualises in graph. enter image description here

as can see, there 8 groups of data, indicated on x-axis. each group contains 90 values.

right now, place of labels hard-coded this: 

axis(1, @ = c(31.25, 93.75, 156.25, 218.75, 281.25, 343.75, 406.25, 468.75),      labels = c("ss oligo 1", "ss oligo 2", "ss oligo 3", "ss oligo 4",                 "ss oligo 4", "ss oligo 5", "ss oligo 6", "ss oligo 7"))

it works fine, wondering if there way more dynamically, telling r assign label each set of 90 values.

example: 

# generate data ################################################################################### x <- vector() y <- vector()  # y[length(y)+1] <- sample(10:12, 1, replace = true) oligo_1 <- runif(62, 10.5, 11.5) oligo_2 <- runif(62, 14, 15) oligo_3 <- runif(62, 17, 18) oligo_4a <- runif(64, 20.5, 22) oligo_4b <- runif(64, 20.5, 22) oligo_5 <- runif(62, 24, 25) oligo_6 <- runif(62, 27, 28) oligo_7 <- runif(62, 30, 31)  y <- c(oligo_1, oligo_2, oligo_3, oligo_4a, oligo_4b, oligo_5, oligo_6, oligo_7) x <- c(1:500)  example <- data.frame(x, y)  # define variables ################################################################################ xmin      <- 10 xmax      <- 36  # generate graph ################################################################################### png(filename = "graph.png", width = 1500, height = 833)  plot(x = example[,2], type="l",      xlim = c(0, nrow(example)), ylim = c(xmin, xmax),      xaxt="n", yaxt="n",      xlab = "", ylab = "")  rect(xleft = par("usr")[1], ybottom =  9.8, xright = par("usr")[2], ytop = 12.2, border = "lightgrey", col = "lightgrey") rect(xleft = par("usr")[1], ybottom = 13.2, xright = par("usr")[2], ytop = 15.5, border = "lightgrey", col = "lightgrey") rect(xleft = par("usr")[1], ybottom = 16.5, xright = par("usr")[2], ytop = 18.9, border = "lightgrey", col = "lightgrey") rect(xleft = par("usr")[1], ybottom = 19.9, xright = par("usr")[2], ytop = 22.3, border = "lightgrey", col = "lightgrey") rect(xleft = par("usr")[1], ybottom = 23.3, xright = par("usr")[2], ytop = 25.5, border = "lightgrey", col = "lightgrey") rect(xleft = par("usr")[1], ybottom = 26.5, xright = par("usr")[2], ytop = 28.7, border = "lightgrey", col = "lightgrey") rect(xleft = par("usr")[1], ybottom = 29.7, xright = par("usr")[2], ytop = 32.1, border = "lightgrey", col = "lightgrey")  axis(1, @ = c(31.25, 93.75, 156.25, 218.75, 281.25, 343.75, 406.25, 468.75),       labels = c("ss oligo 1", "ss oligo 2", "ss oligo 3", "ss oligo 4",                  "ss oligo 4", "ss oligo 5", "ss oligo 6", "ss oligo 7")) axis(2, @ = c(11, 14.35, 17.7, 21.1, 24.4, 27.6, 30.9), las = 1)  lines(x = example[,2]) box()  mtext(paste("qc-check for", "test"), side = "3", line = 1, cex = 2, font = 2) mtext("samples"                       , side = "1", line = 3, cex = 1, font = 1)  legend(x = par("usr")[1]+10, y = par("usr")[4]-1, legend = c("cq", "ccq"), cex=1.5, lwd = 2, col = c("black","red"))  dev.off()

why not:

axis(1, at=31.25+(0:7)*62.5, labels=paste("ss oligo",1:8) )

Comments

Post a Comment

Popular posts from this blog

c# - How to get the current UAC mode -

-
can check current uac type set in device before running application through c#. i got know approach check whether current user administrator or not. but, need know uac type set (i.e. default/never notify) internal logic. is possible? please help. this article self-elevation , checking level program runs in. due security reasons i'm afraid won't able retrieve uac level set. http://code.msdn.microsoft.com/windowsdesktop/csuacselfelevation-644673d3
Read more

postgresql - Lazarus + Postgres: incomplete startup packet -

-
i've been building lazarus pascal program using postgresql on back-end. used work fine following lines opening db. dbconn:= tpqconnection.create(nil); dbconn.hostname := 'localhost'; dbconn.databasename:= 'dbhrs'; dbconn.username:='mizk'; dbconn.password:='123'; dbconn.open; if dbconn.connected openhrsdb := true else openhrsdb := false; but moment change localhost server's ip (on lan), program stops. no errors or warnings. have no clue happening. here related details may narrow down problem. the program i'm running workstation on same lan server. both machines run ubuntu 12.04 desktop i can access postgres db on server using pgadmin iii, guess it's not firewall issue. i have allowed postgresql ufw , ufw disabled. strangely, when firewall disabled, can ssh server terminal using workstation. when it's enabled, can't my biggest concern right why pascal program not working when i...
Read more

javascript - Ajax jqXHR.status==0 fix error -

-
$.ajax({ url: urlstring, datatype: "json", type: "get", success: function (data) { alert(data); }, error: function (jqxhr, exception) { if (jqxhr.status === 0) { alert('not connect.\n verify network.'); } else if (jqxhr.status == 404) { alert('requested page not found. [404]'); } else if (jqxhr.status == 500) { alert('internal server error [500].'); } else if (exception === 'parsererror') { alert('requested json parse failed.'); } else if (exception === 'timeout') { alert('time out error.'); } else if (exception === 'abort') { alert('ajax request aborted.'); } else { alert('uncaught error.\n' + jqxhr.responsete...
Read more