# \donttest{
# estimate a binomial logit model using example data
library(UPG)
data(titanic)
y = titanic[,1]
Ni = titanic[,2]
X = titanic[,-c(1,2)]
results.binomial = UPG(y = y, X = X, Ni = Ni, model = "binomial")
# plot the results and sort coefficients by effect size
plot(results.binomial, sort = TRUE)
# plot only variables 1 and 3 with custom names, credible intervals and axis labels
plot(results.binomial,
include = c(1,3),
names = c("Custom 1", "Custom 2"),
q = c(0.1, 0.9),
xlab = c("Custom X"),
ylab = c("Custom Y"))
# }
Run the code above in your browser using DataLab