# Script 22: Human Capital Index

>  This file creates the parameter for preschool attendance. This is a scenario parameter and filled with some default values only

## File output

The code below generates one model parameter stored in a Modgen .dat file

- Regional school quality parameters
- Coefficioents for HCI formula


## Code

```{r, message=FALSE, warning=FALSE}
####################################################################################################
# 
#  DYNAMIS-POP Parameter Generation File 22 - Human Capital Index
#  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 2019-11-25
#
####################################################################################################

####################################################################################################
# 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")
dat             <- g_residents_dat

# Constants
n_maxreg        <- max(dat$M_ROB)

####################################################################################################
# WRITE PARAMETERS
####################################################################################################

# Set Parameter Output File

parafile <- file(g_para_hci, "w")

cat("parameters { \n", file=parafile)
cat("//EN Schoo Quality\n", file=parafile, append=TRUE)

cat("double	SchoolQuality[REGION_NAT][SCHOOL_QUALITY_PARA] = {(", file=parafile, append=TRUE)
cat(n_maxreg,") {	0.75, 0.1 } \n	};\n", file=parafile, append=TRUE)
cat("double HCICoefficients[HCI_COEF] = {0.08, 0.65, 0.35};\n", file=parafile, append=TRUE)
cat("\n};\n", file=parafile, append=TRUE)
close(parafile)

```