## Not run: ------------------------------------
# #?MASS::bandwidth.nrd
# bandwidth_rot <- function (x) {
# r <- quantile(x, c(0.25, 0.75))
# h <- (r[2L] - r[1L])/1.34
# 4 * 1.06 * min(sqrt(var(x)), h) * length(x)^(-1/5)
# }
## ---------------------------------------------
Run the code above in your browser using DataLab