data <- rnorm(100)
shorth(data)
shorth(data, sample(1:100, size = 20))
# Calculating a CI for shorth using [mboot()]
data <- rnorm(100)
boot.out <- mboot(data, shorth, m = sqrt(length(data)))
basic.ci <- mboot.ci(boot.out, conf =0.95, tau = function(n) return(n^(1/3)), types = "basic")$basic
Run the code above in your browser using DataLab