1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
data(coalition2)
user.prompt()
z.out <- zelig(duration ~ invest + fract + polar + numst2 + crisis +
tag(1 | country),
data=coalition2, model="gamma.mixed", method="PQL",
family=Gamma(link="log"))
user.prompt()
summary(z.out)
## Setting the explanatory variables at their default values
## (mode for factor variables and mean for non-factor variables),
## comparing the ruling coalition in the minority to the ruling
## coalition in the majority.
user.prompt()
x.high <- setx(z.out, numst2 = 1)
x.low <- setx(z.out, numst2 = 0)
## Simulating draws using the default bootstrap method.
user.prompt()
s.out <- sim(z.out, x = x.low, x1=x.high)
user.prompt()
## Viewing the simulated quantities of interest, for every
## observation:
summary(s.out)
user.prompt()
|
