summary(los) # quite skewed, with median(.) = 8
plot(table(los))
boxplot(los, horizontal=TRUE, add=TRUE, col = "red", axes=FALSE)
##-> "outliers" instead of "just skewed"
hist(log(los))
boxplot(log(los), add=TRUE, col=2, border=2, horizontal = TRUE, at = -1)
if(R.version$arch == "x86_64") { ## FIXME -- bug in C code
## Hubert and Vandervieren (2006), p. 15, Fig. 11.
adjbox(los, col = "gray", outcol = "red")
boxplot(los, add = TRUE, staplewex = 0.2, outcex = 0.5, outpch = 4)
title("adjbox(los) and boxplot(los) superimposed")
}
Run the code above in your browser using DataLab