Preview

Проблемы Эндокринологии

Расширенный поиск

Математические модели прогнозирования конечного роста и его коэффициента стандартного отклонения у детей с дефицитом гормона роста в российской популяции

Содержание

Перейти к:

Аннотация

Актуальность. Прогнозирование конечного роста (КР), а следовательно и эффективности длительной терапии рекомбинантным гормоном роста (рГР) у пациентов с СТГ-дефицитом, — важное условие работы клинициста. Это позволяет осуществлять персонализированный подход к лечению пациентов с СТГ-дефицитом: контролировать регулярность введения и правильность дозирования препарата и идентифицировать факторы, от которых зависит КР.


Цель исследования — разработать математические модели прогнозирования КР и его коэффициента стандартного отклонения (standard deviation score — SDS) у пациентов с СТГ-дефицитом в российской популяции.


Материал и методы. Для моделирования использованы данные 121 пациента, наблюдавшегося в Институте детской эндокринологии ФГБУ ЭНЦ в период с 1978 по 2016 г. с диагнозом СТГ-дефицит и получавшего рГР от момента диагностики до достижения КР. В качестве предикторов использовались 11 признаков: пол, хронологический возраст на момент диагностики СТГ-дефицита, пубертатный статус, форма заболевания, регулярность терапии рГР, SDS роста при рождении, SDS роста на момент диагностики СТГ-дефицита, костный возраст на момент диагностики СТГ-дефицита, коэффициент костный возраст/хронологический, SDS генетически прогнозируемого роста, максимальный уровень СТГ в ходе теста с клофелином. Для построения моделей применялись методы множественной линейной регрессии и искусственные нейронные сети (ИНС).


Результаты. Разработанные ИНС продемонстрировали высокую точность прогнозирования КР (среднеквадратическая ошибка — 4,4 см, доля объяснимой дисперсии — 76%) и меньшую точность прогнозирования SDS КР (среднеквадратическая ошибка — 0,601, доля объяснимой дисперсии — 42%). Линейные регрессионные модели, построенные только на количественных предикторах, имели существенно худшее качество. Для наилучшей ИНС разработана программная реализация, доступная в Интернет.


Заключение. Программно реализованная модель прогнозирования КР на основе ИНС использует в качестве предикторов показатели, доступные для любого врача-клинициста, и может применяться для индивидуального прогнозирования КР. В перспективе анализ более крупных баз данных позволит улучшить качество прогнозирования эффективности терапии рГР.

Для цитирования:


Гаврилова А.Е., Нагаева Е.В., Реброва О.Ю., Ширяева Т.Ю. Математические модели прогнозирования конечного роста и его коэффициента стандартного отклонения у детей с дефицитом гормона роста в российской популяции. Проблемы Эндокринологии. 2017;63(5):282-290.

For citation:


Gavrilova A.E., Nagaeva E.V., Rebrova O.Yu., Shiryaeva T.Yu. Mathematical models for predicting final achieved height and its standard deviation rate in children with growth hormone deficiency in the Russian population. Problems of Endocrinology. 2017;63(5):282-290.

Соматотропная недостаточность (СТГ-дефицит) — патологическое состояние, вызванное нарушением синтеза, секреции, регуляции или биологического эффекта СТГ. Выделяют изолированный СТГ-дефицит (ИДГР) и множественный дефицит гормонов аденогипофиза (МДГА). При ИДГР единственным компонентом гипопитуитаризма является дефицит ГР, а при МДГА, помимо СТГ-дефицита, имеет место недостаточность еще одного или нескольких тропных гормонов гипофиза (адренокортикотропного гормона (АКТГ), тиреотропного гормона (ТТГ), лютеинизирующего гормона (ЛГ), фолликулостимулирующего гормона (ФСГ), пролактина).

С 1985 г. общепризнанным средством нормализации темпов роста и физического развития ребенка, обусловленной СТГ-дефицитом, является рекомбинантный гормон роста (рГР), высокая эффективность которого многократно доказана [1—4]. В то же время ответ на лечение у детей значительно различается, а терапия является длительной и дорогостоящей. В последние годы появилась возможность анализировать результаты лечения нескольких поколений детей, длительно получавших терапию рГР и достигших конечного роста (КР), что дает основания для поиска предикторов эффективности лечения. Особенно актуальным является персонифицированный подход к назначению рГР, оптимизация эффекта и стоимости терапии.

Уже на заре применения рГР исследователи пытались разработать математические модели индивидуального прогнозирования эффективности лечения СТГ-дефицита, которые можно использовать на ранних этапах терапии рГР для расчета долговременного ростового ответа и прогнозирования КР пациента. В настоящее время разработано несколько математических моделей, две из которых прогнозируют скорость роста в первый год терапии [5, 7], одна — суммарный ростовой эффект от начала пубертата до достижения КР [8] и три — КР и коэффициент его стандартного отклонения (standard deviation score — SDS) [6, 9, 10]. В качестве предикторов в этих моделях использовались как доступные к началу терапии ауксологические и лабораторные показатели, так и данные, получаемые после начала терапии — длина голени после 3 мес лечения [5], SDS скорости роста до и в первые два года лечения, уровень дезоксипиридинолина (маркера резорбции кости) в моче через 1 мес после начала терапии, скорость роста через 3 мес после начала терапии [7]. Доля объяснимой моделями дисперсии варьирует от 37 [9] до 89% [7]. Основными недостатками разработанных ранее моделей являются ограниченные горизонты прогнозирования, недостаточная точность и необходимость использования специфических предикторов, не всегда доступных в рутинной практике.

Цель исследования — разработать математические модели прогнозирования КР и его SDS у пациентов с СТГ-дефицитом в российской популяции.

Материал и методы

Дизайн исследования

Проанализированы данные 141 пациента (64% мальчики, n=90), наблюдавшихся в Институте детской эндокринологии ФГБУ ЭНЦ в период с 1978 по 2016 г. по поводу СТГ-дефицита и достигших КР. Был использован сплошной метод формирования выборки. Период наблюдения составил от 4 до 38 лет. Информация о 23 (16%) пациентах собрана ретроспективно, о 118 (84%) — проспективно. К началу терапии рГР 114 (81%) человек находились в допубертатном возрасте и 27 (19%) — в пубертатном.

Критерии соответствия

Критерии включения в исследуемую группу:

— диагноз СТГ-дефицита, подтвержденный пробой с клофелином (максимальный уровень СТГ <10 нг/мл);

— достижение КР;

— регулярная терапия рГР в течение минимум 3 лет до достижения КР;

— клинико-лабораторная компенсация вторичных гормональных недостаточностей в течение периода наблюдения.

Критерии исключения:

— наличие объемных образований гипофиза;

— маленькие для гестационного возраста рост и масса тела при рождении (<2 SD для срока гестации) [15];

— синдромальная задержка роста;

— костно-хрящевые дисплазии;

— наличие тяжелых сопутствующих заболеваний;

— декомпенсация вторичных недостаточностей в течение периода наблюдения.

Продолжительность исследования

Период наблюдения за пациентами составил от 4 до 15 лет. Частота клинико-лабораторного обследования — 1 раз в год.

Описание медицинского вмешательства

После подтверждения диагноза СТГ-дефицита все пациенты получали заместительную терапию рГР в дозе 0,033 мг/кг (ежедневные подкожные инъекции в вечернее время).

Методы регистрации исходов

Клинико-лабораторные показатели:

  1. Пол (м/ж).
  2. Хронологический возраст на момент диагностики СТГ-дефицита (с точностью до дня).
  3. Пубертатный статус (допубертатный/пубертатный) по классификации Таннера. Объем тестикул при осмотре мальчиков исследовался с помощью орхидометра Прадера.
  4. Форма заболевания (ИДГР/МДГА) устанавливалась на основании лабораторных исследований.
  5. Максимальный уровень СТГ в ходе теста с клофелином (нг/мл). Исследование крови проводилось в утренние часы натощак. Клофелин (клонидин) вводился однократно внутрь в дозе 0,15 мкг/м2 поверхности тела. Уровень СТГ исследовался до и через 30, 60, 90, 120 мин после введения препарата.
  6. Регулярность терапии рГР (да/нет) оценивалась путем опроса пациентов. Перерыв в лечении препаратами рГР суммарно не более 1 мес в год оценивался как регулярная терапия, суммарно более 1 месяца — как нерегулярная.

Ауксологические показатели:

  1. SDS роста при рождении — рассчитывался по формуле: SDS роста=(х–Х)/SD, где х — рост ребенка, Х — средний рост для данного хронологического возраста и пола, SD – стандартное отклонение роста для данного хронологического возраста и пола (при рождении для мальчиков российской популяции SD=2,02 см, Х=54,79 см, для девочек SD=2,02 см, Х=53,71 см [12]).
  2. SDS роста для хронологического возраста и пола на момент диагностики СТГ-дефицита — рассчитывался для оценки степени отклонения роста пациента от среднего в британской популяции с помощью приложения Auxology [17] по результатам измерения длины тела пациента при помощи механического ростомера с точностью до 0,1 см (норма ± 2SD).
  3. SDS генетически прогнозируемого роста (SDS ГПР) — рассчитывался с использованием данных о росте родителей пациента по формуле Таннера: ГПР для мальчиков = [(рост отца + рост матери + 13 см)/2] ±7 см; ГПР для девочек = [рост отца + рост матери – 13 см)/2] ± 7 см [7].
  4. Костный возраст пациента (КВ) на момент диагностики СТГ-дефицита (с точностью до месяца). Оценка степени дифференцировки скелета («костного возраста») проводилась методом Greulich & Pyle [16] по рентгенограмме кистей рук с лучезапястными суставами.
  5. Отношение «костный возраст/хронологический возраст» (КВ/ХВ) на момент диагностики СТГ-дефицита — рассчитывалось арифметически.
  6. КР (в см) фиксировался при достижении пациентом хронологического возраста больше 16 лет и/или костного возраста более 16 лет у мальчиков и более 14 лет у девочек, а также при снижении скорости роста до менее 2 см в год на фоне терапии рГР.
  7. SDS КР рассчитывали с помощью приложения Auxology (для британской популяции) [17].

Этическая экспертиза

Исследование одобрено локальным этическим комитетом ФГБУ «Эндокринологический научный центр» (протокол ЛЭК №12 от 22.10.14).

Статистический анализ

Анализ данных проводился с использованием пакета прикладных программ Statistica v.13 («StatSoft Inc.», США). В качестве параметров распределений количественных признаков использованы медиана (Me), 1-й и 3-й квартили (Q1 и Q3), в случае их нормальных распределений — средние и среднеквадратические отклонения; для бинарных признаков — абсолютные и относительные частоты. Оценка соответствия распределения нормальному осуществлялась с применением критерия Шапиро—Уилка. Для сравнения независимых групп по количественному признаку использовался метод Манна—Уитни. Для анализа связи количественных признаков применялся непараметрический корреляционный анализ по Спирмену. Для нивелирования проблемы множественных сравнений применялась поправка Бонферрони.

Применялись два метода математического моделирования — пошаговая множественная линейная регрессия (МЛР) и искусственные нейронные сети (ИНС). В обоих случаях откликом служили непрерывные признаки — КР и SDS КР. МЛР допускает в качестве предикторов только непрерывные признаки. Пошаговый вариант этого анализа позволяет оптимизировать число и набор предикторов в окончательной модели таким образом, чтобы достичь максимальной эффективности модели при минимальном наборе предикторов. ИНС — класс методов математического моделирования, в основе которого лежит бионический принцип сети узлов (нейроны) с множественными входными сигналами (дендритами) и одним выходным сигналом (аксон). На входной слой сети подаются данные, выходным сигналом является моделируемый отклик (рис. 1). Преимуществом данного класса методов является возможность использовать сочетание любых признаков как в качестве предикторов, так и в качестве откликов моделей, а также высокая эффективность результирующих моделей. Ранее было показано, что ИНС являются ценным инструментом создания математических моделей прогнозирования при решении медицинских задач, особенностью которых служит малое число наблюдений и необходимость совместного использования качественных и количественных предикторов [7]. В качестве топологий сети тестировались многослойные (3- и 4-слойные) персептроны (МСП), радиальные базисные функции, байесовские и линейные нейронные сети. Оба метода моделирования исключают наблюдения с пропусками из анализа, поэтому количество наблюдений при различных наборах предикторов могло различаться. Для оценки качества модели использовались среднеквадратическая ошибка модели (root mean square error — RMSE), коэффициент детерминации (R2), коэффициенты корреляции наблюдаемых и прогнозируемых величин.

Рис. 1. Схематическое изображение 3-слойного персептрона — одной из наиболее популярных топологий искусственных нейронных сетей. Между входным слоем (из 10 нейронов — слева) и выходным слоем (из 1 искусственного нейрона — справа) расположен слой из 6 скрытых искусственных нейронов.

Результаты

Объекты (участники) исследования

Характеристика выборки пациентов, достигших КР, представлена в табл. 1. Как можно видеть, данные содержат до 9% пропусков (максимальное количество пропусков — по максимальному стимулированному уровню СТГ). Поскольку многомерный анализ исключает наблюдения с пропусками (либо требует заполнения пропусков, что обычно приводит к смещенным оценкам и завышению качества получаемых моделей), нами были вынужденно удалены все наблюдения с пропусками, после чего осталось 122 наблюдения для проведения регрессионного анализа (в таком анализе участвуют только количественные признаки) и 121 — для применения методов ИНС.

Таблица 1. Характеристика группы пациентов с КР (n=141)

Параметр

Абс.

Описательная статистика, абс. частоты и % или Me (Q1; Q3)

Пол (м/ж), абс. (%)

141

90/51 (64/36)

Возраст на момент установки СТГ-дефицита, годы

140

9,62 (6,82; 12,91)

Пубертатный статус (допуб/пуб.), абс. (%)

140

114/26 (81/19)

Диагноз (ИДГР/МДГА), абс. (%)

141

35/106 (25/75)

Регулярность терапии (да/нет), абс. (%)

141

99/42 (70/30)

SDS роста при рождении

132

0,56 (–0,38; 1,24)

SDS роста на момент диагностики СТГ-дефицита

139

–3,39 (–4,16; –2,83)

КВ на момент диагностики СТГ-дефицита, годы

133

6,0 (3,0; 8,0)

КВ/ХВ

133

0,58 (0,48; 0,68)

КР, см:

141

166,0 (160,0; 173,1)

у девочек

 

158,0 (153,8; 162,3)

у мальчиков

 

171,0 (165,0; 176,0)

Доля пациентов, не достигших социально приемлемой нормы [8], абс. (%)

141

32 (23)

девочек (<154 см)

13 (25)

мальчиков (<165 см)

19 (21)

SDS КР

141

–0,67 (–1,37; 0,15), min= –4,00, max=1,80

SDS ГПР

131

–0,03 (–0,56; 0,42)

Максимальный стимулированный уровень СТГ, нг/мл

129

0,66 (0,20; 2,20)

Основные результаты исследования

Результаты анализа статистической связи отдельных предикторов с КР и SDS КР для 121 полного наблюдения приведены в табл. 2. Наборы признаков, статистически значимо связанные с КР и SDS КР, несколько различаются. Наличие или отсутствие таких связей может использоваться в качестве аргументов при формировании наборов предикторов для построения многофакторных моделей, однако такой подход не является ни единственным, ни заведомо эффективным. Другими возможными способами снижения размерности признакового пространства являются пошаговые процедуры в МЛР, генетический алгоритм, эвристические аргументы (например, сложность или трудоемкость оценки некоторых предикторов).

Таблица 2. Характеристика одномерных статистических связей потенциальных предикторов с моделируемыми откликами (n=121)

Потенциальный предиктор

Отклик: КР

Отклик: SDS КР

Мe (Q1; Q3)

p, тест Манна–Уитни

Мe (Q1; Q3)

p, тест Манна—Уитни

Пол:

муж. (n=78)

171 (166; 176)

<0,001*

–0,45 (–1,31; 0,27)

0,649

жен. (n=43)

159 (154; 163)

–0,67 (–1,34; 0,16)

Диагноз:

МДГА (n=89)

167 (160; 175)

0,323

–0,34 (–1,26; 0,48)

0,006

ИДГР (n=32)

165 (159; 171)

–0,93 (–1,52; –0,49)

Пубертатный статус:

пуб. (n=24)

162 (158; 170)

0,069

–0,71 (–1,39; –0,27)

0,344

допуб. (n=97)

168 (161; 175)

–0,45 (–1,31; 0,23)

Регулярность терапии:

да (n=88)

169 (164; 176)

<0,001

–0,29 (–1,08; 0,42)

<0,001

нет (n=33)

158 (154; 161)

–1,22 (–1,85; –0,67)

Потенциальные предикторы

Отклик: КР

Отклик: SDS КР

коэффициенты корреляции

p, тест

Спирмена

коэффициенты

корреляции

p, тест

Спирмена

SDS роста при рождении

0,105

0,249

0,078

0,392

SDS целевого роста

0,340

<0,001

0,365

<0,001

Возраст при постановке диагноза, годы

–0,208

0,022

–0,216

0,017

SDS роста при постановке диагноза

0,457

<0,001

0,295

0,001

Костный возраст, годы

–0,202

0,026

–0,262

0,004

Отношение костного возраста к хронологическому

–0,105

0,251

–0,185

0,042

Максимальный СТГ, нг/мл

0,127

0,164

–0,016

0,866

           

Поскольку МЛР допускает применение в качестве предикторов лишь количественных признаков, были использованы следующие показатели: SDS роста при рождении, SDS ГПР, хронологический возраст на момент диагностики СТГ-дефицита, SDS роста на момент диагностики СТГ-дефицита, КВ на момент установки СТГ-дефицита, КВ/ХВ, максимальный стимулированный уровень СТГ. В выборке было 122 наблюдения с отсутствием пропусков в этих признаках; соответственно, модели строились на этих 122 наблюдениях. Разработаны две модели — для прогнозирования КР и SDS КР, их характеристики приведены в табл. 3.

Таблица 3. Характеристики множественных линейных регрессионных моделей (n=122)

Характеристика

Моделирование КР

Моделирование SDS КР

Коэффициенты уравнений регрессии для предикторов:

SDS роста при рождении

SDS генетически прогнозируемого роста

2,112

0,313

Возраст на момент установки СТГ-дефицита, годы

SDS роста на момент установки СТГ-дефицита

3,160

0,268

Костный возраст на момент установки СТГ-дефицита, годы

–0,493

–0,082

Коэффициент КВ/ХВ

Максимальный стимулированный уровень СТГ, нг/мл

–0,077

Константа

180,5

0,989

Стандартная ошибка

7,79 см

–0,93 SDS

p (уровень статистической значимости модели)

<0,001

0,006

Доля объяснимой дисперсии R2, %

27

26

В модели МЛР, рассчитывающей КДР, статистически значимыми параметрами были SDS роста пациента при диагностике СТГ-дефицита, КВ на момент диагностики СТГ-дефицита и SDS генетически прогнозируемого роста. В модель, прогнозирующую SDS КР, дополнительно включен максимальный стимулированный уровень СТГ. Остальные 3 предиктора (SDS роста при рождении, возраст на момент диагностики СТГ-дефицита, коэффициент КВ/ХВ) оказались статистически незначимыми. Несмотря на статистическую значимость обеих моделей, они характеризуются низкой долей объяснимой дисперсии.

Модели прогнозирования КР и SDS КР строились также с применением метода ИНС. Однако теперь в качестве предикторов использовались и 4 бинарных признака. Применялись два набора признаков:

— все признаки (полный набор предикторов);

— признаки, имеющие статистически значимую одномерную связь с прогнозируемыми показателями (сокращенный набор предикторов) (см. табл. 2).

Исходная выборка была случайным образом разделена на обучающую и контрольную в соотношении 7:3 (85 и 36 наблюдений соответственно). Тестовая выборка вынужденно не создавалась в связи с крайне малым количеством полных наблюдений, а также с учетом того, что ошибка на ней обычно приближенно равна таковой на контрольной.

Наилучшие результаты при прогнозировании КР были получены на моделях многослойного персептрона при применении полного набора предикторов (см. табл. 3). Общая ошибка на лучшей ИНС (МСП 11:15-4—1:1) составила 4,417 см, а доля объяснимой дисперсии — 76%, что представляется неплохим результатом, однако и в этом случае проявился типичный недостаток ИНС — их переобученность, особенно при малом числе наблюдений: ошибки на контрольных выборках выше ошибок на обучающих, либо не согласуются друг с другом.

На рис. 2 приведена диаграмма рассеяния фактических и прогнозируемых значений КР для лучшей полученной ИНС. Видно, что эти значения расположены вблизи линии их совпадения, корреляция этих показателей является сильной (r=0,873; p<0,001, критерий Спирмена). Наблюдения, расположенные над прямой, имеют переоценку КР, под прямой — недооценку. Гистограмма распределения остатков модели (разностей между прогнозируемыми и фактическими величинами) приведена на рис. 3. Остатки со средним и стандартным отклонением (–0,08)±4,44 нормально распределены (p=0,916; критерий Шапиро—Уилка).

Рис. 2. Диаграмма рассеяния фактических значений конечного достигнутого роста (КДР) (ось абсцисс) и прогнозируемых с помощью трехслойного персептрона 11:15—4-1:1 значений (ось ординат).

Рис. 3. Гистограмма распределения остатков модели ИНС для прогнозирования КР (разностей между фактическими и прогнозируемыми величинами).

Модели прогнозирования SDS КР с помощью ИНС оказались менее эффективными. Лучшая из них, построенная по полному набору предикторов, описана в табл. 4. Доля объяснимой дисперсии была высокой на обучающей выборке, но крайне низкой — на контрольной. Модели, построенные по сокращенному набору предикторов, также оказались неэффективными, что подтверждает нецелесообразность использования критерия наличия одномерных связей предиктора и отклика при формировании набора потенциальных предикторов моделей.

Таблица 4. Характеристики разработанных моделей искусственных нейронных сетей (n=121)

Характеристика

Моделирование КР

Моделирование SDS КР

полный набор предикторов

сокращенный набор предикторов

полный набор предикторов

сокращенный набор предикторов

Количество предикторов

11

4

11

4

Архитектура сети

МСП (3)

11:15—4-1:1

МСП (3)

4:4—4-1:1

МСП (3)

11:15—9-1:1

МСП (3)

4:4—10-1:1

Среднеквадратическая ошибка (RMSE):

все наблюдения

4,417 см

5,546 см

0,601 SDS

0,823 SDS

обучающая выборка

4,129 см

5,495 см

0,404 SDS

0,840 SDS

контрольная выборка

5,033 см

5,664 см

0,911 SDS

0,784 SDS

Доля объяснимой дисперсии (R2):

все наблюдения, %

75,9

62,0

42,4

43,7

обучающая выборка, %

77,6

60,1

86,0

35,7

контрольная выборка, %

72,4

65,5

28,0

59,8

Коэффициент корреляции Спирмена:

все наблюдения

0,873

0,784

0,824

0,577

обучающая выборка

0,882

0,766

0,933

0,509

контрольная выборка

0,823

0,830

0,545

0,709

Уравнения моделей ИНС являются громоздкими, поэтому в статье не приводятся, однако лучшая из моделей (МСП (3) 11:15—4—1:1) была программно реализована: разработан свободно доступный в Интернет программный калькулятор (http://alfa-endo.ru/page/programmnii_kalkulyator_rascheta_konechnogo_dostig).

Обсуждение

Возможность прогнозировать КР — важный инструмент клинициста. С его помощью можно осуществлять персонализированный подход к терапии рГР у пациентов с СТГ-дефицитом: подбирать дозу рГР, контролировать соблюдение регулярности введения и дозирования препарата, оценивать эффективность терапии, наглядно демонстрировать факторы, от которых зависит КР.

Первые прогностические модели использовались для прогнозирования скорости роста в первый год терапии рГР. Считается, что именно этот год максимально отражает восприимчивость организма пациента к лечению. Однако после того как появилась возможность анализировать поколение детей, регулярно и длительно получавших терапию рГР и достигших роста близкого к конечному (до закрытия зон роста), стали появляться модели, прогнозирующие КР и его SDS [9, 10]. Стоит отметить, что КР как у здоровых людей, так и у пациентов с СТГ-дефицитом различается в зависимости от популяции [2, 4]. Для российской популяции социально приемлемым ростом считается рост выше 154 см у женщин и выше 165 см у мужчин [8]. Значительным недостатком большинства известных в настоящее время математических моделей прогнозирования КР и его SDS являются низкая доля объяснимой дисперсии и достаточно высокая погрешность вычисления. К тому же в этих моделях используется информация, доступная на момент начала терапии рГР, а для расчета КР детей в зависимости от пубертатного статуса необходимо использовать две разные формулы [9]. К недостаткам моделей МЛР (такими являются 5 из 6 известных моделей) относится невозможность использования бинарных предикторов, таких как пол, вид СТГ-дефицита (ИДГР/МДГА), регулярность терапии рГР, пубертатный статус пациентов. Модели ИНС лишены этого недостатка. Первая ИНС для прогнозирования эффективности терапии рГР у пациентов с ИДГР была разработана J. Smyczynska и соавт. [10] в 2015 г. Она обладает достаточно высокой прогностической точностью при прогнозировании КР (СКО=3,8 см, R2=77,7%) и более низкой — при прогнозировании SDS КР (СКО=0,59 SDS, R2=43,0%). Наши модели характеризуются сходной прогностической точностью (СКО=4,4 см, R2=75,9% при прогнозировании КР; СКО=0,60 SDS, R2=42,4% при прогнозировании SDS КР), однако учитывают более широкий набор предикторов — форму заболевания, пубертатный статус, регулярность терапии. Учитывая, что СТГ-дефицит является хроническим заболеванием и требует длительного лечения, а регулярность этого лечения влияет на его успешность [13, 14], последний предиктор представляется очень важным.

Все предикторы, используемые в наших моделях, доступны для оценки практическим врачом-эндокринологом и не требуют малодоступных исследований и частого динамического наблюдения пациента. Для удобства применения разработанных моделей на практике создан доступный в Интернет программный продукт.

Ограничения исследования

Основным ограничением исследования является недостаточное число наблюдений. С этим связана малая эффективность и переобученность полученных моделей, хотя методы ИНС позволили получить модель приемлемого качества. Попытки снизить число предикторов не привели к улучшению результатов.

При проведении многомерного анализа исключаются наблюдения с пропусками, поэтому мы включали в рассмотрение только такие показатели, по которым доля пропусков не превышала 10%. В частности, поэтому в число рассматриваемых признаков не были включены результаты анализа мутаций в генах, ответственных за эмбриогенез гипофиза (ARNT2, GH1, GHRH, GHRHR, GHSR, GLI2, HESX1, LHX3, LHX4, OTX2, PAX6, POU1F1, PROP1, SHH, SOX2, SOX3).

Заключение

ИНС является ценным инструментом для создания эффективных моделей в медицине, в частности, при прогнозировании эффективности терапии рГР у пациентов с СТГ-дефицитом. В качестве исходных данных для разработанных нами моделей используются доступные для любого врача-клинициста показатели. Лучшая из разработанных моделей продемонстрировала достаточно высокую прогностическую точность при прогнозировании КР (ошибка составила 4,417 см, а доля объяснимой дисперсии 75,9%). Модель для прогнозирования SDS КР была менее качественной (ошибка составила 0,601 SDS, а доля объяснимой дисперсии 42,4%). В перспективе использование большего объема данных для моделирования позволит улучшить качество прогнозирования эффективности терапии рГР.

Дополнительная информация

Конфликт интересов. Авторы декларируют отсутствие явных и потенциальных конфликтов интересов, связанных с выполнением данного исследования.

Благодарности. Коллектив авторов благодарит компанию «StatSoft Russia» за помощь в генерации программного кода и Фонд поддержки развития филантропии «КАФ» за программную реализацию калькулятора и его размещение в Интернет.

Список литературы

1. Darendeliler F, Lindberg A, Wilton P. Response to growth hormone treatment in isolated growth hormone deficiency versus multiple pituitary hormone deficiency. Horm Res Paediatr. 2011;76 Suppl 1:42-46. doi: 10.1159/000329161

2. Blethen SL, Baptista J, Kuntze J, et al. Adult height in growth hormone (GH)-deficient children treated with biosynthetic GH. The Genentech Growth Study Group. J Clin Endocrinol Metab. 1997;82(2):418-420. doi: 10.1210/jcem.82.2.3734

3. Rachmiel M, Rota V, Atenafu E, et al. Final height in children with idiopathic growth hormone deficiency treated with a fixed dose of recombinant growth hormone. Horm Res. 2007;68(5):236-243. doi: 10.1159/000101427

4. Thomas M, Massa G, Bourguignon JP, et al. Final height in children with idiopathic growth hormone deficiency treated with recombinant human growth hormone: the Belgian experience. Horm Res Paediatr. 2001;55(2):88-94. doi: 10.1159/000049976

5. Mortensen HB, Main K, Michaelsen KF, et al. Predicting and monitoring of growth in children with short stature during the first year of growth hormone treatment. Acta Paediatr. 1991;80(12):1150-1157. doi: 10.1111/j.1651-2227.1991.tb11803.x

6. Wikland KA, Kriström B, Rosberg S, et al. Validated multivariate models predicting the growth response to GH treatment in individual short children with a broad range in GH Secretion Capacities. Pediatr Res. 2000;48(4):475-484. doi: 10.1203/00006450-200010000-00010

7. Schonau E, Westermann F, Rauch F, et al. A new and accurate prediction model for growth response to growth hormone treatment in children with growth hormone deficiency. Eur J Endocrinol. 2001;144(1):13-20. doi: 10.1530/eje.0.1440013

8. Ranke MB, Martin DD, Lindberg A, on behalf of the Kigs International Board o. prediction model of total pubertal growth in idiopathic growth hormone deficiency: analysis of data from KIGS. Horm Res Paediatr. 2004;60(1):58-59. doi: 10.1159/000071227

9. de Ridder MAJ, Stijnen T, Hokken-Koelega ACS. Prediction of adult height in growth-hormone-treated children with growth hormone deficiency. J Clin Endocr Metab. 2007;92(3):925-931. doi: 10.1210/jc.2006-1259

10. Smyczynska J, Hilczer M, Smyczynska U, et al. Neural network models — a novel tool for predicting the efficacy of growth hormone (GH) therapy in children with short stature. Neuro Endocrinol Lett. 2015;36(4):348-353.

11. Реброва О.Ю. Применение методов интеллектуального анализа данных для решения задачи медицинской диагностики. // Новости искусственного интеллекта. — 2004. — №3. — С. 76-80. [Rebrova OYu. Primenenie metodov intellektual’nogo analiza dannykh dlya resheniya zadachi meditsinskoy diagnostiki. Novosti iskusstvennogo intellekta. 2004;(3):76-80.(In Russ.)].

12. Нагаева Е.В. Федеральные клинические рекомендации по диагностике и лечению гипопитуитаризма у детей и подростков // Проблемы Эндокринологии. — 2013. — Т. 59. — №6. — C. 27—43. [Nagaeva EV. Russian association of endocrinologists national practice guidelines(clinical signs, diagnosis, differential diagnosis, treatment). Hypopituitarism in children and adolescents. Problems of Endocrinology. 2013;59(6):27-43. (In Russ.)].doi: 10.14341/probl201359627-43

13. Blethen SL, Baptista J, Kuntze J, et al. Adult height in growth hormone (GH)-deficient children treated with biosynthetic GH. The Genentech Growth Study Group. J Clin Endocrinol Metab. 1997;82(2):418-420. doi: 10.1210/jcem.82.2.3734

14. Rohayem J, Drechsel H, Tittel B, et al. Long-term outcomes, genetics, and pituitary morphology in patients with isolated growth hormone deficiency and multiple pituitary hormone deficiencies: a single-centre experience of four decades of growth hormone replacement. Horm Res Paediatr. 2016;86(2):106-116. doi: 10.1159/000448098

15. WHO Child Growth Standards. Acta Paediatrica. 2006;95(suppl 1):5-101. Available from: http://www.who.int/childgrowth/standards/ru

16. Greulich WW, Pyle SI. Radiographic atlas of skeletal development of the hand and wrist. The American Journal of the Medical Sciences. 1959;238(3):393.

17. Kromeyer-Hauschild K, Jaeger U. Growth studies in Jena, Germany: changes in body size and subcutaneous fat distribution between 1975 and 1995. Am J Hum Biol. 1998;10(5):579-587. doi: 10.1002/(sici)1520-6300(1998)10:5<579::aid-ajhb4>3.0.co;2-9


Об авторах

Анна Евгеньевна Гаврилова

ФГБУ «Национальный медицинский исследовательский центр эндокринологии» Минздрава России


Россия

клин.аспирант



Елена Витальевна Нагаева

ФГБУ «Национальный медицинский исследовательский центр эндокринологии» Минздрава России


Россия

к.м.н., ведущий научный сотрудник 



Ольга Юрьевна Реброва

ФГБУ «Национальный медицинский исследовательский центр эндокринологии» Минздрава России


Россия

д.м.н., проф. каф. эндокринологии



Татьяна Юрьевна Ширяева

ФГБУ «Национальный медицинский исследовательский центр эндокринологии» Минздрава России


Россия

к.м.н., зав. отд. тиреодологии репродуктивного и соматического развития



Дополнительные файлы

1. Математические модели для прогнозирования конечного достигнутого роста и его коэффициента стандартного отклонения у детей с дефицитом гормона роста российской популяции.
Тема
Тип Прочее
Скачать (161KB)    
Метаданные
2. Рис. 1. Схематическое изображение 3-слойного персептрона — одной из наиболее популярных топологий искусственных нейронных сетей. Между входным слоем (из 10 нейронов — слева) и выходным слоем (из 1 искусственного нейрона — справа) расположен слой из 6 скрытых искусственных нейронов.
Тема
Тип Исследовательские инструменты
Посмотреть (41KB)    
Метаданные
3. Рис. 2. Диаграмма рассеяния фактических значений конечного достигнутого роста (КДР) (ось абсцисс) и прогнозируемых с помощью трехслойного персептрона 11:15—4-1:1 значений (ось ординат).
Тема
Тип Исследовательские инструменты
Посмотреть (44KB)    
Метаданные
4. Рис. 3. Гистограмма распределения остатков модели ИНС для прогнозирования КР (разностей между фактическими и прогнозируемыми величинами).
Тема
Тип Исследовательские инструменты
Посмотреть (29KB)    
Метаданные

Рецензия

Для цитирования:


Гаврилова А.Е., Нагаева Е.В., Реброва О.Ю., Ширяева Т.Ю. Математические модели прогнозирования конечного роста и его коэффициента стандартного отклонения у детей с дефицитом гормона роста в российской популяции. Проблемы Эндокринологии. 2017;63(5):282-290.

For citation:


Gavrilova A.E., Nagaeva E.V., Rebrova O.Yu., Shiryaeva T.Yu. Mathematical models for predicting final achieved height and its standard deviation rate in children with growth hormone deficiency in the Russian population. Problems of Endocrinology. 2017;63(5):282-290.

Просмотров: 6758


Creative Commons License
Контент доступен под лицензией Creative Commons Attribution 4.0 License.


ISSN 0375-9660 (Print)
ISSN 2308-1430 (Online)