library(Zelig)
data(coalition)
z.out <- zelig(duration ~ fract + numst2 + crisis, model = "gamma", data = coalition[1:100, ])
x.low <- setx(z.out, fract=300, numst2 = 0, crisis=200)
x.high <- setx(z.out, fract=300, numst2 = 1, crisis=200)
s.out <- sim(z.out, x = x.low, x1 = x.high, num = 10, bootstrap=TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.