How can I plot square wave data in R?

Question

Consider the following data, where the left column represents the bit (1 or 0), and the right column represents the number of microseconds that we observe the bit.

0 664 1 63 0 404 1 544 0 651 1 686 0 507 1 1155 0 664 1 271 0 456 1 2763 0 664 1 115 0 456 1 4010 0 664 1 63 0 351 1 3855 

I would like to build this data so that there is a horizontal line at 0 with a width of 664, followed by a rise to a horizontal line by 1 with a width of 63, followed by a fall on a horizontal line at 0 with a width of 404, etc.

Is there an efficient and direct way to build this in R that does not require manual comparison with constraints?

Here is my current code for this, which is extremely inefficient and naive, so I hope there is a better way.

 args <- commandArgs(trailingOnly = TRUE) data = read.table(args[1]) current = 1 sumA = 0 pf = function(x) { if (x < sumA) { return(data[current,1]) } for (i in current: length(data[,1])) { sumA <<- sumA + data[i,2] if (x < sumA) { current <<- i + 1 return(data[i,1]) } } return("OUT OF BOUNDS") } cumSum = colSums(data)[[2]] print(cumSum - 1); h = Vectorize(pf) plot(h, 1, cumSum-1, n=cumSum-1, lwd=0.001, xlim=c(0,cumSum-1)) 

Answer 1

As mentioned in my comment, the plot command with the type flag set to s should do the trick.

For example, for the first 10 samples:

 x <- c(0,664,63,404,544,651,686,507,1155,664,271) xC <- cumsum(x) y <- c(0,1,0,1,0,1,0,1,0,1,0) plot(xC,y,type='s')