if (require("TH.data") && require("tram")) {
data("bodyfat", package = "TH.data")
### estimate unconditional model
m_mlt <- BoxCox(DEXfat ~ 1, data = bodyfat, prob = c(.1, .99))
### get corresponding in-sample log-likelihood
logLik(m_mlt)
### estimate conditional transformation model
bm <- ctmboost(m_mlt, formula = DEXfat ~ ., data = bodyfat,
method = quote(mboost::mboost))
### in-sample log-likelihood (NEEDS TUNING OF mstop!)
logLik(bm)
### evaluate conditional densities for two observations
predict(bm, newdata = bodyfat[1:2,], type = "density")
}
Run the code above in your browser using DataLab