suppressMessages(library(cubfits, quietly = TRUE))
# fit Shah & Gilchrist (2011)
init.function(model = "roc")
fitlist <- fitMultinom(ex.train$reu13.df, ex.train$phi.Obs,
ex.train$y, ex.train$n)
ret.fit <- prop.model.roc(fitlist, phi.Obs.lim = range(ex.train$phi.Obs))
aa.names <- names(ex.train$reu13.df)
# plot.
par(mfrow = c(1, 3))
for(i.aa in 1:length(aa.names)){
plotmodel(ret.model = ret.fit[[i.aa]], main = aa.names[i.aa])
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.