4.4 Estimación del coeficiente de variación

Este bloque de código tiene como objetivo calcular la desviación estándar (sd) y el coeficiente de variación (cv) de los parámetros estimados del modelo para las áreas de observación (theta_obs) y las áreas de predicción (theta_pred). En primer lugar, se convierten los valores de theta_obs y theta_pred en matrices y se les asignan nombres de columnas correspondientes a los tres estados laborales: Ocupado_mod, Desocupado_mod e Inactivo_mod. A continuación, se calcula la desviación estándar de cada uno de estos parámetros para las áreas de observación y se crea una nueva matriz theta_obs_ordenado_sd. De manera similar, se calcula la desviación estándar de cada parámetro para las áreas de predicción y se crea una matriz theta_pred_ordenado_sd. Posteriormente, se calcula el coeficiente de variación para cada parámetro en las áreas de observación y predicción, dividiendo la desviación estándar entre la media de cada parámetro y se crean las matrices theta_obs_ordenado_cv y theta_pred_ordenado_cv, respectivamente. Finalmente, se combinan todas las matrices en una sola llamada estimaciones_obs mediante la función cbind().

theta_obs_ordenado_sd <- matrix(theta_obs[,"sd"], 
                             nrow = D,
                             ncol = P,byrow = TRUE) 

colnames(theta_obs_ordenado_sd) <- c("Ocupado_mod_sd", "Desocupado_mod_sd", "Inactivo_mod_sd")
theta_obs_ordenado_sd%<>% as.data.frame()
theta_obs_ordenado_sd <- cbind(dam2 = indicador_dam1$dam2,
                            theta_obs_ordenado_sd)
theta_obs_ordenado_cv <- theta_obs_ordenado_sd[,-1]/theta_obs_ordenado[,-1]

colnames(theta_obs_ordenado_cv) <- c("Ocupado_mod_cv", "Desocupado_mod_cv", "Inactivo_mod_cv")

theta_obs_ordenado_cv <- cbind(dam2 = indicador_dam1$dam2,
                               theta_obs_ordenado_cv)

theta_pred_ordenado_sd <- matrix(theta_pred[,"sd"], 
                              nrow = D1,
                              ncol = P,byrow = TRUE)

colnames(theta_pred_ordenado_sd) <- c("Ocupado_mod_sd", "Desocupado_mod_sd", "Inactivo_mod_sd")
theta_pred_ordenado_sd%<>% as.data.frame()
theta_pred_ordenado_sd <- cbind(dam2 = dam_pred$dam2, theta_pred_ordenado_sd)

theta_pred_ordenado_cv <- theta_pred_ordenado_sd[,-1]/theta_pred_ordenado[,-1]

colnames(theta_pred_ordenado_cv) <- c("Ocupado_mod_cv", "Desocupado_mod_cv", "Inactivo_mod_cv")

theta_pred_ordenado_cv <- cbind(dam2 = dam_pred$dam2, theta_pred_ordenado_cv)