library("topmodels")
library("topmodels")
data("SerumPotassium", package = "topmodels")
## Figure 9.3a-c from Rice (2007), and actual hanging rootogram
## (note that Rice erroneously refers to suspended rootograms as hanging)
sp <- lm(SerumPotassium ~ 1)
br <- 32:54/10 - 0.05
rootogram(sp, scale = "raw", style = "standing",
breaks = br, col = "transparent")


rootogram(sp, breaks = br)