# NOT RUN {
data(Roaddamage)
y <- Roaddamage$damage
x <- Roaddamage$population
discmethod <- c("equal","natural","quantile","geometric","sd")
discitv <- c(3:7)
optidisc(y, x, discmethod, discitv)
# }
Run the code above in your browser using DataLab