####################################################################################################
#
# DYNAMIS-POP Parameter Generation File 15 Primary Education Grade Tracking
# This file is generic and works for all country contexts.
# Input file: globals_for_analysis.RData (To generate such a file run the setup script)
# Last Update: Martin Spielauer 2018-06-29
#
####################################################################################################
####################################################################################################
# Clear work space, load required packages and the input object file
####################################################################################################
rm(list=ls())
library(haven)
library(dplyr)
library(data.table)
library(sp)
library(maptools)
library(survival)
library(fmsb)
library(eha)
load(file="globals_for_analysis.RData")
parafile <- file(g_para_primaryeductrack, "w") # Output Parameter File
cat("parameters {\n", file=parafile)
cat(" double StartSchoolOneYear = ",g_school_season_start,"; //EN Start at school season\n", file=parafile, append=TRUE)
cat(" double EndSchoolOneYear = ",g_school_season_end,"; //EN End of school season\n", file=parafile, append=TRUE)
cat(" //EN Grade repetition rate\n double SchoolOneRepetitionRate[EDUC_ONE_GEO][EDUC_ONE_GROUP][SIM_YEAR_RANGE] = { (101) 0.1, };\n", file=parafile, append=TRUE)
cat(" //EN School interruption rate\n double SchoolOneInterruptionRate[EDUC_ONE_GEO][EDUC_ONE_GROUP][SIM_YEAR_RANGE] = { (101) 0.1, };\n", file=parafile, append=TRUE)
cat(" //EN Distribution entry age\n cumrate EducOneEntryAge[EDUC_ONE_GEO][EDUC_ONE_GROUP][YOB_GRAD_PRIMARY][EDUC_ONE_ENTRY_AGE] = { (67) { 0.05, 0.85, (2) 0.05, }, };\n", file=parafile, append=TRUE)
cat(" //EN Distribution dropout grade \n cumrate EducOneDropoutGrade[EDUC_ONE_GEO][EDUC_ONE_GROUP][YOB_GRAD_PRIMARY][EDUC_ONE_GRADE] = { (67) { 4, 5, (4) 1, }, };\n", file=parafile, append=TRUE)
# end and close parameter output writing
cat("\n\n};", file=parafile, append=TRUE)
close(parafile)