Instantly share code, notes, and snippets.
library(data.table) q <- c(0.001, 0.002, 0.003, 0.003, 0.004, 0.004, 0.005,