From b09c9a95e2d7d8d1700130b1752fa80cfe3a8c5f Mon Sep 17 00:00:00 2001
From: Fabian Marder <fabian.marder@idiv.de>
Date: Fri, 24 Mar 2023 15:51:46 +0000
Subject: [PATCH] Upload New File
---
R_Scripts/main.R | 47 +++++++++++++++++++++++++++++++++++++++++++++++
1 file changed, 47 insertions(+)
create mode 100644 R_Scripts/main.R
diff --git a/R_Scripts/main.R b/R_Scripts/main.R
new file mode 100644
index 0000000..f52e5c4
--- /dev/null
+++ b/R_Scripts/main.R
@@ -0,0 +1,47 @@
+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))
+
--
GitLab