install.packages("flexdashboard")
install.packages("tidyverse")
source("Functions/ReadData.R")
#backcolor="azure2" #background color for plots
#gridcolor="azure1"
fillcolor = "skyblue4"
source("Functions/WelfareFunctions.R")
printFigures = T
today = Sys.Date()
install future.apply
install magrittr
install.packages(c("magrittr", "xtable", "future.apply", "forcats"))
source("Functions/TreatmentAssignmentFunctions.R")
#backcolor="azure2" #background color for plots
#gridcolor="azure1"
fillcolor = "skyblue4"
source('~/Dropbox (Personal)/policyexperiments/R/replication_files/replication_files_kasy_sautmann/simulations/CalibratedSimulations.R')
source("Functions/WelfareFunctions.R")
setwd("~/Dropbox (Personal)/policyexperiments/R/replication_files/replication_files_kasy_sautmann/simulations")
#backcolor="azure2" #background color for plots
#gridcolor="azure1"
fillcolor = "skyblue4"
DataList = ReadData(printFigures)
DataList = ReadAllData(printFigures)
DoHistograms = T
printFigures = T
source("Functions/ReadData.R")
source("Functions/ReadData.R")
source("Functions/calibratedSimulationFunctions.R")
source("Functions/WelfareFunctions.R")
source("Functions/TreatmentAssignmentFunctions.R")
#backcolor="azure2" #background color for plots
#gridcolor="azure1"
fillcolor = "skyblue4"
today = Sys.Date()
printFigures = T
DoHistograms = T
DataList = ReadAllData(printFigures)
#number of replications
MC_replicates = 100000
# parallelize simulations
plan(multiprocess)
#which methods of treatment assignment to simulate, and in which order to display them
methods = c(4,3,1)
#multiples of original sample size
multiplicities = c(1) #c(.5, 1, 1.5)
#number of waves
wavenumbers = c(2, 4, 10)
#columnames for tables
columnames = paste(wavenumbers, "waves")
for (Mult in multiplicities) {
for (application in 1:length(DataList)) {
applicationData = rep(DataList[application], length(wavenumbers))
for (i in 1:length(wavenumbers)) {
applicationData[[i]]$wavesizes =
rep(floor(applicationData[[i]]$N * Mult / wavenumbers[i]), wavenumbers[i])
applicationData[[i]]$waves = wavenumbers[i]
}
filename = paste(today, "_CalibratedSimulations_",
DataList[[application]]$filename,
"_" , Mult, sep = "")
DesignTable(applicationData, methods, MC_replicates,
columnames, filename, DoHistograms)
}
}
