«В новый век – с доказательной биомедициной!»Возможности и практика применения статистических методов для экспериментальной би

Практика студентов кафедры «Биофизика» МГУ им.М.В.Ломоносова на Белом море(1967г.)

Гален (Galen)
«Эмпирики говорят, что нечто не может быть ни принято, ни расценено как истинное, если оно замечено однократно или если оно было замечено только несколько раз. Они полагают, что нечто может быть принято и считаться верным, только если это было замечено очень много раз и каждый раз в том же самом виде»
On Medical Experience. 11-столетие н.э.

По вине медиков в год в России гибнут 50 тыс. чел.(Лига защиты пациентов); в ДТП–30 тыс.чел.(Росстат)
ЗАДАЧИ ПРЕЗЕНТАЦИИ:
Анализ публикаций с применением статистических методов.
Практика наиболее используемых статистических методов.
Ошибки при интерпретации результатов исследований.
Научная значимость некоторых малоиспользуемых статистических методов.

А.Н.Мамаев д.м.н.,зав.отделением
Алтайского филиала ГНЦ РАМН (2011г.)
«Неправильное применение методов математической статистики в медицине на сегодняшний день в России приобретает характер эпидемии, поскольку ошибки исследователей при анализе медицинских данных можно встретить повсеместно и число их неуклонно растет»

Леонов В.П. Анализ публикаций по экспериментальной биомедицине (www.biometrica.tomsk.ru)

Результаты анализа публикаций
Использование стат.методов: 82% в R-статьях и 87% в E-статьях.
Использование описательной статистики и t-критерия (ч.50-85%) среди 53-х методов и критериев!!!
В 52% R-статей нет описания использованных методов и точного значения уровня достоверности «р».
Отсутствие проверки условий применимости критериев и методов.
Неверная интерпретация результатов исследования

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

Статистические методы
Описательная статистика
Анализ взаимосвязей (корреляционный анализ)
Регрессионный анализ
Анализ временных рядов
Анализ выживаемости
Принцип Парето (ABC анализ)

Значения медицинского показателя до и после лечения

Описательная статистика (R-ч. 82%)
Mean – среднее арифметическое
Ассиметрия (скошенность), эксцесс (пикообразность) распределения
Медиана разбивает выборку на две равные части (50-й процентиль); 25-й и 75-й процентили (квартили)
Мода – наиболее часто встречающееся значение в выборке
+SE – стандартная ошибка среднего
Дисперсия S2– среднее значение квадратов отклонений значений показателя от среднего значения
СО –стандартное отклонение S
1,96*SE (ДИ) – доверительный интервал показывает, что с вероятностью 0,95 показатель X будет находиться в интервале значений X+1,96*SE.
Коэффициент вариации (%) – отношение стандартного отклонения к среднему значению CV=СО/X*100% -воспроизводимость метода исследования (сравнение разных методик).CV>10% -плохая воспроизводимость метода иссл-я
Размах – самое большое и самое малое значение

Показатели описательной статистики (Excel)

Ящичковая диаграмма Тьюки
(«ящик с усами»-«вне отечественной традиции»)

Доктор L.Brown (1941г.)
«Затруднения в медицине исходят не от статистического метода, а от медицинских работников (и работниц), которые не знают, как его использовать….Не поймите меня неправильно. Это не аргумент в пользу сухих статистических статей, которые все мы не читаем. Но если я смогу заставить вас увидеть, насколько важно для нас прекратить использовать любимую фразу «мой личный опыт» кроме тех случаев, когда у нас есть достаточно данных, чтобы подтвердить ее, я буду считать, что достиг того, на что я надеялся»

Типичные ошибки интерпретации результатов исследований
Уровень значимости “р” t-критерия равен вероятности ошибочно отвергнуть гипотезу о равенстве средних двух выборок. Приняли, что р<0,05!
Неверно: «Уровень значимости р<0,05»
Правильно: р=0,04 (или р=0,06). «Наблюдаемое различие в эксперименте статистически значимо на уровне р=0,04 или р=0,06».
Неверно: Среднее = (72+6)ед.
Правильно: Среднее (СО)= 72ед. (8ед.)
Среднее (при 95%ДИ)=72ед.(от 66 до 78)ед.

Проверка гипотез о равенстве средних t-критерием (R-ч. 58%; Е-ч.8%)
Требования для применения t-критерия.
Нормальное распределение показателя
Равенство дисперсий сравниваемых совокупностей
Практика: В выборке из 200-х публикаций в 50-70% реальных данных нормальное распределение не выполняется! (Леонов В.П.)
В 2500 биомед. массивов данных только в 8% было нормальное распределение (г.Тарту)
Непараметрические методы (R-ч. 2%),
т.е. в 48-68% публикациях неверные выводы!

Нормальное распределение
- фундаментальный закон природы
Мировая практика: нормальное распределение 20-25% !!!

Проверка на нормальность(R-ч.0,06%)
Визуальный метод
График вероятности
Критерий Колмогорова-Смирнова.
Если K-S d>0,0895, распределение не соответствует нормальному на уровне значимости 0,05.
Критерий Шапиро-Уилка
Совпадение значений моды и медианы
Равенство нулю показателей асимметрии и эксцесса (таблицы критических значений для заданной выборки)

Проверка на нормальность с помощью графика вероятности

Пример невыполнения нормального распределения (до лечения)

Пример невыполнения нормального распределения (после лечения)

Пример нормального распределения (Результаты голосования в г.Москве 4 марта 2012г.)

Статистическая «Развлекуха» или «чуровское распределение»

Графики медицинского показателя

Гистограмма медицинского показателя до и после лечения (пример отсутствия нормального распределения)

Выявление различий показателей в группе до и после лечения (связанные выборки)
Пример 1: Значение медианы в группе до лечения = 36,5ед.; а после лечения = 34,2ед.
Параметрический способ: t-тест: p=0,108
Различие незначимо. Эффекта лечения нет.
Непараметрический метод сравнения 2-х связанных выборок с помощью критерия Вилкоксона (R-ч. 0,7%): W=57.
Критическое значение (для n=20) Wкр.=53.
Так как W > Wкр., принимается гипотеза о наличии различий между значениями показателя в группе
(до лечения) и (после лечения), т.е. выявлено влияние выбранного метода на результат лечения!

Выявление эффекта от лечения в 2-х группах (несвязанные выборки)
Пример 2:
Значение медианы в 1-ой группе (лек-во 1) - 29,6ед.;
во 2-ой группе (лек-во 2) - 35,9ед.
Параметрический способ: t-тест (p=0,13).
Различие незначимо. Разницы в эффективности лекарств нет.
Непараметрический критерий Манна-Уитни для сравнения 2-х несвязанных выборок (R-ч. 2%): U=66.
Критическое значение Uкр.=70 (для n1=17 и n2=13).
Так как U < Uкр., принимается гипотеза о наличии различий между двумя выборками, т.е. выявлена разница в эффективности от использования двух лекарств!

Вместо заключения о медицинской описательной статистике
Все результаты и выводы в статьях, отчетах и диссертациях, полученные на основе данных описательной статистики с использованием параметрических методов (среднее, СО,ДИ и t-тесты) в группах больных не верны ввиду отсутствия нормального распределения данных для многих медицинских показателей!!!

Использование непараметрических критериев

2. Анализ связей и корреляций (R-ч. 8%)
K.Peаrson (1901г.) - обработка биомедицинских данных
1. Взаимосвязь качественных переменных –критерий Пирсона хи-квадрат (R-ч. 7%)
2. Взаимосвязь количественных переменных –критерий корреляции Пирсона “r” (R-ч. 8%)
Принято считать, «что при r < 0,25 корреляция слабая, при 0,25 < r < 0,75 умеренная, при r > 0,75 – сильная».
Коэффициент детерминации: R=r2 –доля изменения Y, вызванная изменением X.
Пример: Если r=0,75, то R=0,56. Cильная связь?

Критерий Пирсона «Хи-квадрат».
Таблица сопряженности 2х2 (R-ч.7%; E-ч.19%)

Пример «умеренной» связи показателей

Корреляционная матрица биохимических показателей до лечения

Корреляционная матрица биохимических показателей после лечения

Краткие научные выводы:
Измерена теснота связи между различными медицинскими показателями до и после лечения (слабая, умеренная, сильная).
Выявлена тесная и умеренная связь между отдельными показателями.
Обнаружено изменение тесноты связи (критерий r) и ковариации (+r) между показателями до и после лечения.
Вопрос для науки: В чем причина произошедших изменений?

3.Регрессионный анализ (R-ч. 4%)

Парная регрессионная зависимость между медицинскими показателями

Парная регрессия после верификации (искл. одной точки)

Стоимость строительства кв.м жилья в России

Результаты прогнозирования

4. Анализ временных рядов (R-ч.0,1%;
E-ч.3%) (наличие сезонной компоненты)

Прогнозные значения показателя с доверительным интервалом

Потребление пива в России (данные Росстата)

Динамика заболевания алкоголизмом и наркоманией в России (данные Росстата)
1990
2000
2010

5. Анализ выживаемости - оценка эффективности лечения (R-ч.0,3%;Е-ч.9%)
Оценка вероятности пережить пациенту произвольный интервал времени, измеренный от даты стартового до мониторируемого события.
Мониторируемое событие – рецидив, ремиссия заболевания, летальный исход и т.д.
Стартовое событие – начало лечения, оперативное вмешательство, поступление в реанимацию (в клинику) и т.д.

Градация участников исследования

Таблица времен жизни

Функция риска
Число смертей/число доживших

Kaplan-Meier product-limit method
Кривая Каплана-Мейера - частота выживания больных раком в двух группах (R-ч.0,3%; E-ч.9% )

Доказательный вывод: Хирургическое лечение более эффективное для продления жизни по сравнению с медикоментозным лечением
Медиана времени выживания составила 6,3 года для группы хирургического лечения и 3,8 года для группы медикоментозного лечения.
Лог-ранговый критерий выявил статистически значимое различие по частоте выживания с течением времени (р<0,001).

6. Правило Парето (правило «20 на 80») (ч.0,01%) «За большинство возможных результатов отвечает относительно небольшое число причин».
Научно-методический центр по обеспечению качественной медицинской помощи ЦНИИ ОИЗ
Доллары США Процент

Диаграмма Парето

Jean-Etienne D.Esquirol (французский психиатр, 1838г.)
«Думали ли они, что науки, основанные на наблюдении, могут стимулироваться только статистикой?...Если бы медицина не пренебрегла этим инструментом, это означало бы прогресс, она обладала бы большим числом реальных истин, стала бы менее подверженной обвинению в том, что является наукой нетвердых принципов, неуловимых и предположительных»

ПРИЧИНЫ СТАТИСТИЧЕСКИХ ОШИБОК В ПУБЛИКАЦИЯХ ПО БИОМЕДИЦИНЕ
В.П.Леонов (www.biometrica.tomsk.ru)
Низкая статистическая культура исследователей + бедность литературы
Отсутствие в структурах биомедицинских НИИ и вузов специализированных лабораторий и кафедр биостатистики.
Отсутствие отраслевой нормативной базы (отраслевые стандарты, стандарты предприятий и т.п.), регламентирующей этап статистического анализа.
Отсутствие квалифицированной статистической экспертизы в редакциях журналов, в диссертационных и экспертных советах ВАК

XVII век, Голландия
Или: «Врач изучает колбу с мочой, принесённою молодой женщиной»
Или:
Статистические методы
MS Excel + StatPlus
SPSS
STATISTICA 6.0

Предложения (Или что делать?)
1. Курс лекций «Основы медицинской статистики» для слушателей послевузовского образования.
2. Планирование биомедицинского исследования.
3. «Статистическое» рецензирование результатов биомедицинских исследований для отчетов, публикаций и диссертаций

Рекомендации исследователям
Шире применяйте различные статистические методы
Правильно применяйте статистические методы и критерии.
Знайте, любые биомедицинские данные можно статистически обработать.
Проводите многократные расчеты статистических показателей и цензурирование биомедицинских данных.
Корректно выполняйте интерпретацию результатов статистических расчетов
Помните, кроме вас самих никто не сможет сделать правильные научные выводы по результатам исследования!