inst/replication_code/table_code/tab9_3.r

In psre: Presenting Statistical Results Effectively

## The psre package must be installed first.  
## You can do this with the following code
# install.packages("remotes")
# remotes::install_github('davidaarmstrong/psre')
## load packages
library(tidyverse)
library(psre)
library(stargazer)
library(car)

## load data from psre package
data(wvs)

## recode civilization into three categories
## make the democratic history factor
wvs1 <- wvs %>% 
  mutate(
    civ2 = case_when(
      civ == 6 ~ "Latin American", 
      civ == 9 ~ "Western", 
      TRUE ~ "Other"), 
    civ2 = factor(civ2, levels=c("Western", "Latin American", "Other")), 
    democrat= factor(democrat, levels=1:2, 
                     labels=c("New Democracy", "Established Democracy")))

## estimate interactive model
intmod2 <- lm(resemaval ~ civ2 * pct_secondary, data=wvs1)

## print model 
stargazer(intmod2, digits=2, type="text")

## incremental F-test for interaction
Anova(intmod2, type="II", test="F")