Hello Alex, Here is one way to do it. It works but it's not pretty :)
interp <- approx(c(2, 62), c(3, 34), method="linear", xout=2:62) m <- matrix(c(interp$x, round(interp$y)), ncol=2) tie <- m[,2] == c(-Inf, m[-nrow(m),2]) m <- m[ !tie, ] You might want to examine the result like this... plot(m) # plots points lines(c(2,26), c(3, 34)) # overlay line for comparison Michael ______________________________________________ R-help@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.
