# NOT RUN {
library(dplyr) # contains pipe operator %>%
data(turnout)
# plot first differences
zelig(Fertility ~ Education, data = swiss, model = 'ls') %>%
      setx(z4, Education = 10) %>%
      setx1(z4, Education = 30) %>%
      sim() %>%
      plot()
# }
Run the code above in your browser using DataLab