Nothing
source("setup.R")
data("db", package = "gamlss.data")
db$lage <- with(db, age^(1/3))
mod <- ctm(bbs(head, df = 2.5) ~ bbs(lage, df = 3, knots = 50), data = db,
family = Binomial(link = "probit"), monotone = FALSE, ngrid = 50,
control = boost_control(nu = .5, trace = TRUE))
l <- with(db, seq(from = min(lage), to = max(lage), length = 100))
pr <- predict(mod, newdata = data.frame(lage = l), type = "response",
annotate = TRUE)
pr$lage <- l^3
pr$cut <- factor(pr$lage > 2.5)
levels(pr$cut) <- c("Age < 2.5 yrs", "Age > 2.5 yrs")
save(db, pr, file = "ex_heads.Rda")
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.