Validación de datos

Se crea en una lista que contiene el resumen del método DP2 para todos los años

for(i in tablas){
List_DP2 <- mget(paste0("ind_", tablas))
}

Orden de entrada de las variables

Se crea un data.frame de acuerdo al orden de importancia de las variables para todos los años Dependiendo del proceso, el DP2 adoptará diferentes valores. Por lo tanto, es importante que el método de como resultado un orden de entrada único de los indicadores parciales.

Variables_sort <- NULL
for(i in 1:3){
Variables_sort[i] <- as.data.frame(lapply(List_DP2, function(x) get(paste0("ind_", tablas[i]))$variables_sort))
}

Variables_sort <- do.call(cbind.data.frame, Variables_sort)
colnames(Variables_sort) <- c("2010", "2015", "2020")
Variables sort
2010 2015 2020
ANALF ANALF ANALF
SBASC OVPT SBASC
OVPT SBASC OVPT
PO2SM PO2SM OVSAE
OVSAE OVSAE VHAC
PL.5000 VHAC PL.5000
VHAC PL.5000 PO2SM
OVSEE OVSEE OVSDE
OVSDE OVSDE OVSEE

Coeficiente de correlación

El coeficiente de correlación es una medida que se utiliza para jerarquizar a los indicadores simples de acuerdo con el grado de correlación absoluta con respecto al indicador sintético resultante, es decir, ayuda a visualizar que variable tiene un mejor o peor apego con el fenómeno.

Se crea un data.frame de acuerdo a la correlación de cada variable con el indicador sintético (DP2) para todos los años

cor.coeff <- NULL
cor.coeff <- lapply(1:3, function(i) data.frame(rownames(get(paste0("ind_",tablas[i]))[["cor.coeff"]]),
                                                get(paste0("ind_",tablas[i]))$cor.coeff))
cor.coeff <- do.call(cbind.data.frame, cor.coeff)
colnames(cor.coeff) <- c("Ind_2010", "2010", "Ind_2015", "2015", "Ind_2020", "2020")
Correlation Coefficient
Ind 2010 Ind 2015 Ind 2020
ANALF 0.964 ANALF 0.961 ANALF 0.954
SBASC 0.891 SBASC 0.882 SBASC 0.887
OVSDE 0.728 OVSDE 0.696 OVSDE 0.723
OVSEE 0.746 OVSEE 0.745 OVSEE 0.682
OVSAE 0.869 OVSAE 0.825 OVSAE 0.851
VHAC 0.817 VHAC 0.824 VHAC 0.829
OVPT 0.876 OVPT 0.892 OVPT 0.887
PL.5000 0.840 PL.5000 0.814 PL.5000 0.828
PO2SM 0.869 PO2SM 0.836 PO2SM 0.745

Factor de corrector

El factor corrector, como se mencionó anteriormente, indica la proporción de información con la que contribuye el indicador simple al nuevo índice sintético, además, evita la duplicidad e incorpora información útil que retiene cada indicador simple.

Se crea un data.frame de acuerdo al factor corrector de cada indicador parcial para todos los años

correction_factors <- NULL
for(i in 1:3){
correction_factors[i] <- as.data.frame(lapply(List_DP2, function(x) get(paste0("ind_",tablas[i]))$correction_factors))
}
###Nombres de la columna 
correction_factors <- do.call(cbind.data.frame, c(Variables_sort, correction_factors)) %>%
                        subset(., select = c(1, 4, 2, 5, 3, 6))
colnames(correction_factors) <- c("Ind_2010", "2010", "Ind_2015", "2015", "Ind_2020", "2020")
Correction Factors
Ind 2010 Ind 2015 Ind 2020
ANALF 1.0000 ANALF 1.0000 ANALF 1.0000
SBASC 0.2042 OVPT 0.2076 SBASC 0.1378
OVPT 0.2139 SBASC 0.1618 OVPT 0.2053
PO2SM 0.2254 PO2SM 0.2937 OVSAE 0.1586
OVSAE 0.1433 OVSAE 0.2066 VHAC 0.2797
PL.5000 0.2420 VHAC 0.3062 PL.5000 0.2717
VHAC 0.2565 PL.5000 0.2763 PO2SM 0.4052
OVSEE 0.4599 OVSEE 0.4322 OVSDE 0.5728
OVSDE 0.5405 OVSDE 0.5950 OVSEE 0.5042

Coeficiente de Discriminación

El coeficiente de discriminación de Ivanovic mide el poder discriminante de la variable \(j\) en el conjunto de observaciones \(i\).

\[CD_{j}=\frac{2}{m\left(m-1\right)}\sum_{i,l>i}^{k_{j}}m_{ij}m_{lj}\left|\frac{x_{ij}-x_{lj}}{{\overline{X}}_{i}}\right|\]

donde:

\(\circ\:m_{ij}\): El número de observaciones de la variable \(x_{j}\)
\(\circ\:k_{j}\): El número de diferentes valores que toma \(x_{i}\) en el conjunto \(j\).

Esta medida está comprendida entre \([0, 2]\). Si una variable toma el mismo valor para todos los estados, el CD vale cero, indicando que posee un valor nulo de poder discriminante. Por el contrario, si una variable toma el valor teórico de máximo poder discriminante, el discriminante de la variable es total.

Se crea un data.frame de acuerdo al Coeficiente de discriminación (CD) de cada indicador parcial para todos los años

discrimination_coefficient <- NULL
discrimination_coefficient <- lapply(1:3, function(i) data.frame(names(get(paste0("ind_",tablas[i]))[["discrimination.coefficient"]]),
                                                                 get(paste0("ind_",tablas[i]))$discrimination.coefficient))

discrimination_coefficient <- do.call(cbind.data.frame, discrimination_coefficient)
colnames(discrimination_coefficient) <- c("Ind_2010", "2010", "Ind_2015", "2015", "Ind_2020", "2020")
Discrimination Coefficient
Ind 2010 Ind 2015 Ind 2020
ANALF 0.668 ANALF 0.699 ANALF 0.703
SBASC 0.225 SBASC 0.242 SBASC 0.276
OVSDE 0.937 OVSDE 1.016 OVSDE 1.012
OVSEE 0.769 OVSEE 0.854 OVSEE 0.842
OVSAE 0.940 OVSAE 0.883 OVSAE 0.938
VHAC 0.210 VHAC 0.243 VHAC 0.338
OVPT 0.772 OVPT 0.925 OVPT 0.906
PL.5000 0.624 PL.5000 0.624 PL.5000 0.670
PO2SM 0.344 PO2SM 0.319 PO2SM 0.164

“Cantidad de Información Global de Ivanovic Pena Relativa Individual”

\[\alpha_{i}=\frac{CD_{i}\left(1-R^{2}_{i,i-1,...,1} \right)}{\sum_{i=1}^{n}CD_{i} \left(1-R^{2}_{i,i-1,...,1} \right)}\] Esta medida, comprendida entre 0 y 1, combina la información útil y el poder discriminante de cada indicador simple y mide la cantidad de información (combinada) relativa que aporta individualmente cada indicador simple, cuando entra de forma ordenada a formar parte del indicador sintético DP2. La suma de todos los valores de \(\alpha_{i}\) es la unidad. [Zarsosa 1996, págs 158-174]

## son 9 indicadores simples
alpha <- NULL
for(i in 1:3){
alpha[[i]] <- sapply(1:length(Indicadores), function(x)(get(paste0("ind_",tablas[i]))[["correction_factors"]][x] * get(paste0("ind_",tablas[i]))[["discrimination.coefficient"]][x]) /                                                         sum(get(paste0("ind_",tablas[i]))[["correction_factors"]] * get(paste0("ind_",tablas[i]))[["discrimination.coefficient"]]))
}
Cantidad de Información Global de Ivanovic Pena Relativa Individual
Ind 2010 Ind 2015 Ind 2020
ANALF 0.344 ANALF 0.327 ANALF 0.310
SBASC 0.024 OVPT 0.024 SBASC 0.017
OVPT 0.103 SBASC 0.077 OVPT 0.092
PO2SM 0.089 PO2SM 0.117 OVSAE 0.059
OVSAE 0.069 OVSAE 0.085 VHAC 0.116
PL.5000 0.026 VHAC 0.035 PL.5000 0.040
VHAC 0.102 PL.5000 0.120 PO2SM 0.162
OVSEE 0.148 OVSEE 0.126 OVSDE 0.169
OVSDE 0.096 OVSDE 0.089 OVSEE 0.036