Upload New File

R_Scripts/main.R

+rm(list = ls())
+## Load  libraries
+library(apollo)
+library(ggplot2)
+library(ggrepel)
+library(ggpubr)
+library(dplyr)
+library(tidyr)
+library(tidylog)
+library(SciViews)
+library(texreg)
+library(lme4)
+library(haven)
+library(stargazer)
+library(sjPlot)
+library(sjmisc)
+library(viridis)
+library(hrbrthemes)
+source("R_Functions/helpfunctions.R")
+## Data preperation 
+source("R_Scripts/data_preperation.R")
+
+## Estimations
+source("R_Scripts/mnl_PMI.R") # Model 1
+source("R_Scripts/mnl_CMI.R") # Model 2
+source("R_Scripts/mnl_PMI_CMI.R") # Model 3
+source("R_Scripts/mnl_CMI_id_cont.R") # # Model 4
+source("R_Scripts/mnl_CMI_id_cont_sce.R") # Model 5
+
+## Outputs
+model_PMI_out <- quicktexregapollo(model = model_PMI)
+model_CMI_out <- quicktexregapollo(model = model_CMI)
+model_PMI_CMI_out <- quicktexregapollo(model = model_PMI_CMI)
+model_CMI_Q35_1_out <- quicktexregapollo(model = model_CMI_Q35_1)
+model_CMI_id_sce_out <- quicktexregapollo(model = model_CMI_id_sce)
+
+coefnames =c("ASC_A", "ASC_B", "ln(Profit)", "Hold", "Visible", "Community", "Utility", "Hold*PMI", "Hold*CMI","Hold*ID", "Hold*CMI*ID", "Hold*Mun", "Hold*CMI*Mun", "Hold*ID*Mun","Hold*CMI*ID*Mun", "Hold*Country", "Hold*CMI*Country", "Hold*ID*Country",  "Hold*CMI*ID*Country")
+## Table IV
+wordreg(list(m1=model_PMI_out, m2=model_CMI_out, m3=model_PMI_CMI_out, m4=model_CMI_Q35_1_out, m5=model_CMI_id_sce_out ),file = "Tables/mnl_models.docx", custom.coef.names = coefnames, digits = 3,  no.margin = TRUE ,stars = c(0.01, 0.05, 0.10))
+
+
+## Model 6 
+source("R_Scripts/lmer_Q64_1.R") 
+
+## Table VI
+wordreg(fit, file = "Tables/lmer_model.docx",  digits = 3,  no.margin = TRUE ,stars = c(0.01, 0.05, 0.10))
+