Segmentación de países según infraestructura y tecnología: Análisis de Componentes Principales y clusterización de Ward

Segmentación de países por infraestructura y tecnología con ACP y clusterización de Ward

Universidad del Valle · Facultad de Ingeniería · Gestión de Datos

Aprendizaje no supervisado (Taller 3) · Profesor: Orlando Joaqui Barandica

Banco Mundial (WDI) · 114 países · año 2019 · 17 variables cuantitativas

Integrantes:Miguel Fabricio Arguello Cedeño, Juan Camilo Alvarado Gutierrez, Santiago Grijalba Hernandez, Jose Manuel Toro Rodriguez · 7 de noviembre de 2026

114
Países con datos completos
17
Variables cuantitativas
2019
Año de análisis
74,4 %
Varianza explicada por 3 componentes
5
Clusters de países (Ward)

Resumen

Este trabajo caracteriza y segmenta a 114 países a partir de 17 indicadores de los World Development Indicators del Banco Mundial para 2019, agrupados en cuatro dimensiones: energía, conectividad digital, transporte aéreo y economía y estructura del país. Como los indicadores están muy relacionados entre sí, se aplicó primero un Análisis de Componentes Principales (ACP) sobre las variables estandarizadas (con logaritmo decimal en 10 variables muy asimétricas) y después el método jerárquico de Ward sobre los puntajes factoriales.

El ACP retuvo 3 componentes que explican 74,4 % de la variabilidad: un eje de desarrollo y conectividad (CP1, 52,6 %), un eje de tamaño del país (CP2, 14,2 %) y un eje de perfil energético y tecnológico (CP3, 7,6 %). La clusterización de Ward produjo 5 grupos (silueta promedio de 0,339): infraestructura incipiente, emergentes grandes y medianas, pequeñas de conectividad media-alta, alto consumo energético (hidrocarburos) y economías avanzadas. Los grupos se ordenan según el nivel de desarrollo y solo en parte coinciden con los continentes, lo que indica que la similitud en infraestructura y tecnología depende más del perfil de la economía que de su ubicación geográfica.

Palabras clave: Análisis de Componentes PrincipalesMétodo de WardInfraestructuraConectividad digitalBanco MundialDesarrollo

Cómo leer este reporte. Los gráficos son interactivos: pase el cursor sobre un punto para ver el país y sus valores, haga clic en la leyenda para ocultar o mostrar grupos y arrastre para ampliar una zona. Las tablas se pueden ordenar, buscar y descargar (Copiar, CSV, Excel). El botón ↑ de la esquina inferior derecha lleva al inicio y el botón ☰ despliega o contrae los apartados desplegables.

Introducción

Contexto

La infraestructura energética, digital y de transporte de un país condiciona su capacidad de producir, de conectarse con el resto del mundo y de ofrecer servicios básicos a su población. El acceso a la electricidad, la penetración de internet y de la banda ancha, la intensidad del transporte aéreo o la capacidad de exportar bienes de alta tecnología son, en conjunto, una radiografía del nivel de desarrollo y de la brecha tecnológica entre países. Por esta razón, estos temas ocupan un lugar central en la Agenda 2030 (Naciones Unidas, 2015), en particular en los objetivos de energía asequible y no contaminante (ODS 7) y de industria, innovación e infraestructura (ODS 9), y su comparación entre países es un insumo habitual de la economía del desarrollo y de la política pública.

Las diferencias son muy grandes. En la muestra de este estudio, el porcentaje de la población que usa internet va de 9,3 % (Burkina Faso) a 99,7 % (Bahrain), y el consumo de electricidad por habitante pasa de 56 kWh (Rwanda) a 21.814 kWh (Bahrain), es decir, 388 veces más. Sin embargo, estas brechas no se distribuyen igual en todas las dimensiones ni coinciden siempre con el nivel de ingreso o con la región geográfica: un país puede tener un ingreso alto y una matriz eléctrica casi totalmente fósil, o un ingreso medio y una conectividad digital elevada.

Los Indicadores del Desarrollo Mundial (World Development Indicators, WDI) del Banco Mundial reúnen, con criterios comparables, cientos de series para más de doscientas economías (Banco Mundial, s. f.). En este trabajo se emplean los datos de 2019 de 217 países y, tras conservar únicamente aquellos con información completa, de 114 países descritos por 17 indicadores agrupados en cuatro dimensiones: energía, conectividad digital, transporte aéreo y economía y estructura del país. Se eligió 2019 por ser el último año anterior a la pandemia de COVID-19, que alteró de forma abrupta indicadores como el transporte aéreo.

Justificación y problema a modelar

Comparar países con una sola variable, como el PIB per cápita, oculta que el desarrollo de la infraestructura es multidimensional. Mirar los 17 indicadores uno por uno tampoco es práctico, porque además están muy relacionados entre sí: el valor absoluto medio de sus correlaciones es de 0,46 y 42 % de los pares supera 0,5. Esa redundancia hace que la información útil esté repartida y repetida entre muchas variables.

El problema es de aprendizaje no supervisado (James et al., 2021): no existe una variable respuesta que predecir, sino un conjunto de indicadores que se quiere resumir y usar para descubrir estructura. Por eso se aplica, en primer lugar, un Análisis de Componentes Principales (ACP) que resume los 17 indicadores en pocas dimensiones no correlacionadas (Jolliffe y Cadima, 2016) y, en segundo lugar, el método jerárquico de Ward sobre los factores estimados (Ward, 1963), que agrupa a los países con perfiles similares. Aplicar la clusterización sobre los componentes y no sobre las variables originales es una estrategia habitual en el análisis multivariado exploratorio (Husson et al., 2017), porque evita que la información redundante pese varias veces en las distancias.

La utilidad de este tipo de análisis es múltiple. El ACP permite construir ejes sintéticos (por ejemplo, un eje de desarrollo y conectividad) que facilitan comparar países y seguir su evolución. La clusterización entrega una tipología que sirve para definir grupos de referencia o benchmarking y para decidir qué política de infraestructura o de digitalización conviene a cada grupo. Y los patrones que aparecen (regiones similares, niveles de desarrollo, perfiles energéticos) ayudan a formular hipótesis sobre las trayectorias de los países. Cabe señalar que la segmentación es descriptiva: agrupa países parecidos, pero no explica por qué lo son ni demuestra causalidad.

Pregunta de investigación

¿Es posible resumir en pocas dimensiones la información de 17 indicadores de infraestructura, tecnología y economía de los países y, a partir de ellas, identificar grupos de países con perfiles similares?

Esta pregunta se desglosa en las tres que plantea la guía del taller: (i) ¿qué componentes principales explican la mayor variabilidad y qué representan?; (ii) ¿cuántos clusters se forman y qué características define a cada grupo?; y (iii) ¿qué patrones aparecen entre los países (regiones, niveles de desarrollo, perfil energético)?

Objetivo general

Aplicar el Análisis de Componentes Principales y la clusterización de Ward para caracterizar y segmentar a los países según su infraestructura, su conectividad tecnológica y su estructura económica en 2019.

Objetivos específicos

  1. Construir la base de análisis a partir de los WDI (solo países, año 2019, 17 variables cuantitativas) y describir sus variables y su relación con el continente.
  2. Examinar las relaciones entre las variables mediante la matriz de correlaciones y verificar la adecuación de los datos para el ACP.
  3. Estimar los componentes principales, evaluar cuánta variabilidad explica cada uno e interpretarlos.
  4. Agrupar a los países con el método de Ward a partir de los factores estimados y determinar el número adecuado de clusters.
  5. Caracterizar cada cluster e identificar patrones entre países (regiones, niveles de desarrollo, perfil energético).
  6. Discutir posibles usos del análisis en economía y políticas públicas, y las limitaciones de los resultados.

Hipótesis de trabajo

Antes de estimar los modelos se plantearon tres expectativas, que se contrastan en los resultados y se retoman en las conclusiones:

  • H1. Por la fuerte relación entre los indicadores, unos pocos componentes principales concentrarán la mayor parte de la variabilidad.
  • H2. El primer componente reflejará el nivel de desarrollo y conectividad (ingreso, energía, internet), y otro componente distinto recogerá el tamaño de los países.
  • H3. Los clusters se ordenarán según el nivel de desarrollo y no coincidirán por completo con los continentes: habrá grupos que mezclen regiones.

Alcance y organización del reporte

El análisis es transversal (un solo año), cubre los 114 países con datos completos y es descriptivo. El reporte sigue la estructura pedida en la guía: la Metodología describe la construcción de la base, las variables y los modelos; los Resultados descriptivos presentan tablas y gráficos dinámicos de las variables; los Resultados del modelo reúnen el ACP y la clusterización de Ward con su interpretación; y el cierre contiene las Conclusiones, las Limitaciones, las Recomendaciones y la Bibliografía.

Metodología

El estudio adoptó un enfoque de aprendizaje no supervisado con dos técnicas encadenadas: el Análisis de Componentes Principales (ACP) para reducir la dimensionalidad y el método jerárquico de Ward para segmentar a los países. El proceso se organizó en seis etapas, que se resumen a continuación y se detallan en las secciones siguientes.

DatosWDI 2019: 217 países y 19 indicadores
Depuración114 países con datos completos y 17 variables
Transformaciónlog10 en 10 variables y estandarización
ACP3 componentes (74 % de la varianza)
Ward5 clusters sobre los puntajes factoriales
InterpretaciónPerfiles de los grupos, patrones y usos

Construcción de la base de datos

Fuente de datos: World Development Indicators (WDI) del Banco Mundial, descargados desde DataBank (https://databank.worldbank.org/source/world-development-indicators).

Unidad de análisis: el país, para el año 2019. Se filtraron únicamente países, sin agregados regionales ni bloques económicos (por ejemplo, “Mundo”, “Unión Europea” o “Ingreso alto”): la base descargada contiene 217 países y 19 indicadores.

Criterio de selección del año. Se eligió 2019 por tres razones: (i) es el último año previo a la pandemia de COVID-19, cuyos efectos en 2020 y 2021 (confinamientos, cierre de fronteras, caída del transporte aéreo) distorsionan la comparación entre países; (ii) la mayoría de los indicadores seleccionados tiene en ese año una cobertura razonable (Tabla 1); y (iii) un solo año evita mezclar valores de momentos distintos, de modo que cada país se describe por un único perfil.

Criterio de selección de las variables. Se descargaron 19 indicadores relacionados con cuatro dimensiones de la infraestructura y la tecnología: energía (acceso a electricidad, consumo de electricidad, uso de energía y participación de la electricidad renovable), conectividad digital (usuarios de internet, servidores seguros, banda ancha fija, telefonía celular y telefonía fija), transporte aéreo (pasajeros y salidas de aerolíneas) y economía y estructura del país (PIB per cápita y total, exportaciones de alta tecnología, población, urbanización y superficie). Para decidir qué indicadores conservar se aplicaron cuatro criterios: que fueran cuantitativos y continuos, que fueran comparables entre países (definiciones y unidades homogéneas del Banco Mundial), que tuvieran cobertura suficiente y que describieran aspectos distintos de la infraestructura. Se incluyen las variables de tamaño porque la infraestructura de un país depende tanto de su nivel de desarrollo como de su escala. Dos indicadores se descartaron por tener muy pocos datos: los kilómetros de ferrocarril (74 países) y el gasto en investigación y desarrollo (98 países). Quedaron así 17 variables cuantitativas, por encima del mínimo de 15 que exige la guía.

Depuración. El ACP exige observaciones completas, por lo que se conservaron solo los países con información en todas las variables del análisis, sin imputar valores. La base final quedó con 114 países (52,5 % de los 217 descargados). Los 103 países excluidos son, en su mayoría, economías pequeñas: su población mediana es de 0,9 millones de habitantes, frente a 14,6 millones en los países retenidos. Por tanto, los resultados describen a los países con información completa y no a todo el universo mundial.

Variables derivadas. Los totales del transporte aéreo (pasajeros y salidas) dependen del tamaño del país, por lo que se expresaron por habitante (pasajeros por habitante y salidas por millón de habitantes). Así, esas variables miden la intensidad del transporte aéreo y no el tamaño del país, que ya está representado por la población, el PIB total y la superficie.

Transformación logarítmica. Varias variables presentan una asimetría fuerte a la derecha, que hace que unos pocos países extremos dominen la varianza. Se midió la asimetría de cada variable con el coeficiente

\[g_1 = \frac{m_3}{m_2^{3/2}}, \qquad m_k = \frac{1}{n}\sum_{i=1}^{n}(x_i - \bar{x})^k\]

y se aplicó el logaritmo decimal a las 10 variables con \(g_1 > 1\): Consumo eléctrico pc, Uso de energía pc, Servidores seguros, Pasajeros aéreos pc, Salidas aéreas / millón, PIB pc PPA, PIB total, Exp. alta tecnología, Población y Superficie. Todas toman valores estrictamente positivos (mínimo observado de 0,0005), por lo que el logaritmo está definido sin sumar constantes. Las demás variables se dejaron en su escala original.

Descripción de las variables

Las variables se agrupan en cuatro dimensiones que sirven de guía para interpretar los resultados:

  • Energía. Mide la cobertura y la intensidad del sistema eléctrico. El acceso a electricidad es una condición básica de desarrollo; el consumo de electricidad y el uso de energía por habitante reflejan la intensidad de la actividad productiva y del consumo; la participación de renovables describe la composición de la matriz eléctrica y no necesariamente el nivel de desarrollo.
  • Conectividad digital. Los usuarios de internet, la banda ancha fija y los servidores seguros miden la adopción y la calidad de la infraestructura digital; la telefonía celular mide la cobertura móvil, y la telefonía fija, una tecnología más antigua.
  • Transporte aéreo. Los pasajeros y las salidas por habitante aproximan la integración del país con el resto del mundo y su capacidad logística.
  • Economía y estructura. El PIB per cápita (en paridad de poder adquisitivo) es una medida del ingreso medio; las exportaciones de alta tecnología reflejan la sofisticación productiva; el PIB total, la población y la superficie miden el tamaño; y la urbanización, el grado de concentración de la población en ciudades.

Preprocesamiento de variables

Las variables están expresadas en unidades muy distintas (porcentajes, kWh, miles de millones de dólares, suscripciones por cada 100 habitantes). Como el ACP sobre la matriz de covarianzas daría un peso desproporcionado a las variables de mayor escala, las 17 variables (ya transformadas con logaritmo cuando correspondía) se estandarizaron restando la media y dividiendo por la desviación estándar:

\[z_{ij} = \frac{x_{ij} - \bar{x}_j}{s_j}\]

Con esta estandarización el ACP equivale a trabajar con la matriz de correlaciones y todas las variables tienen el mismo peso inicial.

Modelos utilizados

Análisis de Componentes Principales (ACP)

El ACP busca combinaciones lineales de las variables originales, llamadas componentes principales, que sean no correlacionadas entre sí y que capturen la mayor variabilidad posible (Jolliffe y Cadima, 2016). Se obtienen de la descomposición espectral de la matriz de correlaciones \(R\):

\[R\, v_k = \lambda_k\, v_k, \qquad C_k = Z\, v_k\]

donde \(\lambda_k\) es el autovalor del componente \(k\) (la varianza que explica), \(v_k\) su autovector y \(C_k\) los puntajes de los países en ese componente. La proporción de varianza explicada es \(\lambda_k / \sum_j \lambda_j\). Para interpretar cada componente se usan las cargas \(r_{jk} = v_{jk}\sqrt{\lambda_k}\), que son las correlaciones entre cada variable y el componente; se consideran relevantes las cargas con valor absoluto de 0,5 o más (Hair et al., 2019). Antes de extraer los componentes se verifica la adecuación de los datos con la medida de Kaiser-Meyer-Olkin (KMO; Kaiser, 1974) y la prueba de esfericidad de Bartlett (Bartlett, 1954).

El número de componentes a retener se decide con tres criterios: el criterio de Kaiser (autovalor mayor que 1; Kaiser, 1960), la varianza acumulada (se busca una referencia del orden de 70 %) y el gráfico de sedimentación. Como el signo de un componente es arbitrario, se fijó una convención: el CP1 es positivo cuando aumenta el PIB per cápita, el CP2 cuando aumenta la población y el CP3 cuando aumenta la participación de electricidad renovable.

Clusterización jerárquica de Ward

Con los puntajes de los países en los componentes retenidos se aplicó el método de Ward, un algoritmo jerárquico aglomerativo que parte de cada país como un grupo y, en cada paso, une los dos grupos cuya fusión produce el menor aumento de la suma de cuadrados dentro de los grupos (Ward, 1963):

\[\Delta(A, B) = \frac{n_A\, n_B}{n_A + n_B}\, \lVert \bar{c}_A - \bar{c}_B \rVert^2\]

donde \(n_A\) y \(n_B\) son los tamaños de los grupos y \(\bar{c}_A\), \(\bar{c}_B\) sus centros en el espacio de los componentes. Se usó la distancia euclidiana y la implementación ward.D2 de R, que corresponde al criterio original de Ward (Murtagh y Legendre, 2014). Trabajar con los puntajes factoriales, y no con las variables originales, evita que la información redundante entre indicadores correlacionados pese varias veces en las distancias.

El número de clusters se determinó combinando tres elementos: el dendrograma y las alturas de fusión, el ancho de silueta promedio (Rousseeuw, 1987) y la interpretabilidad y el tamaño de los grupos resultantes. El ancho de silueta de un país es

\[s(i) = \frac{b(i) - a(i)}{\max\{a(i),\, b(i)\}}\]

donde \(a(i)\) es la distancia media al resto de su cluster y \(b(i)\) la distancia media al cluster más cercano; valores cercanos a 1 indican países bien asignados y valores negativos, países que quizá pertenecen a otro grupo. Según Kaufman y Rousseeuw (1990), un promedio por encima de 0,5 indica una estructura razonable y entre 0,25 y 0,5, una estructura débil. Una vez fijado el número de clusters, estos se caracterizaron con las medianas de las variables en sus unidades originales, con las medias estandarizadas (para ver qué tan por encima o por debajo del promedio está cada grupo) y con su composición por continente. Los clusters se numeraron de menor a mayor CP1 medio.

Supuestos y alcance de los métodos

  • El ACP capta relaciones lineales entre variables y es sensible a las escalas y a los valores extremos, de ahí la estandarización y la transformación logarítmica.
  • El método de Ward tiende a formar grupos compactos y de tamaño comparable, y sus fusiones son irreversibles: una vez que dos países se unen, no se separan en pasos posteriores.
  • Ambos son métodos exploratorios y descriptivos: no contrastan hipótesis causales ni predicen, sino que revelan estructura en los datos.

Resumen de las decisiones metodológicas

Herramientas y reproducibilidad

El análisis se realizó con los paquetes tidyverse, readxl, plotly, DT y cluster. El ACP se estimó con la función prcomp() y la clusterización con hclust() (método ward.D2). Ninguno de los dos procedimientos usa números aleatorios, por lo que el código reproduce exactamente los mismos resultados cada vez que se ejecuta con el mismo archivo de datos.

Resultados descriptivos

Esta sección describe las 17 variables de los 114 países de la base de análisis (las mismas del ACP), en sus unidades originales y antes de aplicar la transformación logarítmica.

Distribución de países por continente

La base está compuesta principalmente por países de Asia (38 países, 33,3 %) y Europa (32, 28,1 %); Oceanía tiene solo 2. Por eso las comparaciones por continente de las secciones siguientes se limitan a los continentes con al menos 10 países.

Resumen estadístico de las variables

Los valores de la Tabla 4 muestran escalas muy distintas entre variables, lo que justifica la estandarización. En varias de ellas la media supera con claridad a la mediana (por ejemplo, el PIB total y la población), señal de distribuciones asimétricas a la derecha; la última columna recoge el coeficiente de asimetría con el que se decidió aplicar el logaritmo a 10 variables.

Efecto de la transformación logarítmica

Figura 1. Histograma del consumo eléctrico per cápita (escala original).

Figura 2. Histograma del consumo eléctrico per cápita (escala log10).

En escala original (Figura 1) la mayoría de los países se concentra en valores bajos y unos pocos, de consumo muy alto, alargan la cola derecha (asimetría de 1,55). Tras el logaritmo (Figura 2) la distribución es mucho más simétrica, de modo que ningún país extremo domina la varianza del ACP.

Medianas por continente

Entre los continentes con al menos 10 países, la mediana de usuarios de internet es más alta en Europa (82,9 %) y más baja en África (30,9 %); en el PIB per cápita la mediana más alta es la de Europa (41.480 USD) y la más baja la de África (5.102 USD). El consumo eléctrico per cápita sigue el mismo orden: Europa tiene la mediana más alta (5.328 kWh) y África la más baja (361 kWh).

Variables estandarizadas

Figura 3. Distribución de las variables estandarizadas (diagramas de caja).

La Figura 3 muestra, en escala original, cuántos países se alejan más de dos desviaciones estándar del promedio en cada variable. Las variables de tamaño (población, PIB total y superficie) y los servidores seguros presentan los valores atípicos más extremos, lo que confirma la necesidad de la transformación logarítmica antes del ACP.

Resultados del modelo

Análisis de Componentes Principales

En esta sección se aplica el ACP a las 17 variables de 114 países, ya transformadas con logaritmo decimal en 10 de ellas y estandarizadas, tal como se describió en la metodología. El análisis sigue el orden habitual de un ACP exploratorio (Husson et al., 2017): (i) examinar la matriz de correlaciones, (ii) verificar que los datos son adecuados, (iii) decidir cuántos componentes retener, (iv) interpretar los componentes a partir de sus cargas y de las contribuciones, y (v) representar a las variables y a los países en los planos factoriales.

0,84
KMO global (adecuación del ACP)
3
Componentes retenidos (autovalor > 1)
52,6 %
Varianza explicada por el CP1
74,4 %
Varianza acumulada con 3 componentes
74,4 %
Comunalidad media de las variables

Matriz de correlaciones

Figura 4. Matriz de correlaciones de Pearson entre las 17 variables (variables transformadas, ordenadas por dimensión).

La Figura 4 muestra las correlaciones entre las 17 variables, que son altas en su conjunto: el valor absoluto medio es de 0,46 y 42 % de los 136 pares supera 0,5 en valor absoluto. Esto indica que hay mucha información compartida y que tiene sentido resumirla en pocos componentes. Al mismo tiempo, 38 pares tienen correlación negativa, lo que anticipa que no todas las variables se alinean en un único eje.

Los pares más correlacionados (Tabla 6) son Pasajeros aéreos pc con Salidas aéreas / millón (r = 0,97), Consumo eléctrico pc con Usuarios internet (r = 0,92) y Consumo eléctrico pc con Uso de energía pc (r = 0,92). El primero es esperable, porque ambas variables miden la intensidad del transporte aéreo en relación con la población. Los siguientes muestran que el PIB per cápita, el consumo eléctrico, el uso de energía y el uso de internet forman un bloque muy cohesionado, con correlaciones entre 0,87 y 0,92: los países con mayores ingresos por habitante consumen más energía y están más conectados.

Un segundo grupo de variables se comporta de forma distinta. La población, la superficie y el PIB total (las variables de tamaño) se relacionan entre sí (por ejemplo, población y PIB total, r = 0,70), pero tienen correlaciones bajas, o incluso negativas, con las variables de nivel de vida: la población se relaciona con el PIB per cápita con r = -0,29, la superficie con r = -0,30 y la electricidad renovable con r = -0,33. Las variables que se mueven en sentido contrario al PIB per cápita son Electricidad renovable, Población y Superficie. Por último, las variables menos relacionadas con el conjunto son Exp. alta tecnología, Electricidad renovable y Superficie (correlación absoluta media de 0,23, 0,23 y 0,25), de modo que aportan información propia y no repetida. Esto anticipa que el tamaño de los países y la matriz renovable se separarán del eje principal de desarrollo.

Adecuación de los datos para el ACP

Adecuación global. La medida KMO global es de 0,840, que según Kaiser y Rice (1974) se califica de meritoria (valores superiores a 0,6 se consideran aceptables). La prueba de esfericidad de Bartlett (1954) rechaza que la matriz de correlaciones sea la identidad (chi-cuadrado = 2.452,6; 136 g.l.; p < 0,001). Ambas pruebas indican que los datos son adecuados para aplicar el ACP.

Adecuación por variable. El KMO individual va de 0,437 (Población) a 0,960 (Celulares). Solo Población (0,437) queda por debajo de 0,5, el valor mínimo que suele recomendarse para cada variable (Hair et al., 2019), y PIB total (0,575) se ubica entre 0,5 y 0,6. Se trata de variables de tamaño: se relacionan entre sí pero poco con el resto, y por eso su adecuación individual es menor. Se mantienen en el análisis porque describen una dimensión propia (el CP2) y porque la adecuación global es alta; a cambio, hay que tener presente que el CP2 es un eje de tamaño que se separa del resto de los indicadores.

Varianza explicada y número de componentes

Figura 5. Gráfico de sedimentación: autovalores y varianza acumulada de los componentes.

El primer componente explica 52,6 % de la variabilidad total, el segundo 14,2 % y el tercero 7,6 %. El autovalor del CP1 (8,9) equivale a la varianza de casi 9 de las 17 variables estandarizadas, lo que muestra qué tan redundante es el conjunto de indicadores. Para decidir cuántos componentes retener se usaron los tres criterios de la metodología:

  • Criterio de Kaiser (autovalor mayor que 1; Kaiser, 1960): retiene 3 componentes, con autovalores de 8,95, 2,41 y 1,30.
  • Varianza acumulada: con 2 componentes se alcanza 66,8 %, por debajo de la referencia de 70 %, y con 3 se llega a 74,4 %, que la supera. Es decir, con 3 dimensiones se conserva cerca de tres cuartas partes de la información de las 17 variables originales.
  • Gráfico de sedimentación (Cattell, 1966): hay una caída muy fuerte entre el CP1 y el CP2, y a partir del CP4 la curva se aplana, con autovalores por debajo de 1 (Figura 5).

El autovalor del CP4 (0,90) queda por debajo de 1 y aporta 5,3 % adicional; se descarta porque no alcanza el umbral de Kaiser, y porque los 3 primeros componentes ya tienen una interpretación clara, como se verá a continuación. Como señalan Jolliffe y Cadima (2016), la decisión sobre el número de componentes combina criterios estadísticos con el juicio del analista, de modo que otra elección razonable podría ser distinta.

Cargas e interpretación de los componentes

Las cargas son las correlaciones entre cada variable y cada componente (Jolliffe y Cadima, 2016), y se consideran relevantes las de 0,5 o más en valor absoluto (Hair et al., 2019). La última columna de la Tabla 9 es la comunalidad: la proporción de la varianza de cada variable que recogen en conjunto los 3 componentes retenidos, \(h_j^2 = \sum_{k=1}^{3} r_{jk}^2\).

Figura 6. Mapa de calor de las cargas de las variables en los componentes retenidos.

CP1: nivel de desarrollo y conectividad (52,6 % de la varianza). 12 de las 17 variables tienen una carga de 0,5 o más en valor absoluto, todas con signo positivo. Las más altas son las del PIB per cápita (0,97), el consumo eléctrico (0,94), los usuarios de internet (0,94), el uso de energía (0,88), los servidores seguros (0,86) y la banda ancha fija (0,84), seguidas de las dos variables de transporte aéreo (0,81 y 0,80), la telefonía fija (0,79), la urbanización (0,79) y el acceso a electricidad (0,76). La telefonía celular tiene una carga más moderada (0,54), porque su cobertura está casi saturada en buena parte de los países. Las únicas variables con signo negativo son Electricidad renovable (-0,37), Superficie (-0,28) y Población (-0,25), todas con cargas débiles (menores que 0,5 en valor absoluto). En el caso de la electricidad renovable, su relación negativa con el ingreso (r = -0,33 con el PIB per cápita) indica que una matriz renovable no es sinónimo de mayor desarrollo; una explicación posible, que estos datos no permiten comprobar, es que en los países de menor ingreso la matriz eléctrica, pequeña, descanse más en la hidroelectricidad.

Este componente es, por tanto, un eje de desarrollo y conectividad: reúne ingreso, energía, tecnologías digitales y transporte. Un país con un CP1 alto es rico en infraestructura y tecnología; uno con un CP1 bajo es pobre en ella. Los países con los puntajes más altos son Hong Kong SAR, China, Singapore, Luxembourg, United Arab Emirates, Netherlands, y los más bajos son Uganda, Madagascar, Burkina Faso, Mozambique, Ethiopia.

CP2: tamaño del país (14,2 %). Este componente es casi independiente del nivel de desarrollo. Está definido por la población (0,92), el PIB total (0,83) y la superficie (0,81); las exportaciones de alta tecnología tienen una carga más débil (0,37) y las demás variables tienen cargas cercanas a cero (a lo sumo 0,20 en valor absoluto). El PIB total también carga en el CP1 (0,49), porque depende a la vez del tamaño del país y de su ingreso. Es un eje de escala y tamaño: separa a los países grandes en población, territorio y economía total de los pequeños, sin importar si son ricos o pobres. Los países con puntajes más altos son China, United States, India, Brazil, Russian Federation, y los más bajos son Malta, Montenegro, Eswatini, Mauritius, Trinidad and Tobago.

CP3: renovables y exportaciones de alta tecnología (7,6 %). Es el componente más difícil de leer, porque explica menos varianza y reúne dos rasgos que se relacionan muy poco entre sí (r = 0,05). Las cargas positivas más altas son las de la electricidad renovable (0,63) y las exportaciones de alta tecnología (0,63), seguidas de la banda ancha fija (0,38) y los servidores seguros (0,26); las cargas negativas son débiles (por ejemplo, telefonía celular, -0,22, y urbanización, -0,19). El componente separa, dentro de un mismo nivel de desarrollo, a los países con alguna de estas dos características (una matriz eléctrica renovable o una oferta exportadora tecnológica) de aquellos que no tienen ninguna. Los países con los puntajes más altos son Lithuania, Luxembourg, Lao PDR, Uganda, Estonia, y los más bajos son Qatar, Kuwait, Saudi Arabia, Libya, Botswana: estos últimos tienen una participación renovable de a lo sumo 0,25 % y exportaciones de alta tecnología de a lo sumo 0,9 % de las manufacturas.

Calidad de la representación. Con 3 componentes la comunalidad media es de 74,4 %. Las variables mejor representadas son Población (0,94), PIB pc PPA (0,93) y PIB total (0,93), y las peor representadas son Celulares (0,38), Electricidad renovable (0,53) y Acceso electricidad (0,58). Parte de la información de estas últimas queda en componentes descartados: por ejemplo, la telefonía celular tiene una carga de 0,43 y la electricidad renovable de 0,45 (en valor absoluto) en el CP4. Es el costo de resumir 17 indicadores en 3 dimensiones.

Contribución de las variables y de los países

La contribución mide qué parte de cada componente se debe a una variable o a un país concreto (Abdi y Williams, 2010; Husson et al., 2017). Para la variable \(j\) y para el país \(i\) en el componente \(k\) se calcula, respectivamente,

\[\text{Contr}_{jk} = 100\,\frac{r_{jk}^2}{\lambda_k}, \qquad \text{Contr}_{ik} = 100\,\frac{c_{ik}^2}{(n-1)\,\lambda_k}\]

donde \(r_{jk}\) es la carga, \(c_{ik}\) el puntaje del país y \(\lambda_k\) el autovalor. En cada componente las contribuciones suman 100 %.

En el CP1 las contribuciones están repartidas entre muchas variables; las mayores son las de PIB pc PPA (10,4 %), Usuarios internet (10,0 %) y Consumo eléctrico pc (9,8 %). En el CP2 pesan sobre todo Población (35,4 %), PIB total (28,3 %) y Superficie (27,5 %), y en el CP3, Electricidad renovable (30,5 %), Exp. alta tecnología (30,1 %) y Banda ancha fija (11,1 %). Esto confirma la lectura de los componentes: el CP1 resume un conjunto amplio de indicadores de nivel, mientras que el CP2 y el CP3 dependen cada uno de muy pocas variables.

En el CP1 y en el CP2 la contribución está repartida: el país que más aporta es Uganda en el CP1 (4,7 %) y China en el CP2 (7,0 %). En el CP3, en cambio, un solo país, Qatar, aporta 15,7 % de la inercia del componente (Tabla 11). Por eso el CP3 es el más sensible a los valores extremos y su interpretación debe tomarse con más cautela que la de los otros dos.

Representación de las variables y de los países

Figura 7. Círculo de correlaciones: cargas de las variables en CP1 y CP2.

La Figura 7 resume la interpretación. Casi todas las variables de energía, conectividad y transporte aéreo apuntan hacia la derecha (carga alta en CP1), con una carga cercana a cero en CP2, mientras que la población, la superficie y el PIB total apuntan hacia arriba (carga alta en CP2). El PIB total es la variable mixta: su flecha apunta hacia arriba y hacia la derecha, porque depende del tamaño del país y de su ingreso. Los colores muestran que las dimensiones de energía, conectividad y transporte aéreo se alinean sobre el mismo eje, y que la dimensión de economía y estructura se divide en dos: el ingreso y la urbanización acompañan al CP1, y las variables de tamaño forman el CP2. Las flechas de la electricidad renovable, la población y la superficie apuntan algo hacia la izquierda, es decir, en sentido contrario al nivel de desarrollo, aunque con cargas débiles en el CP1. Cuanto más largo es el vector (más cerca del círculo), mejor representada está la variable en este plano.

Figura 8. Países en los planos CP1–CP2 (tamaño) y CP1–CP3 (renovables y alta tecnología), según el continente.

La Figura 8 muestra a los países sobre los ejes de los componentes, coloreados por continente. El eje horizontal (CP1) ordena a los países de menor a mayor nivel de desarrollo y conectividad: a la izquierda quedan casi todos los países africanos (media de CP1 de -3,78) y a la derecha los europeos (2,22) y los dos de Oceanía (3,02), con Asia (0,34) y América (0,04) en posiciones intermedias y muy dispersas. En el panel izquierdo, el eje vertical (CP2) separa a los países por su tamaño y no por el continente: arriba están los grandes (China, United States, India, Brazil, Russian Federation) y abajo los pequeños (Malta, Montenegro, Eswatini, Mauritius, Trinidad and Tobago). Esto confirma que el continente se relaciona con el CP1, pero que el CP2 recoge una diferencia de escala que atraviesa todos los continentes. Como referencia, Colombia queda cerca del centro del eje de desarrollo (CP1 = -0,14) y por encima del promedio en tamaño (CP2 = 1,21).

En el panel derecho (CP3), la media por continente es más alta en Europa (0,78) y más baja en Asia (-0,53). Los países con CP3 muy bajo (Qatar, Kuwait, Saudi Arabia, Libya, Botswana) se ven en la parte inferior del panel, lo que coincide con su peso en este componente (Tabla 11).

Figura 9. Países con mayor y menor puntaje en el CP1 (nivel de desarrollo y conectividad).

La Figura 9 pone nombre a los extremos del eje de desarrollo. En la parte superior aparecen economías de alto ingreso e infraestructura avanzada, como Hong Kong SAR, China, Singapore, Luxembourg, y en la inferior países de baja infraestructura, como Uganda, Madagascar, Burkina Faso: 12 de los 12 países con menor puntaje son africanos. Los puntajes del CP1 van de -6,9 a 5,2, con un promedio de cero por construcción, lo que refleja que este componente concentra la mayor parte de las diferencias entre los países.

Síntesis del ACP

¿Qué componentes explican mayor variabilidad y qué representan? El ACP resumió las 17 variables en 3 componentes que acumulan 74,4 % de la variabilidad. El CP1 (52,6 %) representa el nivel de desarrollo y conectividad: ingreso por habitante, energía, internet, banda ancha y transporte aéreo. El CP2 (14,2 %) representa el tamaño del país: población, PIB total y superficie, y es independiente del desarrollo. El CP3 (7,6 %) distingue el perfil energético y tecnológico: participación de renovables y exportaciones de alta tecnología.

Resultados a tener en cuenta. Los datos son adecuados para el ACP (KMO de 0,840 y prueba de Bartlett con p < 0,001), aunque Población (0,437) tiene un KMO individual inferior a 0,5. El CP3 depende mucho de pocos países (Qatar aporta 15,7 % de su inercia) y las variables peor representadas, como Celulares y Electricidad renovable, pierden parte de su información al pasar de 17 indicadores a 3 componentes.

Clusterización de Ward

Con los 3 componentes retenidos del ACP se aplicó el método jerárquico de Ward, tal como se describió en la metodología: distancia euclidiana entre los puntajes factoriales de los países y enlace ward.D2. Primero se examina el dendrograma y se justifica el número de clusters; después se caracteriza cada grupo y se revisa su relación con los continentes.

Dendrograma y número de clusters

Figura 10. Dendrograma de los países (método de Ward sobre los puntajes factoriales).

El dendrograma (Figura 10) muestra cómo se van uniendo los países según la semejanza de sus puntajes en los 3 componentes. Las fusiones de menor altura agrupan países casi idénticos y las de mayor altura unen grupos muy distintos. La línea punteada corta el árbol en 5 clusters, y la franja de colores bajo las ramas indica a qué grupo pertenece cada país.

Figura 11. Ancho promedio de silueta según el número de clusters.

La silueta promedio con 5 clusters es de 0,339, que según Kaufman y Rousseeuw (1990) corresponde a una estructura débil. El valor más alto de toda la Tabla 13 se da con 2 clusters (0,437), pero una partición en solo dos grupos separa únicamente los países de más y de menos infraestructura y no permite distinguir perfiles intermedios; entre las particiones de 3 o más grupos, el máximo está en 5 clusters (0,339). Se eligió trabajar con 5 clusters porque coincide con el mejor valor de silueta entre las particiones de 3 o más grupos, y produce grupos de tamaño razonable (33 / 28 / 27 / 18 / 8 países) y con una interpretación clara, como se verá a continuación. Como señalan Kaufman y Rousseeuw (1990), la silueta orienta la decisión pero no la determina sola: se combinó con el dendrograma y con la interpretabilidad de los grupos. 4,4 % de los países tiene silueta negativa, es decir, podría estar mejor asignado a otro grupo, lo que es esperable en países de frontera entre perfiles.

Caracterización de los clusters

La Tabla 14 muestra el tamaño de cada cluster y su posición promedio en los tres primeros componentes. Los clusters están numerados de menor a mayor CP1, de modo que C1 reúne a los países con menor desarrollo y conectividad y C5 a los de mayor. Recuérdese que el CP1 mide desarrollo y conectividad, el CP2 tamaño del país y el CP3 el perfil renovable y de exportaciones de alta tecnología; por eso dos clusters con un CP1 parecido pueden diferenciarse por su tamaño o por su matriz energética. Por la silueta media, el grupo más compacto y mejor separado es C1; las siluetas medias por cluster son C1 (0,44), C2 (0,23), C3 (0,38), C4 (0,29) y C5 (0,35).

La Tabla 15 compara los clusters con las medianas de nueve indicadores en sus unidades originales; se usa la mediana porque unos pocos países extremos no distorsionan el valor central del grupo. A continuación se resume cada cluster:

  • C1. Infraestructura incipiente (27 países, 24 % del total). Mediana de PIB per cápita de 4.687 USD, 29,7 % de usuarios de internet, 411 kWh por habitante, 57,6 % de electricidad renovable y una población mediana de 25,5 millones. Países más cercanos al centro del grupo: Cote d’Ivoire, Cameroon, Myanmar, Kenya y Nepal. Continentes: África (17), Asia (8) y América (2).
  • C2. Emergentes grandes y medianas (33 países, 29 % del total). Mediana de PIB per cápita de 15.960 USD, 69,7 % de usuarios de internet, 2.359 kWh por habitante, 16,2 % de electricidad renovable y una población mediana de 36,2 millones. Países más cercanos al centro del grupo: Ukraine, Peru, Thailand, Morocco y Colombia. Continentes: Asia (15), América (10), África (6) y Europa (2).
  • C3. Pequeñas de conectividad media-alta (18 países, 16 % del total). Mediana de PIB per cápita de 24.040 USD, 73,6 % de usuarios de internet, 3.730 kWh por habitante, 39,5 % de electricidad renovable y una población mediana de 2,9 millones. Países más cercanos al centro del grupo: Estonia, Costa Rica, Moldova, Georgia y Lithuania. Continentes: Europa (11), América (3), Asia (3) y África (1).
  • C4. Alto consumo energético (hidrocarburos) (8 países, 7 % del total). Mediana de PIB per cápita de 58.080 USD, 97,4 % de usuarios de internet, 12.895 kWh por habitante, 0,1 % de electricidad renovable y una población mediana de 3,5 millones. Países más cercanos al centro del grupo: Kuwait, Oman, Bahrain, Trinidad and Tobago y United Arab Emirates. Continentes: Asia (7) y América (1).
  • C5. Economías avanzadas (28 países, 25 % del total). Mediana de PIB per cápita de 51.261 USD, 89,5 % de usuarios de internet, 7.022 kWh por habitante, 27,5 % de electricidad renovable y una población mediana de 10,7 millones. Países más cercanos al centro del grupo: Sweden, Belgium, Netherlands, Ireland y Finland. Continentes: Europa (19), Asia (5), América (2) y Oceanía (2).

Figura 12. Perfil estandarizado de los clusters (media de cada variable, en desviaciones estándar).

La Figura 12 muestra, para cada cluster, cuánto se aleja el promedio de cada variable (ya transformada con logaritmo cuando correspondía) del promedio general: los valores positivos indican que el grupo está por encima del promedio de los 114 países y los negativos, por debajo. El gradiente de C1 a C5 es muy claro en las variables de ingreso, energía, internet y transporte aéreo (en C1 el PIB per cápita promedia -1,39 desviaciones estándar y en C5 1,02), mientras que las variables de tamaño (población, superficie y PIB total) y la electricidad renovable no siguen ese orden, lo que coincide con lo observado en el ACP.

Figura 13. Países en el plano CP1–CP2 según su cluster.

En el plano CP1–CP2 (Figura 13) los clusters se ordenan de izquierda a derecha según el nivel de desarrollo y conectividad. El eje vertical separa los países grandes de los pequeños dentro de cada nivel de desarrollo, por lo que los grupos se solapan en esa dirección: el tamaño del país no es lo que define la pertenencia a un cluster, sino que la distingue dentro de él. Algunos grupos parecen superponerse porque Ward trabaja con los 3 componentes a la vez y esta figura solo muestra dos de ellos.

Figura 14. Países en el plano CP1–CP3 según su cluster.

La Figura 14 añade la dimensión del perfil energético y tecnológico (CP3). Países que están cerca en el CP1 pueden quedar separados en el CP3, por ejemplo por tener una matriz eléctrica renovable o una oferta exportadora de alta tecnología, y eso explica por qué algunos grupos se distinguen entre sí aunque tengan un nivel de desarrollo parecido. Esto muestra la utilidad de clusterizar sobre varios componentes y no solo sobre el primero.

Clusters y continentes

Figura 15. Composición de cada cluster por continente (% de países).

La ubicación geográfica no se usó para formar los grupos, pero se relaciona con ellos. C1 está formado sobre todo por países de África (63 %), C2 está formado sobre todo por países de Asia (45 %), C3 está formado sobre todo por países de Europa (61 %), C4 está formado sobre todo por países de Asia (88 %) y C5 está formado sobre todo por países de Europa (68 %). Sin embargo, 4 de los 5 clusters incluyen países de tres o más continentes (Tabla 16), de modo que pertenecer al mismo continente no implica el mismo perfil de infraestructura y tecnología.

Figura 16. Distribución mundial de los clusters.

El mapa (Figura 16) permite ver si hay concentraciones regionales. Los países en gris no hacen parte de la base de análisis (103 países sin información completa). La Tabla 17 lista cada país con su cluster.

Síntesis de la clusterización

¿Cuántos clusters se formaron y qué características define a cada grupo? Ward produjo 5 clusters (silueta promedio de 0,339, una estructura débil): C1. Infraestructura incipiente, C2. Emergentes grandes y medianas, C3. Pequeñas de conectividad media-alta, C4. Alto consumo energético (hidrocarburos) y C5. Economías avanzadas. Se ordenan por nivel de desarrollo y conectividad: la mediana de usuarios de internet va de 29,7 % en C1 a 89,5 % en C5, y la del PIB per cápita de 4.687 a 51.261 USD.

¿Qué patrones aparecen entre los países? El nivel de desarrollo, y no el tamaño ni la ubicación, es lo que más separa a los grupos. El continente se asocia con el cluster (C1 es sobre todo de África y C5 de Europa), pero 4 de los 5 clusters mezclan tres o más continentes, y dentro de un mismo continente hay países de perfiles distintos.

Conclusiones

El ACP resumió las 17 variables de 114 países en 3 componentes que acumulan 74,4 % de la variabilidad: un eje de desarrollo y conectividad (CP1, 52,6 %), un eje de tamaño del país (CP2, 14,2 %) y un eje de perfil renovable y de exportaciones de alta tecnología (CP3, 7,6 %). Sobre esos puntajes, el método de Ward formó 5 clusters con una silueta promedio de 0,339, que ordenan a los países por nivel de desarrollo, desde infraestructura incipiente (C1) hasta economías avanzadas (C5).

Las hipótesis de trabajo planteadas en la introducción se contrastan así:

  • H1 se confirma: 3 componentes concentran 74,4 % de la variabilidad de 17 indicadores.
  • H2 se confirma: el CP1 recoge el nivel de desarrollo y conectividad, y el CP2, de forma casi independiente, el tamaño de los países.
  • H3 se confirma: los clusters se ordenan por desarrollo y 4 de los 5 incluyen países de tres o más continentes, de modo que la similitud en infraestructura y tecnología depende más del perfil de la economía que de su ubicación geográfica.

El método de Ward resultó útil para organizar la información y para identificar perfiles de países que facilitan las comparaciones internacionales.

Limitaciones

  • Cobertura. El ACP exige datos completos, así que se trabajó con 114 de 217 países; los 103 excluidos son sobre todo economías pequeñas, y los resultados no describen a todo el universo mundial.
  • Un solo año. Los resultados reflejan 2019 y no necesariamente la situación actual.
  • Sensibilidad. Los grupos dependen de las variables elegidas y del número de clusters. La silueta con 5 clusters es de 0,339 (una estructura débil), y las fusiones de Ward son irreversibles.
  • CP3. Depende mucho de pocos países (Qatar aporta 15,7 % de su inercia), por lo que su lectura debe tomarse con cautela.
  • Alcance. El análisis es descriptivo: identifica semejanzas y diferencias entre países, pero no establece relaciones de causa y efecto.

Recomendaciones

  • Este tipo de análisis puede usarse para comparar países e identificar necesidades de infraestructura, acceso a tecnologías, energía y desarrollo económico, y para definir grupos de referencia (benchmarking).
  • En ingeniería industrial, los resultados pueden apoyar el análisis de diferencias en infraestructura y recursos entre países, y orientar estudios de inversión y planificación.
  • Conviene repetir el análisis con otros años para ver cómo cambian los indicadores y si algunos países pasan de un cluster a otro.
  • Conviene complementar los resultados con variables sociales y ambientales para tener una visión más completa del desarrollo de los países.

Bibliografía

Abdi, H., & Williams, L. J. (2010). Principal component analysis. WIREs Computational Statistics, 2(4), 433–459. https://doi.org/10.1002/wics.101

Banco Mundial. (s. f.). World Development Indicators. https://databank.worldbank.org/source/world-development-indicators

Bartlett, M. S. (1954). A further note on the multiplying factors for various χ² approximations in factor analysis. Journal of the Royal Statistical Society: Series B, 16(2), 296–298. https://doi.org/10.1111/j.2517-6161.1954.tb00174.x

Cattell, R. B. (1966). The scree test for the number of factors. Multivariate Behavioral Research, 1(2), 245–276. https://doi.org/10.1207/s15327906mbr0102_10

Hair, J. F., Black, W. C., Babin, B. J., & Anderson, R. E. (2019). Multivariate data analysis (8.ª ed.). Cengage.

Husson, F., Lê, S., & Pagès, J. (2017). Exploratory multivariate analysis by example using R (2.ª ed.). Chapman & Hall/CRC.

James, G., Witten, D., Hastie, T., & Tibshirani, R. (2021). An introduction to statistical learning: With applications in R (2.ª ed.). Springer.

Jolliffe, I. T., & Cadima, J. (2016). Principal component analysis: A review and recent developments. Philosophical Transactions of the Royal Society A, 374(2065), 20150202. https://doi.org/10.1098/rsta.2015.0202

Kaiser, H. F. (1960). The application of electronic computers to factor analysis. Educational and Psychological Measurement, 20(1), 141–151. https://doi.org/10.1177/001316446002000116

Kaiser, H. F., & Rice, J. (1974). Little Jiffy, Mark IV. Educational and Psychological Measurement, 34(1), 111–117. https://doi.org/10.1177/001316447403400115

Kaufman, L., & Rousseeuw, P. J. (1990). Finding groups in data: An introduction to cluster analysis. Wiley.

Murtagh, F., & Legendre, P. (2014). Ward’s hierarchical agglomerative clustering method: Which algorithms implement Ward’s criterion? Journal of Classification, 31(3), 274–295. https://doi.org/10.1007/s00357-014-9161-z

Naciones Unidas. (2015). Transformar nuestro mundo: la Agenda 2030 para el Desarrollo Sostenible (A/RES/70/1).

Rousseeuw, P. J. (1987). Silhouettes: A graphical aid to the interpretation and validation of cluster analysis. Journal of Computational and Applied Mathematics, 20, 53–65. https://doi.org/10.1016/0377-0427(87)90125-7

Ward, J. H. (1963). Hierarchical grouping to optimize an objective function. Journal of the American Statistical Association, 58(301), 236–244. https://doi.org/10.1080/01621459.1963.10500845