data(coalition)
# Fit the statistical model
z.out <- zelig(duration ~ fract + numst2, model = "gamma", data = coalition)
## Setting the explanatory variables at their default values
## (mode for factor variables and mean for non-factor variables),
## with numst2 set to the vector 0 = no crisis, 1 = crisis.
x.low <- setx(z.out, numst2 = 0)
x.high <- setx(z.out, numst2 = 1)
## Simulating draws using the default bootstrap method.
s.out <- sim(z.out, x = x.low, x1 = x.high)
# Summary of fitted model
summary(z.out)
# Summary of simulated quantities of interest
summary(s.out)
# Plot of simulated quantities of interest
plot(s.out)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.