Zelig4: Everyone's Statistical Software

library(Zelig)

#####  Example 1: Basic Example with First Differences  #####

# Attach sample data and variable names:  
data(macro)

# Estimate model and present a summary:

z.out1 <- zelig(unem ~ gdp + capmob + trade, model = "normal", data = macro)


# Set explanatory variables to their default (mean/mode) values, with
# high (80th percentile) and low (20th percentile) values:

x.high <- setx(z.out1, trade = quantile(macro$trade, 0.8))
x.low <- setx(z.out1, trade = quantile(macro$trade, 0.2))

# Generate first differences for the effect of high versus low trade on
# GDP:

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

# Summarize the fitted model

summary(z.out1)

# Summarize the simulated quantities of interest

summary(s.out1)

# Plot the simulated quantities of interest

plot(s.out1)