In altaf-ali/regsim: Regression Model Simulations

library(foreign)
library(regsim)

bes <- read.dta("http://uclspp.github.io/PUBLG100/data/bes.dta")

bes$Gender <- factor(bes$Gender, levels = c(0, 1), labels = c("Female", "Male"))

model <- glm(Turnout ~ Income + polinfoindex + Influence + Gender + Age + 
               edu15 + edu17 + edu18 + edu19plus + in_school + in_uni,
             family = binomial,
             data = bes)

summary(model)

set.seed(666)

x <- list(Gender = levels(bes$Gender))

sim <- regsim(model, x)

summary(sim)

x <- list(Gender = levels(bes$Gender), Income = 1:13)

sim <- regsim(model, x)

plot(sim, ~Income + Gender, 
     ylab = "Predicted Probability of Voting",
     main = "Effect of Income by Gender")

altaf-ali/regsim documentation built on May 9, 2019, 9:57 p.m.