Run »
values <- c(13,21,21,40,42,48,55,72) IQR(values)
[1] 28.75