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
PO2SM SBASC OVPT
SBASC VHAC SBASC
VHAC OVPT VHAC
OVPT PO2SM PO2SM
OVSEE OVSEE OVSEE
OVSAE PL.5000 PL.5000
PL.5000 OVSAE OVSAE
OVSDE OVSDE OVSDE

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.844 ANALF 0.852 ANALF 0.848
SBASC 0.787 SBASC 0.784 SBASC 0.769
OVSDE 0.535 OVSDE 0.522 OVSDE 0.566
OVSEE 0.654 OVSEE 0.646 OVSEE 0.637
OVSAE 0.631 OVSAE 0.594 OVSAE 0.582
OVPT 0.676 OVPT 0.741 OVPT 0.772
VHAC 0.708 VHAC 0.746 VHAC 0.731
PL.5000 0.602 PL.5000 0.597 PL.5000 0.587
PO2SM 0.794 PO2SM 0.678 PO2SM 0.699

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
PO2SM 0.4767 SBASC 0.4083 OVPT 0.5840
SBASC 0.3014 VHAC 0.5778 SBASC 0.3895
VHAC 0.5572 OVPT 0.5484 VHAC 0.5041
OVPT 0.6520 PO2SM 0.5099 PO2SM 0.4435
OVSEE 0.7916 OVSEE 0.7370 OVSEE 0.7405
OVSAE 0.7086 PL.5000 0.6009 PL.5000 0.5678
PL.5000 0.5245 OVSAE 0.7419 OVSAE 0.7915
OVSDE 0.7524 OVSDE 0.7574 OVSDE 0.6382

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.743 ANALF 0.781 ANALF 0.792
SBASC 0.278 SBASC 0.310 SBASC 0.348
OVSDE 1.170 OVSDE 1.271 OVSDE 1.251
OVSEE 1.146 OVSEE 1.181 OVSEE 1.159
OVSAE 1.200 OVSAE 1.211 OVSAE 1.281
OVPT 0.961 OVPT 1.063 OVPT 1.093
VHAC 0.316 VHAC 0.352 VHAC 0.443
PL.5000 0.504 PL.5000 0.504 PL.5000 0.537
PO2SM 0.362 PO2SM 0.351 PO2SM 0.159

“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.178 ANALF 0.172 ANALF 0.188
PO2SM 0.032 SBASC 0.028 OVPT 0.048
SBASC 0.085 VHAC 0.162 SBASC 0.115
VHAC 0.153 OVPT 0.143 VHAC 0.138
OVPT 0.188 PO2SM 0.136 PO2SM 0.135
OVSEE 0.182 OVSEE 0.172 OVSEE 0.192
OVSAE 0.054 PL.5000 0.047 PL.5000 0.060
PL.5000 0.063 OVSAE 0.082 OVSAE 0.101
OVSDE 0.065 OVSDE 0.059 OVSDE 0.024