data(war1800)
f1 <- esc + war ~ s_wt_re1 + revis1 | 0 | regime1 | balanc + regime2
m1 <- egame12(f1, data = war1800, boot = 10)
pp1 <- predProbs(m1, x = "balanc", n = 5) ## continuous x
pp2 <- predProbs(m1, x = "regime1") ## discrete x
## if "ask" is FALSE and "which" isn't specified, all plots are printed
op <- par(mfrow = c(2, 2))
plot(pp1)
par(op)
plot(pp1, ask = TRUE)
# displays menu:
# Make a plot selection (or 0 to exit):
# 1: plot: Pr(~esc)
# 2: plot: Pr(esc,~war)
# 3: plot: Pr(esc,war)
# 4: plot all terms
## to change line type for confidence bounds, use argument "lty.ci"
plot(pp1, which = 3, lty.ci = 3)
## all the standard plotting options work too
plot(pp1, which = 3, xlab = "Capabilities", ylab = "Probability", main = "Title")
## discrete x variables are plotted via R's boxplot functionality
plot(pp2, which = 3)
Run the code above in your browser using DataLab