Глобальная заболеваемость раком щитовидной железы растёт с выраженной географической неоднородностью, которую не объяснить только скринингом. Экологическое исследование с использованием машинного обучения интегрировало данные GLOBOCAN 2022, FAOSTAT (113 продуктовых категорий, 2010–2022 гг.), индекс человеческого развития (ИЧР) ООН и распространённость ожирения ВОЗ по 163 странам (155 переменных, 150 числовых предикторов) для выявления ассоциаций между национальной доступностью продуктов питания и возрастно-стандартизованной заболеваемостью.
Стекинг-ансамбль показал лучшую производительность в предсказании заболеваемости: R² = 0.613 (95% ДИ [0.510, 0.726]), RMSE = 3.04, MAE = 1.82, ρ Спирмена = 0.901 (95% ДИ [0.862, 0.924]). Случайный лес (R² = 0.609), бэггинг (R² = 0.599) и KNN (R² = 0.580) следовали за ним; глубокая табличная модель DynaTab была сопоставима (R² = 0.626). Однако предсказательная точность почти полностью определялась социально-экономическими факторами: девять переменных развития дали R² = 0.609, 140 диетических переменных — только 0.284, а инкрементный вклад диетического блока составил –0.048 (95% ДИ от –0.169 до +0.034) в основной модели. При валидации с исключением регионов (leave-one-region-out) ранжирование стран сохранилось (ρ = 0.813), но абсолютный уровень — нет.
В нескорректированном анализе положительные корреляции с заболеваемостью показали индекс доступности йодсодержащих продуктов (ρ = 0.234), общая калорийность (ρ = 0.230), морепродукты (ρ = 0.219) и соотношение животного/растительного белка (ρ = 0.221). После поправки на ИЧР, ожирение, субиндексы ИЧР и региональные эффекты, с коррекцией Бенджамини–Хохберга, доступность морепродуктов (скорректированная r = 0.196, q = 0.031), молочных продуктов (r = 0.198, q = 0.017), индекс йодсодержащих продуктов (r = 0.224, q = 0.026) и соотношение животного/растительного белка остались значимыми; растительный белок — нет (r = –0.095, q = 0.309).
В серии чувствительных анализов ассоциация индекса определялась морепродуктами (r = 0.196, q = 0.017) и молочными продуктами (r = 0.198, q = 0.017), но не яйцами (r = 0.065, q = 0.411); была устойчива к различным схемам взвешивания (ρ = 0.238–0.256), сохранялась при дополнительной поправке на общее потребление животного белка (r = 0.203, q = 0.032) и не воспроизводилась индексом йод-бедных животных продуктов (r = 0.102, q = 0.198). Декомпозиция пути показала, что 50.2% корреляции между морепродуктами и заболеваемостью (95% бутстрэп ДИ 16%–90%) объясняется путём через ИЧР.
Во вторичном анализе смертности диетическое ранжирование сохранилось, но ассоциации с соотношением смертность/заболеваемость были близки к нулю (ИЧР ρ = 0.053; морепродукты ρ = 0.065), что согласуется со значительным вкладом дифференциальной детекции (скрининг-эффект). В классификационном анализе Extra Trees достиг AUC = 0.932 (95% ДИ [0.893, 0.964]).