data(asdgap)
names(asdgap)
X <- asdgap$X
numcol <- which(colnames(X) == "1000" | colnames(X) == "1800")
numcol
plotsp(X, lwd = 1.5)
abline(v = as.numeric(colnames(X)[1]) + numcol - 1, col = "grey", lty = 3)
Run the code above in your browser using DataLab