swapped out dplyr::ntile for Hmisc::hdquantile