library(Zelig)

data(coalition)

z <- zelig(duration ~ fract + numst2, model = "gamma", data = coalition)

x.low <- setx(z, numst2 = 0)
x.high <- setx(z, numst2 = 1)

s <- sim(z, x = x.low, x1 = x.high, num = 10)
