spplinescol<- function(process,color) { par(col=color) tim = process$time x = process$x n = length(tim) tmax = tim[n] * (1 + 1/n) tima = tim[2:n] time = c(tim[1], as.vector(t(cbind(tima, tima))), tmax) sp = as.vector(t(cbind(x, x))) lines(time, sp, type = "l") }