Horizontes Amazónicos e-ISSN: 3121-2581 Vol. 2, Núm. 2 (2026) pp. 1–13 DOI: 10.67321/gxffd102 horizontesamazonicos.org

Artículo de investigación

¿Chocolate o riqueza? Replicación transnacional de la correlación entre consumo de chocolate y Premios Nobel con datos actualizados (1901–2025), inclusión de América Latina y modelos de conteo con exceso de ceros

Chocolate or wealth? A cross-national replication of the chocolate consumption–Nobel Prize correlation with updated data (1901–2025), inclusion of Latin America and zero-inflated count models

Samantha Lisseth Romero Ruiz1[ORCID]
1 Universidad Estatal Amazónica, Puyo, Pastaza, Ecuador
Recibido: 2026-08-04Aceptado: 2026-09-20Publicado: 2026-10-13

Resumen

La correlación entre consumo per cápita de chocolate y laureados Nobel por 10 millones de habitantes (r = 0.791) es un ejemplo clásico de falacia ecológica, rara vez replicado con datos armonizados y modelos para conteos con exceso de ceros. Se buscó replicarla y ampliarla incorporando a América Latina y el Caribe (ALC) y controles macroeconómicos históricos. Se realizó un estudio ecológico, retrospectivo y transnacional con 129 países (19 de ALC) y 991 premios individuales (1901–2025), atribuidos por país de nacimiento y de afiliación. La exposición fue el cacao disponible para consumo (FAOSTAT); los controles, el PIB per cápita de 1960–2022, el gasto en I+D, la matrícula terciaria y las publicaciones científicas. El análisis avanzó de correlaciones simples y parciales a modelos de Poisson, binomial negativa (BN), inflados por ceros (ZIP, ZINB) y hurdle. En los 23 países originales, la correlación no se reprodujo con la medida armonizada (r = 0.19; p = 0.37). En la muestra ampliada, la correlación cruda (r = 0.37; p < 0.001) desapareció al ajustar por PIB (r parcial = 0.10; p = 0.24). En los modelos BN y ZINB ajustados, la razón de tasas de incidencia del chocolate osciló entre 1.09 y 1.15, con intervalos que incluyeron la unidad; el PIB y la I+D fueron predictores robustos. La leche, usada como control de falsificación, mostró una asociación residual mayor que la del chocolate. La riqueza y la capacidad científica, y no el chocolate, explican la distribución de los Premios Nobel

Abstract

The correlation between per capita chocolate consumption and Nobel laureates per 10 million inhabitants (r = 0.791) became a textbook example of the ecological fallacy, yet it has rarely been replicated with harmonized data, unselected samples and models suited to zero-heavy counts. We aimed to replicate and extend that finding by including Latin America and the Caribbean (LAC) and historically lagged macroeconomic controls. We conducted an ecological, retrospective, cross-national study of 129 countries (19 in LAC) and 991 individual prizes awarded between 1901 and 2025, attributed by country of birth and of affiliation. Exposure was cocoa available for consumption (FAOSTAT); controls were mean GDP per capita for 1960–2022, R&D expenditure, tertiary enrollment and scientific output. The analysis moved from simple and partial correlations to Poisson, negative binomial (NB), zero-inflated Poisson and NB (ZIP, ZINB) and hurdle models. In Messerli's 23 countries the correlation was not reproduced with the harmonized measure (r = 0.19; p = 0.37). In the extended sample, the crude correlation (r = 0.37; p < 0.001) vanished after adjustment for GDP (partial r = 0.10; p = 0.24). In adjusted NB and ZINB models, the incidence rate ratio for chocolate ranged from 1.09 to 1.15 with intervals including one, whereas GDP and R&D were robust predictors. Milk consumption, used as a falsification control, showed a stronger residual association than chocolate. Wealth and scientific capacity, not chocolate, explain the distribution of Nobel Prizes.

Palabras clave

falacia ecológica; replicación; regresión binomial negativa inflada por ceros; Premios Nobel; América Latina

Keywords

ecological fallacy; replication; zero-inflated negative binomial regression; Nobel Prize; Latin America


1. Introducción

En 2012, Messerli publicó en el New England Journal of Medicine una nota breve que vinculaba el consumo anual de chocolate per cápita de 23 países con su número de laureados Nobel por cada 10 millones de habitantes. La correlación lineal fue notablemente alta (r = 0.791; p < 0.0001) y el autor propuso, de forma explícitamente hipotética, que los flavanoles del cacao podrían mejorar la función cognitiva de la población y, con ello, la probabilidad de producir científicos galardonados (Messerli, 2012). El ejercicio, posiblemente concebido como una provocación, se difundió en la prensa de forma literal y despojado de su ironía, hasta convertirse en el ejemplo contemporáneo más citado de los riesgos de interpretar causalmente una correlación ecológica (McClintock et al., 2014; Szromek, 2026).

La literatura posterior se organizó en varias líneas. La primera fueron las críticas y replicaciones directas: Maurage et al. (2013) mostraron que la interpretación basada en flavanoles se debilita cuando se consideran otros alimentos ricos en estos compuestos, y que correlaciones igual de altas pueden obtenerse con variables sin plausibilidad biológica; Linthwaite y Fuller (2013) reprodujeron un patrón análogo con el consumo de leche, y Dunstan (2013) subrayó que el diseño era ecológico, generador de hipótesis y temporalmente desalineado, porque el consumo actual no puede explicar descubrimientos realizados décadas atrás. La segunda línea examinó los determinantes económicos: Li (2013) señaló al ingreso nacional como causa común plausible del consumo de chocolate y de los premios, y Doi et al. (2014), con 220 países y partición jerárquica de la varianza, encontraron que el número de publicaciones y el gasto en investigación predicen los premios Nobel mejor que los hábitos dietéticos y que los factores económicos generales. Esta conclusión coincide con la literatura cienciométrica, según la cual la riqueza científica de los países depende de la actividad tecnológica, la política de investigación y la calidad institucional más que del ingreso por sí solo (Pietrucha, 2018; Rodríguez-Navarro & Brito, 2022).

No obstante, la evidencia no es unánime. Prinz (2020), mediante un modelo de selección de Heckman en dos etapas sobre 70 países, informó que una asociación positiva entre chocolate y laureados persistía tras ajustar por PIB per cápita, gasto en I+D, publicaciones y otros consumos, aunque interpretó ese residuo como reflejo de heterogeneidad no observada y no como prueba causal. Así, el desacuerdo en la literatura es estrecho: se discute si la correlación agregada desaparece por completo tras el ajuste, pero existe consenso en que no justifica una inferencia causal (Doi et al., 2014; Prinz, 2020). En paralelo, la evidencia clínica sobre el mecanismo propuesto es, en el mejor de los casos, modesta: algunos ensayos pequeños mostraron mejoras en tareas ejecutivas o de memoria tras 8–12 semanas de cacao rico en flavanoles (Desideri et al., 2012; Mastroiacovo et al., 2015; Brickman et al., 2014), pero el ensayo COSMOS-Mind, con 2262 adultos mayores que recibieron 500 mg/día de flavanoles durante tres años, no encontró efecto sobre la cognición global (Baker et al., 2023), resultado confirmado en la subcohorte clínica del mismo ensayo (Vyas et al., 2024) y matizado por metaanálisis que solo identifican beneficios pequeños y específicos de dominio (Cheng et al., 2022; Arabpour et al., 2026).

Persisten, sin embargo, tres vacíos metodológicos que justifican una replicación actualizada. Primero, el sesgo de selección: el análisis original incluyó sobre todo países ricos de Europa Occidental y América del Norte que ya contaban con laureados, truncando la varianza de las variables y excluyendo a la mayoría de naciones con cero premios (Li, 2013; Prinz, 2020). América Latina y el Caribe, región que reúne países productores de cacao y sistemas científicos emergentes, aparece en la literatura solo como contexto general y no como objeto de replicación directa (Rizzuto et al., 2017; Cortés & Andrade, 2022). Segundo, la atribución nacional de los laureados: contar por país de nacimiento o por afiliación institucional responde a preguntas distintas, porque los galardonados son altamente móviles y la mayoría realiza su trabajo premiado en unos pocos países, en especial Estados Unidos (Schlagberger et al., 2016; Heinze et al., 2020; Heinze & Fuchs, 2022). Tercero, la medición y el modelado: los datos de consumo de chocolate provienen de fuentes industriales heterogéneas y no armonizadas (Apgar & Tarka, 2019; Del Prete & Samoggia, 2020), la variable dependiente es un conteo con exceso de ceros que viola los supuestos de la correlación de Pearson, y la latencia entre las condiciones nacionales y el reconocimiento Nobel exige promedios históricos de las covariables (Dunstan, 2013; Li, 2025). Más allá de su carácter anecdótico, el caso chocolate–Nobel funciona como un banco de pruebas metodológico para la cienciometría y la epidemiología nutricional: permite cuantificar cuánto de una correlación ecológica sobrevive cuando se corrigen la selección de la muestra, la medición de la exposición, la estructura de confusión y la distribución de la variable respuesta (Loney & Nagelkerke, 2014; Pollet et al., 2014). Para una región como América Latina, donde varios países son grandes productores de cacao fino de aroma y la mayoría invierte menos del 1% de su PIB en I+D, la pregunta tiene además una lectura de política científica: si el chocolate no explica los premios, ¿qué los explica y por qué la región aparece rezagada?

Por lo tanto, el objetivo de este estudio es replicar la correlación entre consumo de chocolate y laureados Nobel con datos actualizados a 2025 y ampliarla a una muestra transnacional que incluye estratégicamente a América Latina y el Caribe, evaluando si la asociación persiste tras controlar por el PIB per cápita histórico, el gasto en I+D, la matrícula terciaria y la producción científica, mediante un plan escalonado que va desde la correlación simple hasta modelos de conteo con exceso de ceros, y con desagregación de los laureados por país de nacimiento y de afiliación.

2. Método

2.1. Diseño del estudio

Se realizó un estudio observacional, ecológico, retrospectivo y transnacional, en el que la unidad de análisis fue el país. El diseño replicó la especificación original de Messerli (2012) —laureados por cada 10 millones de habitantes frente a consumo per cápita de chocolate— y la amplió en cuatro dimensiones: (a) una muestra no restringida a países con laureados; (b) una exposición armonizada internacionalmente; (c) controles macroeconómicos y científicos medidos como promedios históricos, y (d) modelos de conteo que reconocen la sobredispersión y el exceso de ceros. El marco causal se formalizó mediante un grafo acíclico dirigido (Fig. 1), en el que la riqueza histórica y la geografía cultural actúan como causas comunes del consumo aparente de chocolate y de la producción de laureados, mientras que la I+D, la educación terciaria y la migración científica son mediadores o determinantes directos de los premios.

Fig. 1 Grafo acíclico dirigido que representa la estructura causal supuesta. La riqueza histórica y la geografía cultural actúan como causas comunes del consumo aparente de chocolate y de la producción de laureados; la I+D, la educación terciaria, las publicaciones y la migración científica operan como determinantes o mediadores. Recuadros azules: confusores; naranja: exposición; gris: mediadores y desenlace. La flecha discontinua representa el efecto causal puesto a prueba.
Fig. 1. Grafo acíclico dirigido que representa la estructura causal supuesta. La riqueza histórica y la geografía cultural actúan como causas comunes del consumo aparente de chocolate y de la producción de laureados; la I+D, la educación terciaria, las publicaciones y la migración científica operan como determinantes o mediadores. Recuadros azules: confusores; naranja: exposición; gris: mediadores y desenlace. La flecha discontinua representa el efecto causal puesto a prueba.
Fuente: elaboración propia.

2.2. Población, muestra y entorno

La población de referencia fueron los Estados con datos de población en 2023 (n = 159 países con al menos un millón de habitantes). Se excluyeron los microestados (< 1 millón) porque un único laureado genera tasas extremas e inestables (por ejemplo, Santa Lucía con dos laureados); su inclusión se evaluó en el análisis de sensibilidad. La muestra analítica de casos completos estuvo formada por 129 países con información disponible para la exposición y las cuatro covariables principales, de los cuales 19 pertenecen a América Latina y el Caribe (Argentina, Brasil, Chile, Colombia, Costa Rica, Cuba, Ecuador, El Salvador, Guatemala, Honduras, Jamaica, México, Nicaragua, Panamá, Paraguay, Perú, Trinidad y Tobago, Uruguay y Venezuela). Para la replicación directa se conservó además el subconjunto de los 23 países analizados por Messerli (2012). La recolección y el análisis de datos se efectuaron entre septiembre y octubre de 2026 a partir de fuentes secundarias de acceso abierto; el periodo de observación de los premios abarcó desde 1901 hasta 2025.

2.3. Variables y recolección de datos

Variable dependiente. Se construyó una base de 991 premios individuales (1901–2025) en Física (230), Química (200), Fisiología o Medicina (232), Literatura (122), Paz (108) y Ciencias Económicas (99); se excluyeron las organizaciones galardonadas. Los registros de 1901–2016 se obtuvieron de la base pública derivada de la API de la Fundación Nobel (rfordatascience, 2019) y los de 2017–2025 se completaron manualmente a partir de nobelprize.org (Nobel Prize Outreach, 2025). Cada laureado se atribuyó a su país de nacimiento según las fronteras actuales y, para las ciencias, al país de la primera institución declarada en el momento del premio. Se definieron tres conteos: Y1, todos los premios por país de nacimiento (especificación de Messerli); Y2, premios de ciencias (Física, Química y Medicina) por nacimiento, y Y3, premios de ciencias por afiliación. La tasa de laureados se calculó como:

Ti=NiPi×107(1)

donde Nᵢ es el número de laureados atribuidos al país i y Pᵢ su población en 2023.

Exposición principal. Se utilizó la cantidad de cacao en grano y productos derivados disponible para consumo humano (kg/hab/año) de las Hojas de Balance de Alimentos de FAOSTAT, difundida por Our World in Data (FAO, 2025; Our World in Data, 2025a), promediando los años 2015 y 2022. A diferencia de las cifras industriales empleadas originalmente, esta medida es armonizada, pública y cubre a casi todos los países, aunque representa disponibilidad aparente y no ingesta individual.

Controles. El PIB per cápita real (dólares internacionales de 2011) se promedió para 1960–2022 a partir del Maddison Project Database 2023 (Bolt & van Zanden, 2024), con el fin de respetar la latencia entre las condiciones nacionales y el reconocimiento de los descubrimientos. El gasto interno bruto en I+D (% del PIB) se promedió con los cortes 2000, 2010, 2015 y 2021, y la tasa bruta de matrícula terciaria con los cortes 1990, 2000, 2010 y 2020, ambos de UNESCO-UIS y del Banco Mundial (Our World in Data, 2025b, 2025c). La producción científica (artículos científicos y técnicos por millón de habitantes, promedio de 2000 y 2020; Banco Mundial, 2025) se empleó como covariable adicional y como predictor del componente de ceros estructurales. Como control de falsificación se incluyó el suministro de leche (kg/hab/año, promedio 1961–2017; FAO vía Our World in Data), variable sin plausibilidad neurocognitiva previamente asociada con los premios (Linthwaite & Fuller, 2013). La Tabla 1 resume la operacionalización.

Tabla 1. Operacionalización de las variables, fuentes y ventanas temporales utilizadas
VariableDefinición operativaFuenteVentanaUso en el análisis
Y1: laureados, todos los premiosLaureados individuales nacidos en el territorio actual del país (Física, Química, Medicina, Literatura, Paz, Economía); se excluyen organizacionesNobel Prize Outreach (2025); rfordatascience (2019)1901–2025Conteo dependiente (réplica de Messerli)
Y2: laureados en ciencias (nacimiento)Laureados en Física, Química y Fisiología o Medicina por país de nacimientoÍdem1901–2025Conteo dependiente
Y3: laureados en ciencias (afiliación)País de la primera institución declarada al recibir el premioÍdem1901–2025Conteo dependiente
Tasa por 10 M hab.Nᵢ / Pᵢ × 10⁷; Pᵢ = población 2023Our World in Data (población)2023Correlaciones; offset ln(Pᵢ/10⁷)
Chocolate (exposición)Cacao en grano y productos disponible para consumo humano (kg/hab/año), media de dos cortesFAOSTAT, Hojas de Balance (FAO, 2025)2015 y 2022Exposición principal
PIB per cápitaPIB real per cápita (USD int. 2011), media anual; en logaritmo naturalMaddison Project 2023 (Bolt & van Zanden, 2024)1960–2022Confusor principal
I+DGasto interno bruto en I+D (% del PIB), media de cortes disponiblesUNESCO-UIS / Banco Mundial2000, 2010, 2015, 2021Confusor
Matrícula terciariaTasa bruta de matrícula en educación terciaria (%), por 10 puntos porcentualesUNESCO-UIS1990, 2000, 2010, 2020Confusor
PublicacionesArtículos científicos y técnicos por millón de habitantes; en lnBanco Mundial (IP.JRN.ARTC.SC)2000 y 2020Covariable (M5) e inflación de ceros (ZINB)
LecheSuministro de leche, excluida mantequilla (kg/hab/año)FAOSTAT vía Our World in Data1961–2017Control de falsificación
ALCIndicador de América Latina y el Caribe (19 países)Clasificación regional—Modelo M4; análisis regional

Fuente: elaboración propia.

2.4. Análisis estadístico

El análisis siguió un plan escalonado de cinco etapas. Etapa 1 (replicación directa): en los 23 países de Messerli se estimaron los coeficientes de Pearson y Spearman entre la exposición armonizada y la tasa de laureados, tanto con los premios hasta 2011 (horizonte original) como hasta 2025. Etapa 2 (ampliación bivariada y parcial): en la muestra de 129 países se estimaron correlaciones simples y correlaciones parciales de Pearson y de Spearman (rangos residualizados) controlando por el logaritmo del PIB per cápita y, sucesivamente, por I+D, matrícula terciaria y publicaciones.

Etapa 3 (modelos de conteo): dado que la variable respuesta es un conteo con alta proporción de ceros (Cameron & Trivedi, 2013), se ajustaron modelos lineales generalizados con enlace logarítmico y desplazamiento (offset) poblacional, de modo que los coeficientes exponenciados se interpretan como razones de tasas de incidencia (RTI) por cada 10 millones de habitantes:

ln⁡μi=ln⁡(Pi107)+xi⊤β(2)

donde xᵢ contiene el consumo de cacao (kg/hab/año), el logaritmo natural del PIB per cápita, la I+D (% del PIB) y la matrícula terciaria (por 10 puntos porcentuales). Se ajustaron cinco familias: Poisson; binomial negativa tipo 2 (BN), con varianza

Var(Yi)=μi+αμi2(3)

Poisson inflada por ceros (ZIP; Lambert, 1992); binomial negativa inflada por ceros (ZINB), que mezcla un proceso de ceros estructurales con probabilidad πᵢ y un proceso BN:

P(Yi=0)=πi+(1−πi)(1+αμi)−1/α(4)
P(Yi=k)=(1−πi)fBN(k;μi,α),k≥1(5)
logit(πi)=γ0+γ1ln⁡Arti(6)

y un modelo hurdle de dos partes (Mullahy, 1986), con una regresión logística para la probabilidad de tener al menos un laureado y una BN truncada en cero para los conteos positivos, análogo conceptual del modelo de selección de Heckman (1979) empleado por Prinz (2020). El componente de inflación se especificó con el logaritmo de las publicaciones por millón, siguiendo el hallazgo de Doi et al. (2014) de que la actividad científica discrimina la entrada de un país en el mapa Nobel. Para cada familia se estimaron cinco especificaciones anidadas: M1 (solo chocolate), M2 (+ PIB), M3 (+ I+D y matrícula terciaria; modelo principal), M4 (+ indicador ALC) y M5 (+ publicaciones).

La sobredispersión se evaluó con el estadístico χ² de Pearson dividido por los grados de libertad del modelo de Poisson y con la prueba de significación de α en la BN. La selección entre modelos se basó en el criterio de información de Akaike (AIC), el criterio bayesiano (BIC) y la prueba de Vuong (1989) para modelos no anidados (ZINB frente a BN). Etapa 4 (sensibilidad): el modelo M3 se reestimó excluyendo productores de cacao y centros de procesamiento/reexportación (cuyo dato FAOSTAT refleja el comercio y no el consumo), excluyendo a Suecia y Suiza (puntos influyentes en el análisis original), excluyendo a Estados Unidos y Reino Unido (principales receptores de talento), usando el PIB medio de 1990–2022 e incluyendo países de cualquier tamaño. Etapa 5 (falsificación y contexto regional): se sustituyó el chocolate por la leche y se añadió un indicador de Europa para evaluar confusión geocultural residual; finalmente, se describieron los laureados nacidos en ALC y su afiliación al recibir el premio. Se adoptó un nivel de significación de 0.05 (bilateral) y se informan intervalos de confianza del 95% (IC 95%). Los análisis se realizaron en Python 3.13 con statsmodels 0.15 (Seabold & Perktold, 2010) y SciPy.

3. Resultados

3.1. Características de la muestra

De los 991 premios individuales, 972 correspondieron a personas nacidas en alguno de los 129 países de la muestra analítica. La distribución fue marcadamente asimétrica y con exceso de ceros: el 45.7% de los países no tenía ningún laureado por nacimiento, el 60.5% ninguno en ciencias por nacimiento y el 79.8% ninguno en ciencias por afiliación; en ALC estas proporciones fueron 47.4%, 78.9% y 94.7%, respectivamente. Los 19 países latinoamericanos sumaron 19 laureados por nacimiento (5 en ciencias) y solo 2 premios de ciencias por afiliación, frente a 655 en el resto del mundo. ALC no se diferenció del resto en el consumo aparente de cacao, el PIB per cápita histórico ni la matrícula terciaria, pero sí en el gasto en I+D, con una mediana equivalente a la mitad (0.24% frente a 0.48% del PIB; p = 0.003), y mostró una tendencia a menor producción científica per cápita (p = 0.071) (Tabla 2).

Tabla 2. Características de los países de la muestra analítica según región (mediana [rango intercuartílico])
VariableALC (n = 19)Resto (n = 110)p
Cacao disponible (kg/hab/año)0.96 [0.53–1.13]0.84 [0.14–2.14]0.850
PIB pc 1960–2022 (USD int. 2011)8,233 [5,665–11,403]8,404 [2,774–18,219]0.839
I+D (% PIB)0.24 [0.08–0.37]0.48 [0.20–1.04]0.003
Matrícula terciaria (%)40.92 [19.88–44.99]29.73 [8.45–61.54]0.683
Artículos por millón hab.68 [23–109]120 [17–671]0.071
Laureados, todos (nacimiento)1.00 [0.00–2.00]; Σ = 191.00 [0.00–4.75]; Σ = 9530.360
Laureados, ciencias (nacimiento)0.00 [0.00–0.00]; Σ = 50.00 [0.00–2.00]; Σ = 6490.035
Laureados, ciencias (afiliación)0.00 [0.00–0.00]; Σ = 20.00 [0.00–0.00]; Σ = 6550.076
Tasa, todos (por 10 M)0.05 [0.00–0.79]0.08 [0.00–2.53]0.396

Nota: p de la prueba U de Mann-Whitney. Σ = total de laureados. ALC = América Latina y el Caribe.

Fuente: elaboración propia con datos de FAOSTAT, Maddison Project, UNESCO-UIS, Banco Mundial y Nobel Prize Outreach.

3.2. Replicación directa y correlaciones ampliadas

Con la medida armonizada de FAOSTAT, la correlación en los 23 países de Messerli no se reprodujo: r = 0.18 (p = 0.42) con los premios hasta 2011 y r = 0.19 (p = 0.37) con los premios hasta 2025 (Fig. 2A). En contraste, en esos mismos países el logaritmo del PIB per cápita histórico (r = 0.59; p = 0.003), el gasto en I+D (r = 0.56; p = 0.005) y el consumo de leche (r = 0.63; p = 0.001) se correlacionaron significativamente con la tasa de laureados. Algunos de los países con mayor tasa, como Suiza y los Países Bajos, mostraron valores bajos de cacao disponible en FAOSTAT, coherentes con su papel de procesadores y reexportadores.

En la muestra ampliada de 129 países la correlación cruda entre chocolate y laureados fue positiva y significativa para los tres desenlaces (r = 0.37, 0.34 y 0.24; p ≤ 0.006), pero de menor magnitud que la de cualquiera de los indicadores de riqueza o capacidad científica. Tras ajustar únicamente por el PIB per cápita, la correlación parcial se redujo a 0.10 (p = 0.24), 0.07 (p = 0.42) y 0.01 (p = 0.91), respectivamente, y se anuló al añadir I+D y matrícula terciaria (Tabla 3; Fig. 2B). El consumo de cacao se correlacionó a su vez con el PIB (r = 0.57), la matrícula terciaria (r = 0.57), las publicaciones (r = 0.55) y la I+D (r = 0.40), lo que confirma su papel de marcador de prosperidad.

Tabla 3. Replicación directa y ampliada: correlaciones simples y parciales con la tasa de laureados Nobel por 10 millones de habitantes
MuestraPredictor / ajusteTodos los premios (nacimiento) r (p)Ciencias (nacimiento) r (p)Ciencias (afiliación) r (p)
Messerli, 23 paísesChocolate, Pearson (premios hasta 2011)0.18 (0.422)——
Chocolate, Pearson0.19 (0.373)0.07 (0.763)0.08 (0.713)
Chocolate, Spearman0.18 (0.404)0.08 (0.720)0.19 (0.385)
ln PIB per cápita, Pearson0.59 (0.003)0.56 (0.005)0.51 (0.012)
I+D, Pearson0.56 (0.005)0.61 (0.002)0.61 (0.002)
Leche, Pearson0.63 (0.001)0.56 (0.005)0.47 (0.024)
Ampliada, 129 paísesChocolate, Pearson0.37 (< 0.001)0.34 (< 0.001)0.24 (0.006)
ln PIB per cápita, Pearson0.53 (< 0.001)0.51 (< 0.001)0.41 (< 0.001)
I+D, Pearson0.69 (< 0.001)0.71 (< 0.001)0.65 (< 0.001)
Matrícula terciaria, Pearson0.54 (< 0.001)0.53 (< 0.001)0.36 (< 0.001)
ln publicaciones, Pearson0.58 (< 0.001)0.57 (< 0.001)0.43 (< 0.001)
Leche, Pearson (n = 121)0.68 (< 0.001)0.66 (< 0.001)0.53 (< 0.001)
Chocolate | ln PIB (parcial, Pearson)0.10 (0.239)0.07 (0.425)0.01 (0.914)
Chocolate | ln PIB (parcial, Spearman)-0.02 (0.807)-0.04 (0.640)-0.12 (0.182)
Chocolate | PIB, I+D, terciaria (parcial)0.07 (0.457)0.02 (0.804)-0.00 (0.981)
Chocolate | PIB, I+D, terciaria, publicaciones (parcial)0.07 (0.460)0.02 (0.807)-0.00 (0.986)

Nota: la exposición es el cacao disponible para consumo de FAOSTAT (2015–2022). Para comparación, Messerli (2012) informó r = 0.791 (p < 0.0001) con cifras industriales de consumo de chocolate.

Fuente: elaboración propia.

Fig. 2 Relación entre el cacao disponible para consumo (FAOSTAT, media 2015–2022) y los laureados Nobel por 10 millones de habitantes (país de nacimiento, 1901–2025). (A) Replicación en los 23 países de Messerli (2012); la línea es el ajuste por mínimos cuadrados con su banda de confianza del 95%. (B) Muestra ampliada de 129 países con América Latina y el Caribe resaltada; eje vertical en escala log(1 + y) para conservar los países sin laureados. Los recuadros muestran las correlaciones de Pearson simples y parciales.
Fig. 2. Relación entre el cacao disponible para consumo (FAOSTAT, media 2015–2022) y los laureados Nobel por 10 millones de habitantes (país de nacimiento, 1901–2025). (A) Replicación en los 23 países de Messerli (2012); la línea es el ajuste por mínimos cuadrados con su banda de confianza del 95%. (B) Muestra ampliada de 129 países con América Latina y el Caribe resaltada; eje vertical en escala log(1 + y) para conservar los países sin laureados. Los recuadros muestran las correlaciones de Pearson simples y parciales.
Fuente: elaboración propia.

3.3. Modelos de conteo con exceso de ceros

Los datos mostraron una sobredispersión importante (χ² de Pearson/gl del modelo de Poisson M3 = 6.25, 4.66 y 2.72) y el parámetro α de la BN fue significativamente mayor que cero en los tres desenlaces (α = 0.78–0.99; p ≤ 0.004), por lo que los modelos de Poisson y ZIP, que sí asignaron significación estadística al chocolate, resultaron inadecuados (AIC entre 1.5 y 1.9 veces superior al de la BN). Para todos los premios por nacimiento, el componente de inflación de la ZINB convergió hacia cero y el modelo colapsó en la BN (Vuong z ≈ 0; p = 0.50); el modelo hurdle obtuvo el menor AIC, aunque la BN tuvo el menor BIC. Para los premios de ciencias, la ZINB mejoró el ajuste por nacimiento (Vuong z = 2.65; p = 0.004), aunque su componente de inflación mostró separación casi completa, y, por afiliación, obtuvo el menor AIC (207.2), con un componente de ceros estructurales explicado por la escasa producción científica (γ₁ = −1.86; p = 0.003) (Tabla 4).

Tabla 4. Comparación de modelos de conteo para el modelo principal M3 (chocolate + ln PIB + I+D + matrícula terciaria)
DesenlaceModeloLog-verosimilitudkAICBICRTI chocolate (IC 95%)p
Todos (nacimiento)Poisson-408.35826.6840.91.23 (1.14–1.32)< 0.001
Binomial negativa-225.36462.6479.81.11 (0.89–1.38)0.347
ZIP-366.87747.5767.51.15 (1.07–1.23)< 0.001
ZINB-225.38466.6489.51.11 (0.89–1.38)0.347
Hurdle (logit + BN truncada)-214.612453.2 *487.51.02 (0.78–1.32)0.908
Diagnósticoχ²/gl Poisson = 6.25α (BN) = 0.78Vuong ZINB vs BN: z = 0.000.501
Ciencias (nacimiento)Poisson-292.15594.2608.51.16 (1.06–1.26)0.001
Binomial negativa-170.66353.2370.31.09 (0.84–1.41)0.520
ZIP-261.27536.4556.41.08 (0.99–1.18)0.075
ZINB-167.18350.1373.01.11 (0.85–1.45)0.448
Hurdle (logit + BN truncada)-163.012350.0 *384.30.98 (0.72–1.32)0.882
Diagnósticoχ²/gl Poisson = 4.66α (BN) = 0.94Vuong ZINB vs BN: z = 2.650.004
Ciencias (afiliación)Poisson-200.65411.1425.41.27 (1.16–1.38)< 0.001
Binomial negativa-100.76213.4230.51.15 (0.85–1.55)0.379
ZIP-156.47326.8346.91.15 (1.06–1.26)< 0.001
ZINB-95.68207.2 *230.01.09 (0.86–1.37)0.465
Hurdle (logit + BN truncada)-94.012211.9246.21.06 (0.83–1.36)0.640
Diagnósticoχ²/gl Poisson = 2.72α (BN) = 0.98Vuong ZINB vs BN: z = 1.280.101

Nota: RTI = razón de tasas de incidencia por 1 kg/hab/año de cacao; en el modelo hurdle corresponde a la parte de conteo. * menor AIC. (—) IC no estimable por inestabilidad de la matriz de información. Todos los modelos incluyen offset ln(población/10⁷); la inflación de ceros se modeló con ln publicaciones por millón.

Fuente: elaboración propia.

En todas las familias con varianza adecuada, la asociación del chocolate se atenuó drásticamente al introducir el PIB. En el modelo principal M3, la RTI por kg/hab/año fue de 1.11 (IC 95%: 0.89–1.38) en la BN para todos los premios, 1.11 (0.85– 1.45) en la ZINB para ciencias por nacimiento y 1.09 (0.86–1.37) en la ZINB para ciencias por afiliación, sin significación estadística (Tabla 5; Fig. 3). Por el contrario, el PIB per cápita histórico fue el predictor más fuerte en todos los desenlaces (RTI por unidad logarítmica entre 2.63 y 30.27; p ≤ 0.005) y la I+D se asoció de forma independiente con más laureados (RTI 1.57–2.83 por punto porcentual del PIB; p ≤ 0.042). Al añadir las publicaciones (M5), estas fueron también un predictor robusto (RTI 1.92, 3.67 y 5.25; p ≤ 0.004) y la RTI del chocolate se aproximó aún más a la unidad (1.02–1.10). En la parte logística del modelo hurdle, el chocolate no predijo la probabilidad de tener algún laureado por nacimiento (OR = 1.11; IC 95%: 0.68– 1.81), aunque mostró una tendencia no significativa para la afiliación en ciencias (OR = 2.07; 0.95–4.51; p = 0.067), en un componente con solo 26 países positivos.

Tabla 5. Razones de tasas de incidencia (IC 95%) del modelo principal M3 según desenlace y familia de modelo
PredictorTodos (nac.) BNTodos (nac.) ZINBCiencias (nac.) BNCiencias (nac.) ZINBCiencias (afil.) BNCiencias (afil.) ZINB
Chocolate (kg/hab/año)1.11 (0.89–1.38)1.11 (0.89–1.38)1.09 (0.84–1.41)1.11 (0.85–1.45)1.15 (0.85–1.55)1.09 (0.86–1.37)
ln PIB per cápita2.63 (1.56–4.42)†2.63 (1.56–4.42)†3.64 (1.84–7.22)†2.95 (1.39–6.26)†11.64 (4.53–29.95)†30.27 (10.94–83.77)†
I+D (% PIB)1.59 (1.13–2.25)†1.59 (1.13–2.25)†1.68 (1.10–2.58)†1.78 (1.14–2.78)†2.83 (1.64–4.87)†1.57 (1.02–2.42)†
Matrícula terciaria (por 10 pp)1.21 (1.02–1.44)†1.21 (1.02–1.44)†1.24 (1.01–1.53)†1.20 (0.97–1.48)1.05 (0.80–1.40)0.79 (0.62–1.02)
Inflación: ln publicaciones (γ₁)—Colapsa en BN—Separación casi completa—-1.86 (p = 0.003)
α (dispersión)0.78—0.94—0.98—
Hurdle, parte logística: OR chocolate1.11 (0.68–1.81)p = 0.6851.24 (0.69–2.22)p = 0.4682.07 (0.95–4.51)p = 0.067

Nota: † p < 0.05. ln PIB per cápita centrado en 9 (≈ USD 8100); la RTI expresa el cambio por unidad logarítmica (≈ 2.7 veces el PIB). "Colapsa en BN": la probabilidad de cero estructural convergió a 0 y la ZINB se redujo a la BN. "Separación casi completa": γ₁ divergió (los ceros estructurales se concentran en los países con menor producción científica), por lo que no tiene error estándar finito. OR = odds ratio de tener al menos un laureado.

Fuente: elaboración propia.

Fig. 3 Razón de tasas de incidencia (RTI, escala logarítmica) del consumo de cacao en especificaciones anidadas (M1–M5) y tres familias de modelos de conteo para los tres desenlaces. La línea discontinua marca RTI = 1 (sin asociación).
Fig. 3. Razón de tasas de incidencia (RTI, escala logarítmica) del consumo de cacao en especificaciones anidadas (M1–M5) y tres familias de modelos de conteo para los tres desenlaces. La línea discontinua marca RTI = 1 (sin asociación).
Fuente: elaboración propia.

3.4. Sensibilidad y control de falsificación

La ausencia de asociación ajustada fue estable al excluir a Suecia y Suiza, a Estados Unidos y Reino Unido, al usar el PIB de 1990–2022 y al incluir países de cualquier tamaño (RTI entre 0.95 y 1.21; todos los IC incluyeron 1). La única señal residual apareció al excluir a los países productores de cacao y a los cuatro centros europeos de procesamiento y reexportación (19 países en total), donde la BN para ciencias por afiliación arrojó una RTI de 1.54 (1.02–2.32; p = 0.041), que no se confirmó en la ZINB (1.25; 0.84–1.86). El control de falsificación resultó informativo: el suministro de leche mantuvo una asociación ajustada significativa con todos los premios (RTI = 1.010 por kg; IC 95%: 1.006–1.014; equivalente a 2.7 veces por cada 100 kg/hab/año) que persistió incluso tras añadir un indicador de Europa (1.005; 1.002–1.009; p = 0.004), mientras que el chocolate no lo hizo (1.10; 0.93–1.31; p = 0.25). El indicador de Europa, por sí mismo, multiplicó por 4.7 la tasa esperada de laureados (Tabla 6).

Tabla 6. Análisis de sensibilidad y control de falsificación: RTI (IC 95%) del chocolate en el modelo M3
EscenarionTodos (nac.) BNTodos (nac.) ZINBCiencias (afil.) BNCiencias (afil.) ZINB
Base1291.11 (0.89–1.38)1.11 (0.89–1.38)1.15 (0.85–1.55)1.09 (0.86–1.37)
Excl. productores y hubs1101.28 (0.96–1.70)1.28 (0.97–1.71)1.54 (1.02–2.32)1.25 (0.84–1.86)
Excl. Suecia y Suiza1271.14 (0.92–1.42)1.14 (0.92–1.42)1.21 (0.88–1.67)1.12 (0.88–1.43)
Excl. EE. UU. y Reino Unido1271.08 (0.86–1.36)1.08 (0.86–1.36)1.00 (0.73–1.36)0.95 (0.76–1.19)
PIB 1990-20221291.13 (0.91–1.41)1.13 (0.91–1.41)1.17 (0.85–1.62)1.10 (0.86–1.42)
Todas las poblaciones1361.17 (0.95–1.43)1.17 (0.95–1.43)1.12 (0.83–1.49)1.04 (0.84–1.30)
+ indicador Europa (chocolate)1291.10 (0.93–1.31)—1.19 (0.96–1.46)1.14 (0.94–1.39)
Falsificación: leche en lugar de chocolate (por kg)1211.010 (1.006–1.014)1.010 (1.007–1.014)——
Falsificación: leche + indicador Europa (por kg)1211.005 (1.002–1.009)1.005 (1.002–1.009)1.004 (0.999–1.010)1.003 (0.998–1.009)

Nota: "Excl. productores y hubs" excluye a productores de cacao y a Bélgica, Países Bajos, Suiza y Alemania (reexportadores). En las filas de falsificación la RTI corresponde a la leche (por kg/hab/año). (—) no estimable o no aplicable.

Fuente: elaboración propia.

3.5. América Latina y el Caribe: atribución por nacimiento y afiliación

Entre 1901 y 2025, 21 laureados nacieron en ALC: ocho en Literatura, siete en Paz, cinco en ciencias y uno en Economía (Tabla 7; Fig. 4A). De los cinco científicos, solo Bernardo Houssay recibió el premio trabajando en su país de nacimiento; César Milstein (Reino Unido), Baruj Benacerraf y Mario Molina (Estados Unidos) y Peter Medawar (Reino Unido) fueron premiados en el extranjero. El único otro premio científico con afiliación latinoamericana, el de Luis F. Leloir (Argentina, 1970), corresponde a un investigador nacido en Francia. Tras el ajuste multivariado, pertenecer a ALC se asoció con una tasa de laureados en ciencias por nacimiento un 81% menor que la esperada según su riqueza, I+D y educación (RTI = 0.19; IC 95%: 0.05–0.63; p = 0.007), mientras que para todos los premios la diferencia no fue significativa (RTI = 0.61; 0.26–1.42). A escala global, Estados Unidos concentró 320 premios científicos por afiliación frente a 217 por nacimiento (saldo neto de +103), mientras que Polonia, Austria, Hungría y China registraron los mayores saldos negativos (Fig. 4B).

Tabla 7. Laureados Nobel nacidos en América Latina y el Caribe (1901–2025) y país de afiliación al recibir el premio
LaureadoAñoCategoríaPaís de nacimientoAfiliación al recibir el premio
Luis F. Leloir*1970QuímicaFranciaArgentina
Mario J. Molina1995QuímicaMéxicoEstados Unidos
Bernardo A. Houssay1947Fisiología o MedicinaArgentinaArgentina
Peter B. Medawar1960Fisiología o MedicinaBrasilReino Unido
Baruj Benacerraf1980Fisiología o MedicinaVenezuelaEstados Unidos
César Milstein1984Fisiología o MedicinaArgentinaReino Unido
W. Arthur Lewis1979Ciencias EconómicasSanta LucíaEstados Unidos
Gabriela Mistral1945LiteraturaChileNo aplica
Miguel Ángel Asturias1967LiteraturaGuatemalaNo aplica
Pablo Neruda1971LiteraturaChileNo aplica
Gabriel García Márquez1982LiteraturaColombiaNo aplica
Octavio Paz1990LiteraturaMéxicoNo aplica
Derek Walcott1992LiteraturaSanta LucíaNo aplica
V. S. Naipaul2001LiteraturaTrinidad y TobagoNo aplica
Mario Vargas Llosa2010LiteraturaPerúNo aplica
Carlos Saavedra Lamas1936PazArgentinaNo aplica
Adolfo Pérez Esquivel1980PazArgentinaNo aplica
Alfonso García Robles1982PazMéxicoNo aplica
Oscar Arias Sánchez1987PazCosta RicaNo aplica
Rigoberta Menchú Tum1992PazGuatemalaNo aplica
Juan Manuel Santos2016PazColombiaNo aplica
María Corina Machado2025PazVenezuelaNo aplica

Nota: * incluido como referencia; nacido en París y nacionalizado argentino, no se contabiliza como nacido en ALC. "No aplica": categorías sin afiliación institucional.

Fuente: elaboración propia con datos de Nobel Prize Outreach (2025).

Fig. 4 (A) Laureados Nobel nacidos en América Latina y el Caribe según país y categoría del premio (1901–2025). (B) Saldo neto de laureados en ciencias (premios por país de afiliación menos premios por país de nacimiento) en países seleccionados; en naranja, los países latinoamericanos. La barra de Estados Unidos está truncada (+103).
Fig. 4. (A) Laureados Nobel nacidos en América Latina y el Caribe según país y categoría del premio (1901–2025). (B) Saldo neto de laureados en ciencias (premios por país de afiliación menos premios por país de nacimiento) en países seleccionados; en naranja, los países latinoamericanos. La barra de Estados Unidos está truncada (+103).
Fuente: elaboración propia con datos de Nobel Prize Outreach (2025).

3.6. Contraste con las replicaciones previas

La Tabla 8 sitúa estos resultados frente a los principales estudios que analizaron la hipótesis chocolate–Nobel. La correlación original solo se observó con cifras industriales en una muestra seleccionada; las críticas posteriores la reprodujeron con otras variables de prosperidad, y los análisis multivariados coinciden en que la actividad científica y la riqueza son los predictores dominantes. El presente estudio es el primero, hasta donde se conoce, que combina una exposición armonizada, una muestra que incluye a ALC, la desagregación por nacimiento y afiliación y modelos de conteo con exceso de ceros.

Tabla 8. Contraste entre el presente estudio y los principales análisis previos de la hipótesis chocolate–Nobel
EstudioMuestraMedida de exposiciónControlesMétodoHallazgo principal
Messerli (2012)23 países, mayoritariamente ricos y con laureadosVentas de chocolate según fuentes industriales, años distintos (2002–2011)NingunoCorrelación de Pearsonr = 0.791 (p < 0.0001); hipótesis de los flavanoles
Maurage et al. (2013)Países del estudio originalChocolate y otros alimentos ricos en flavanolesAlimentos alternativosCorrelaciones comparadasLa explicación por flavanoles se debilita; riesgo de sobreinterpretación
Linthwaite & Fuller (2013); Dunstan (2013)Muestras similaresConsumo de leche—Correlación; crítica metodológicaCorrelación análoga con la leche; diseño ecológico y desfase temporal
Li (2013)Comparación transnacionalChocolate y lechePIB per cápitaCorrelaciónLa riqueza es la causa común plausible
Doi et al. (2014)220 paísesHábitos dietéticos y factores económicosPublicaciones, I+D, PIBPartición jerárquicaPublicaciones e I+D predicen mejor que la dieta
Prinz (2020)70 paísesConsumo de chocolatePIB, I+D, publicaciones, café, tabaco, alcoholSelección de Heckman en dos etapasAsociación residual positiva tras ajuste; sin interpretación causal
Presente estudio129 países (19 de ALC); 991 premios, 1901–2025Cacao disponible, FAOSTAT armonizado (2015–2022)PIB histórico 1960–2022, I+D, matrícula terciaria, publicaciones, ALC, EuropaCorrelación parcial; Poisson, BN, ZIP, ZINB, hurdleSin réplica con datos armonizados (r = 0.19); RTI ajustada 1.09–1.15 (IC incluye 1); residuos inestables, comparables a la leche

Fuente: elaboración propia a partir de los estudios citados.

4. Discusión

La respuesta a la pregunta de investigación es clara: con datos actualizados a 2025, una muestra ampliada que incluye a América Latina y modelos adecuados a la naturaleza de la variable respuesta, el consumo de chocolate no se asocia de manera independiente con la producción de laureados Nobel. La correlación cruda existe en la muestra ampliada, pero es más débil que la original y desaparece en cuanto se ajusta por la riqueza histórica de los países. Este patrón es el que predice el grafo causal de la Fig. 1, en el que el PIB actúa como causa común del consumo de un bien de lujo y de la infraestructura científica (Li, 2013), y concuerda con la conclusión de Doi et al. (2014) de que la actividad científica predice los premios mejor que la dieta.

Un primer hallazgo relevante es que la correlación de Messerli (2012) no se reproduce cuando la exposición se mide con una fuente armonizada. En los mismos 23 países, el coeficiente cae de 0.791 a 0.19. Esta discrepancia confirma las advertencias sobre la debilidad de los datos de consumo: las cifras industriales empleadas en 2012 procedían de años y fuentes distintas (Apgar & Tarka, 2019; Del Prete & Samoggia, 2020), mientras que el balance de FAOSTAT mide disponibilidad aparente y distorsiona a los países procesadores y reexportadores, como Suiza, Bélgica o los Países Bajos, y a los productores, como Côte d'Ivoire. Que la magnitud de la correlación dependa tanto de la fuente elegida es, en sí mismo, un indicio de fragilidad: una relación causal robusta no debería desaparecer al cambiar un indicador de exposición por otro razonable.

El segundo hallazgo concierne a la elección del modelo. La correlación de Pearson trata como continua y normal una variable que es un conteo con un 46–80% de ceros. Los modelos de Poisson y ZIP, que ignoran la sobredispersión, atribuyeron significación al chocolate incluso tras el ajuste, un error que habría conducido a una falsa confirmación. Una vez modelada la sobredispersión mediante la BN, la ZINB o el modelo hurdle, la asociación se volvió nula. La ZINB aportó además una lectura sustantiva: el exceso de ceros en los premios de ciencias por afiliación está gobernado por la escasez de producción científica, lo que coincide con la ecuación de selección de Prinz (2020), en la que el volumen de publicaciones discrimina qué países entran en el mapa Nobel. Para todos los premios, en cambio, el componente de inflación fue innecesario, porque la Literatura y la Paz distribuyen laureados en países sin sistemas científicos consolidados.

Los resultados matizan, sin contradecirlo del todo, el hallazgo de Prinz (2020) sobre una asociación residual. En nuestro análisis aparecieron dos señales de ese tipo: una RTI significativa en la BN tras excluir productores y reexportadores, y una tendencia en la parte logística del modelo hurdle para la afiliación en ciencias. Ninguna fue consistente entre familias de modelos y ambas surgieron en subconjuntos con pocos países positivos. La prueba de falsificación permite interpretarlas: la leche, una variable sin plausibilidad biológica para la excelencia científica, conservó una asociación ajustada más fuerte que la del chocolate, incluso tras controlar por Europa. Esto indica que los modelos ecológicos con unas pocas covariables macroeconómicas dejan confusión residual de tipo geocultural, relacionada con la tradición académica del norte de Europa, la estabilidad institucional y la dieta láctea. Cualquier residuo del chocolate debe leerse, por tanto, como heterogeneidad no observada y no como efecto de los flavanoles, tal como ya habían advertido Maurage et al. (2013) y Pollet et al. (2014).

La plausibilidad biológica tampoco respalda la hipótesis. Si 500 mg diarios de flavanoles purificados durante tres años no mejoraron la cognición global de 2262 adultos (Baker et al., 2023; Vyas et al., 2024), resulta inverosímil que el consumo poblacional de chocolate comercial, con bajo contenido de flavanoles y medido décadas después de la formación de los laureados, explique diferencias nacionales en la producción de premios. Los efectos agudos o específicos de dominio descritos en ensayos pequeños y metaanálisis (Brickman et al., 2014, 2023; Cheng et al., 2022; Arabpour et al., 2026) operan en escalas individuales que no pueden trasladarse a un agregado nacional sin incurrir en la falacia ecológica (Robinson, 1950; Loney & Nagelkerke, 2014; Szromek, 2026).

La inclusión de América Latina aporta la contribución más original. La región no se diferencia del resto del mundo en consumo aparente de cacao ni en PIB per cápita histórico, pero invierte la mitad en I+D y produce menos ciencia per cápita. Sus laureados se concentran en Literatura y Paz, y en ciencias presenta una tasa un 81% menor que la esperada por sus covariables. La desagregación por afiliación muestra el mecanismo: cuatro de los cinco científicos latinoamericanos fueron premiados en el Reino Unido o Estados Unidos, lo que reproduce a escala regional la concentración del trabajo premiado en unos pocos sistemas científicos descrita por Schlagberger et al. (2016) y Heinze et al. (2020). Ecuador y Perú, grandes productores de cacao fino, registran una disponibilidad de cacao baja o moderada y ningún premio científico, y Brasil solo cuenta con Medawar, nacido allí pero formado y premiado en el Reino Unido; es el contraejemplo natural que la literatura señalaba como ausente (Rizzuto et al., 2017; Cortés & Andrade, 2022). Para la política científica regional, la implicación es directa: lo que distingue a los países con laureados no es la dieta, sino la inversión sostenida en investigación, la capacidad de retener talento y la densidad de su producción científica.

Este estudio tiene limitaciones. Primero, como todo diseño ecológico, no permite inferencias individuales ni causales y su potencia es limitada por el número de países, en especial en los modelos con pocos conteos positivos. Segundo, la exposición de FAOSTAT es contemporánea (2015–2022), mientras que los premios abarcan más de un siglo; no se dispuso de series armonizadas de consumo de chocolate anteriores a 2010, por lo que el desfase temporal señalado por Dunstan (2013) solo se corrigió en las covariables. Tercero, la medida de disponibilidad aparente sobrestima o subestima el consumo real en países productores y reexportadores, lo cual se abordó parcialmente mediante el análisis de sensibilidad. Cuarto, la atribución por país de nacimiento según fronteras actuales y por primera afiliación declarada simplifica trayectorias complejas, y no se codificó el país de formación doctoral. Quinto, las covariables de I+D y matrícula tienen series incompletas para algunos países y periodos, lo que obligó a usar promedios de cortes disponibles. Finalmente, la presencia de confusión residual demostrada por el control de leche implica que ningún modelo ecológico de pocas variables puede excluir por completo la heterogeneidad no observada.

5. Conclusiones

Con datos actualizados a 2025 y una muestra de 129 países que incluye a 19 de América Latina y el Caribe, la correlación entre consumo de chocolate y laureados Nobel no se reproduce con una medida de exposición armonizada y desaparece al ajustar por la riqueza histórica. Los modelos de binomial negativa, binomial negativa inflada por ceros y hurdle coinciden en que el PIB per cápita, el gasto en I+D y la producción científica, y no el chocolate, explican la distribución de los premios. Las señales residuales del chocolate son inestables y de magnitud comparable o menor que las de un control de falsificación como la leche, lo que apunta a confusión geocultural y no a un efecto biológico. En América Latina, el rezago en premios científicos se asocia con una baja inversión en I+D y con la emigración de sus científicos más destacados. Como líneas futuras se proponen: (a) modelos de panel con desfase temporal que relacionen series históricas de consumo, PIB e I+D con premios posteriores; (b) la codificación del país de formación doctoral y de trabajo premiado para cuantificar la fuga de cerebros; (c) modelos de ecuaciones estructurales que estimen simultáneamente las rutas riqueza → I+D → premios y riqueza → consumo de bienes de lujo, y (d) la validación cruzada de las medidas de consumo de cacao con encuestas de hogares en países latinoamericanos productores, como Ecuador.

Agradecimientos

La autora agradece a la Universidad Estatal Amazónica por el apoyo institucional y a los equipos de FAOSTAT, UNESCO-UIS, el Banco Mundial, el Maddison Project y Our World in Data por mantener datos abiertos que hicieron posible esta replicación.

Contribución de los autores: S. L. R. R.: Conceptualización, curación de datos, metodología, análisis formal, visualización, redacción—borrador original, redacción—revisión y edición. La autora ha leído y aprobado la versión publicada del manuscrito.

Financiamiento: Los autores no recibieron apoyo financiero para la investigación, redacción y/o publicación de este artículo.

Conflictos de intereses: Los autores declaran no tener conflictos de interés.

Declaración de uso de inteligencia artificial: Se utilizó Claude (Anthropic; modelo claude-opus-5-5, 2026) como asistente para la recopilación de datos abiertos, la programación de los análisis estadísticos en Python, la elaboración de figuras y la redacción del borrador del manuscrito. La autora revisó, verificó y asume la responsabilidad plena sobre el contenido, los datos, los análisis y las conclusiones.

Referencias

Apgar, J., & Tarka, S. (2019). Methylxanthine composition and consumption patterns of cocoa and chocolate products. En G. A. Spiller (Ed.), Caffeine. CRC Press. https://doi.org/10.1201/9780429126789-7

Arabpour, Z., Arabpour, Z., & Jazayeri, S. (2026). Effect of chronic dark chocolate or cocoa intake on cognitive function: A GRADE-assessed systematic review and meta-analysis of randomized controlled trials. Nutrition Reviews. https://doi.org/10.1093/nutrit/nuaf286

Baker, L. D., Manson, J. E., Rapp, S. R., Sesso, H. D., Gaussoin, S. A., Shumaker, S. A., & Espeland, M. A. (2023). Effects of cocoa extract and a multivitamin on cognitive function: A randomized clinical trial. Alzheimer's & Dementia, 19(4), 1308–1319. https://doi.org/10.1002/alz.12767

Banco Mundial. (2025). Scientific and technical journal articles (IP.JRN.ARTC.SC) [Conjunto de datos]. World Development Indicators. https://ourworldindata.org/grapher/scientific-and-technical-journal-articles

Bolt, J., & van Zanden, J. L. (2024). Maddison-style estimates of the evolution of the world economy: A new 2023 update. Journal of Economic Surveys. https://doi.org/10.1111/joes.12618

Brickman, A. M., Khan, U. A., Provenzano, F. A., Yeung, L.-K., Suzuki, W., Schroeter, H., Wall, M., Sloan, R. P., & Small, S. A. (2014). Enhancing dentate gyrus function with dietary flavanols improves cognition in older adults. Nature Neuroscience, 17(12), 1798–1803. https://doi.org/10.1038/nn.3850

Brickman, A. M., Yeung, L.-K., Alschuler, D. M., Ottaviani, J. I., Kuhnle, G. G. C., Sloan, R. P., Luttmann-Gibson, H., Copeland, T., Schroeter, H., Sesso, H. D., Manson, J. E., Wall, M., & Small, S. A. (2023). Dietary flavanols restore hippocampal-dependent memory in older adults with lower diet quality and lower habitual flavanol consumption. Proceedings of the National Academy of Sciences, 120(23), e2216932120. https://doi.org/10.1073/pnas.2216932120

Cameron, A. C., & Trivedi, P. K. (2013). Regression analysis of count data (2.ª ed.). Cambridge University Press. https://doi.org/10.1017/CBO9781139013567

Cheng, N., Bell, L., Lamport, D. J., & Williams, C. M. (2022). Dietary flavonoids and human cognition: A meta-analysis. Molecular Nutrition & Food Research, 66, e2100976. https://doi.org/10.1002/mnfr.202100976

Cortés, J. D., & Andrade, D. A. (2022). The Colombian scientific elite—Science mapping and a comparison with Nobel Prize laureates using a composite citation indicator. PLoS ONE, 17(6), e0269116. https://doi.org/10.1371/journal.pone.0269116

Del Prete, M., & Samoggia, A. (2020). Chocolate consumption and purchasing behaviour review: Research issues and insights for future research. Sustainability, 12(14), 5586. https://doi.org/10.3390/su12145586

Desideri, G., Kwik-Uribe, C., Grassi, D., Necozione, S., Ghiadoni, L., Mastroiacovo, D., Raffaele, A., Ferri, L., Bocale, R., Lechiara, M. C., Marini, C., & Ferri, C. (2012). Benefits in cognitive function, blood pressure, and insulin resistance through cocoa flavanol consumption in elderly subjects with mild cognitive impairment: The Cocoa, Cognition, and Aging (CoCoA) Study. Hypertension, 60(3), 794–801. https://doi.org/10.1161/HYPERTENSIONAHA.112.193060

Doi, H., Heeren, A., & Maurage, P. (2014). Scientific activity is a better predictor of Nobel award chances than dietary habits and economic factors. PLoS ONE, 9(3), e92612. https://doi.org/10.1371/journal.pone.0092612

Dunstan, F. (2013). Nobel prizes, chocolate and milk: The statistical view. Practical Neurology, 13(3), 206–207. https://doi.org/10.1136/practneurol-2013-000588

FAO. (2025). FAOSTAT: Food balances (2010–) [Conjunto de datos]. Organización de las Naciones Unidas para la Alimentación y la Agricultura. https://www.fao.org/faostat/en/#data/FBS

Heckman, J. J. (1979). Sample selection bias as a specification error. Econometrica, 47(1), 153–161. https://doi.org/10.2307/1912352

Heinze, T., & Fuchs, J. E. (2022). National and organizational patterns of Nobel laureate careers in physiology/medicine, physics, and chemistry. Scientometrics, 127, 7273–7288. https://doi.org/10.1007/s11192-021-04250-0

Heinze, T., von der Heyden, M., & Pithan, D. (2020). Institutional environments and breakthroughs in science. Comparison of France, Germany, the United Kingdom, and the United States. PLoS ONE, 15(9), e0239805. https://doi.org/10.1371/journal.pone.0239805

Lambert, D. (1992). Zero-inflated Poisson regression, with an application to defects in manufacturing. Technometrics, 34(1), 1–14. https://doi.org/10.2307/1269547

Li, J. (2013). Economy and Nobel prizes: Cause behind chocolate and milk? Practical Neurology, 14(1), e1. https://doi.org/10.1136/practneurol-2013-000545

Li, Q.-F. (2025). Assessing and adjusting for bias in ecological analysis using multiple sample datasets. BMC Medical Research Methodology, 25. https://doi.org/10.1186/s12874-025-02552-y

Linthwaite, S., & Fuller, G. N. (2013). Milk, chocolate and Nobel prizes. Practical Neurology, 13(1), 63. https://doi.org/10.1136/practneurol-2012-000471

Loney, T., & Nagelkerke, N. J. (2014). The individualistic fallacy, ecological studies and instrumental variables: A causal interpretation. Emerging Themes in Epidemiology, 11, 18. https://doi.org/10.1186/1742-7622-11-18

Mastroiacovo, D., Kwik-Uribe, C., Grassi, D., Necozione, S., Raffaele, A., Pistacchio, L., Righetti, R., Bocale, R., Lechiara, M. C., Marini, C., Ferri, C., & Desideri, G. (2015). Cocoa flavanol consumption improves cognitive function, blood pressure control, and metabolic profile in elderly subjects: The Cocoa, Cognition, and Aging (CoCoA) Study—A randomized controlled trial. The American Journal of Clinical Nutrition, 101(3), 538–548. https://doi.org/10.3945/ajcn.114.092189

Maurage, P., Heeren, A., & Pesenti, M. (2013). Does chocolate consumption really boost Nobel award chances? The peril of over-interpreting correlations in health studies. The Journal of Nutrition, 143(6), 931–933. https://doi.org/10.3945/jn.113.174813

McClintock, S., Stangl, D., & Çetinkaya-Rundel, M. (2014). Taking a chance in the classroom: The real secret to genius? Reading between the lines. CHANCE, 27(1), 41–44. https://doi.org/10.1080/09332480.2014.890869

Messerli, F. H. (2012). Chocolate consumption, cognitive function, and Nobel laureates. The New England Journal of Medicine, 367(16), 1562–1564. https://doi.org/10.1056/NEJMon1211064

Mullahy, J. (1986). Specification and testing of some modified count data models. Journal of Econometrics, 33(3), 341–365. https://doi.org/10.1016/0304-4076(86)90002-3

Nobel Prize Outreach. (2025). All Nobel Prizes [Base de datos]. NobelPrize.org. https://www.nobelprize.org/prizes/lists/all-nobel-prizes/

Our World in Data. (2025a). Cocoa bean consumption per person [Conjunto de datos, a partir de FAOSTAT]. https://ourworldindata.org/grapher/chocolate-consumption-per-person

Our World in Data. (2025b). Research & development spending as a share of GDP [Conjunto de datos, a partir de UNESCO-UIS y Banco Mundial]. https://ourworldindata.org/grapher/research-spending-gdp

Our World in Data. (2025c). Gross enrollment ratio in tertiary education [Conjunto de datos, a partir de UNESCO-UIS]. https://ourworldindata.org/grapher/gross-enrollment-ratio-in-tertiary-education

Pietrucha, J. (2018). Country-specific determinants of world university rankings. Scientometrics, 114, 1129–1139. https://doi.org/10.1007/s11192-017-2634-1

Pollet, T. V., Tybur, J. M., Frankenhuis, W. E., & Rickard, I. J. (2014). What can cross-cultural correlations teach us about human nature? Human Nature, 25(3), 410–429. https://doi.org/10.1007/s12110-014-9206-3

Prinz, A. (2020). Chocolate consumption and Noble laureates [Documento de trabajo]. SSRN. https://doi.org/10.2139/ssrn.3562373

rfordatascience. (2019). Nobel Prize winners [Conjunto de datos, TidyTuesday 2019-05-14]. GitHub. https://github.com/rfordatascience/tidytuesday/tree/master/data/2019/2019-05-14

Rizzuto, L., Anido, J. D., & Azuaje, R. (2017). El consumo de cacao en Venezuela y en el mundo desde una perspectiva sostenible (1960–2014) (pp. 23–49).

Robinson, W. S. (1950). Ecological correlations and the behavior of individuals. American Sociological Review, 15(3), 351–357. https://doi.org/10.2307/2087176

Rodríguez-Navarro, A., & Brito, R. (2022). The link between countries' economic and scientific wealth has a complex dependence on technological activity and research policy. Scientometrics, 127, 2871–2896. https://doi.org/10.1007/s11192-022-04313-w

Schlagberger, E. M., Bornmann, L., & Bauer, J. (2016). At what institutions did Nobel laureates do their prize-winning work? An analysis of biographical information on Nobel laureates from 1994 to 2014. Scientometrics, 109, 723–767. https://doi.org/10.1007/s11192-016-2059-2

Seabold, S., & Perktold, J. (2010). Statsmodels: Econometric and statistical modeling with Python. En Proceedings of the 9th Python in Science Conference (pp. 92–96). https://doi.org/10.25080/Majora-92bf1922-011

Szromek, A. (2026). Cross-level inference errors in scientific interpretation—From ecological correlation to the control of inference and communication levels. Publications, 14, 41. https://doi.org/10.3390/publications14030041

Vuong, Q. H. (1989). Likelihood ratio tests for model selection and non-nested hypotheses. Econometrica, 57(2), 307–333. https://doi.org/10.2307/1912557

Vyas, C. M., Manson, J. E., Sesso, H. D., Rist, P. M., Weinberg, A., Kim, E., Moorthy, M. V., Cook, N. R., & Okereke, O. I. (2024). Effect of cocoa extract supplementation on cognitive function: Results from the clinic subcohort of the COSMOS trial. The American Journal of Clinical Nutrition, 119(1), 39–48. https://doi.org/10.1016/j.ajcnut.2023.10.031