Форум врачей-аспирантов

Здравствуйте, гость ( Вход | Регистрация )

> А работает ли аппарат?, ... или изгнание бесов при помощи статистики
Insaniens
сообщение 12.12.2014 - 12:25
Сообщение #1





Группа: Пользователи
Сообщений: 10
Регистрация: 12.12.2014
Пользователь №: 26886



Прошу помощи в статистической обработке материала.

Почему в подзаголовок попали "бесы": есть аппарат, который имеет признаки "лженауки". Итак, исходные данные. Имеется аппарат, который неизвестным науке способом (по заявлению производителей на основании "информативности температуры в определенных точках тела") определяет "131 параметр жизнедеятельности организма человека". В том числе, данное чудо научной (лженаучной?) мысли "показывает" параметры, которые определяются рутинными лабораторными методами.

Быть может, мой ироничный настрой ошибочен и мы имеем научный прорыв, либо данный агрегат годится только для наукообразного "шаманства".

Задача: подтвердить или опровергнуть работоспособность аппарата.

Исходные данные: имеются результаты полученные при помощи аппарата и рутинными лабораторными методиками. У каждого пациента выбранный параметр определялся как в лаборатории, так и на аппарате. "На глаз" видно, что показатели полученные в лаборатории и при помощи "чудо-агрегата" заметно отличаются, как в сторону завышенных, так и заниженных данных (если за эталон взять лабораторию). Для примера фрагмент с данными по гемоглобину (ниже). Первая колонка - показатель гемоглобина, определенный в лаборатории, вторая колонка - у тех же пациентов на "чудо-аппарате".

Вопрос: как статистически подтвердить (или опровергнуть) наличие существенных отклонений в показателях?



Гемоглобин_лаб Гемоглобин
116---------------------139,2
146---------------------160
149---------------------136
121---------------------132
115---------------------153,79
121---------------------135,35
131---------------------125,48
142---------------------77,74
131---------------------155,15
127---------------------143,38
155---------------------136,21
163---------------------149,47
130---------------------100,12
134---------------------133,47


Сообщение отредактировал Insaniens - 12.12.2014 - 12:32
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
 
Открыть тему
Ответов
100$
сообщение 4.01.2015 - 21:55
Сообщение #2





Группа: Пользователи
Сообщений: 902
Регистрация: 23.08.2010
Пользователь №: 22694



Я исхожу из той вероятностно-статистической модели данных, что у каждого пациента в момент измерения гемоглобина есть его истинное значение (которое каждый прибор определяет с присущей ему степенью точности). Наличие в мат. модели измерения любым прибором ссылки на это самое истинное значение не позволяет считать эти измерения статистически независимыми. Оно и понятно: пациент-то один и тот же, сколькими бы приборами его не облепили.

Далее. Если взять математическое ожидание от моделей измерения, которые описаны в посте 23 (я по-прежнему утверждаю, что они кривые), то:
1. мат. ожидание измерения равно его истинному значению (иначе это не прибор, а металлолом).
2. мат. ожидание аддитивной погрешности равно нулю.
3. мат. ожидание истинного значения равно истинному значению (оно и понятно: истинное значение- детерминированная величина).
Числитель к-та корреляции представляет собой ковариацию сл. величин, т.е.: Истинное значение - Истинное значение= мат. ожидание аддитивной ошибки =0.
Иначе говоря, числитель к-та корреляции в случае зависимых (связанных выборок) равен произведению мат. ожиданий аддитивных погрешностей из обоих измерений, т. е. =0. Что нам и продемонстрировали.
Поэтому я утверждаю, что тестировать связанные выборки к-том корреляции (применительно к данной ситуации) - дурной тон.
Кроме того, известные к-ты корреляции (за исключением где-то встреченного мной параболического) линейны. А если статистическая связь двух измерений разными приборами нелинейна, то Пирсоном, Фехнером, Спирменом и Кендаллом ничего получить нельзя.

Ну и на сладкое цитатка из рускоязычной Вики:

Если X,Y независимые случайные величины, то R{X,Y} = 0. Обратное в общем случае неверно.

Сообщение отредактировал 100$ - 4.01.2015 - 22:06
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
p2004r
сообщение 4.01.2015 - 22:34
Сообщение #3





Группа: Пользователи
Сообщений: 1091
Регистрация: 26.08.2010
Пользователь №: 22699



Цитата(100$ @ 4.01.2015 - 21:55) *
Кроме того, известные к-ты корреляции (за исключением где-то встреченного мной параболического) линейны. А если статистическая связь двух измерений разными приборами нелинейна, то Пирсоном, Фехнером, Спирменом и Кендаллом ничего получить нельзя.

Ну и на сладкое цитатка из рускоязычной Вики:

Если X,Y независимые случайные величины, то R{X,Y} = 0. Обратное в общем случае неверно.


позволю себе быть кратким

1. "Натягивать сову на глобус" надо уметь. Это два прибора _поверенные_ (например тупым разведением образцов (о! кстати еще один способ поверить прибор работающий с кровью вообще одним единственным образцом smile.gif )) на выдачу результатов одного и того же показателя. Все погрешности "на нелинейность шкалы" _обязаны_ быть учтены конструктором прибора. Единственный ожидаемый результат это линейная зависимость показаний обоих методов.

2. R -- _выборочный_ показатель, он может быть _любым_ в пределах некоего своего _распределения_. Получить оценку этого распределения очень просто для экспериментальных данных и невероятно трудно в случае высасывания из пальцатеоретическом случае. Соответственно делать предположения о том, с какого уровня ошибки можно считать наблюдаемое достоверным событием, мы в своем (обычном повседневном) праве.


Signature
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
100$
сообщение 5.01.2015 - 01:53
Сообщение #4





Группа: Пользователи
Сообщений: 902
Регистрация: 23.08.2010
Пользователь №: 22694



Цитата
R -- _выборочный_ показатель


Золотые слова. И выборочный и теоретический. И даже известны условия, при которых выборочный сходится к теоретическому.

Цитата
он может быть _любым_


Вот, поистине, бесценные сведения. Под словой "любым" , надо думать, понимается диапазон - 1 ... 1

Цитата
в пределах некоего своего _распределения_


Асимптотически нормального ~N(0,1)

Цитата
Получить оценку этого распределения очень просто для экспериментальных данных


Ресэмплинговыми техниками. Кто б спорил.

Цитата
и невероятно трудно в случае высасывания из пальцатеоретическом случае.


Матка боска! (польск.) Ужасть-то какая! Тем не менее все это сделано.

Цитата
Соответственно делать предположения о том, с какого уровня ошибки можно считать наблюдаемое достоверным событием, мы в своем (обычном повседневном) праве.


Красотища-то какая! Случайная величина в результате эксперимента (измерения гемоглобина) приняла некое численное значение. Глазам своим не верю и считаю это невозможным событием.

Неотвержение нулевой гипотезы формально не доказывает ее истинность.

Сообщение отредактировал 100$ - 5.01.2015 - 02:11
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
p2004r
сообщение 5.01.2015 - 12:19
Сообщение #5





Группа: Пользователи
Сообщений: 1091
Регистрация: 26.08.2010
Пользователь №: 22699



Цитата(100$ @ 5.01.2015 - 01:53) *
Золотые слова. И выборочный и теоретический. И даже известны условия, при которых выборочный сходится к теоретическому.



Вот, поистине, бесценные сведения. Под словой "любым" , надо думать, понимается диапазон - 1 ... 1



Асимптотически нормального ~N(0,1)



Ресэмплинговыми техниками. Кто б спорил.



Матка боска! (польск.) Ужасть-то какая! Тем не менее все это сделано.



Красотища-то какая! Случайная величина в результате эксперимента (измерения гемоглобина) приняла некое численное значение. Глазам своим не верю и считаю это невозможным событием.

Неотвержение нулевой гипотезы формально не доказывает ее истинность.


1. Нако выряйте еще отдельных букв из моих слов, может из них сложится еще что то ...

Трудность отвержения отсутствия связи между измерениями "вообще" у топикстартера (как и пути её преодоления) я описал еще в первых постах этого треда.

Тем не менее, мы точно так же считаем массу событий (завтра взойдет солнце) совершенно произвольно "достоверными". В данном случае мы примем возможную ошибку выхода за границу оценки доверительного интервала как несущественную, и будем считать событие "величина находится внутри интервала" достоверным для своих дальнейших рассуждений.

Таким образом величина возможного совпадения результатов не заключена в проверке какой то там гипотезы, а в расчете величины доверительного интервала для некоторой характеристики оцениваемого прибора. И коэффициент корреляции (вернее его доверительный интервал) более чем подходит для сравнения двух приборов выдающих на гора величину в заявленных "грамм гемоглобина на объем". (было бы "много гемоглобина", "мало гемоглобина", "вся сила в гемоглобине" было бы что обсуждать ещё)

2. Асимптотически нормального ~N(0,1)

На практике данное распределение может быть каким угодно у конкретной выборки. Поэтому я пишу у "выборочного". А что там, да кабы поверьте абсолютно не интересно... так же как и рассказываемые вами неумные скетчи (над которыми я не смеялся даже когда слышал их в первый раз).


Signature
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 

Сообщений в этой теме
- Insaniens   А работает ли аппарат?   12.12.2014 - 12:25
- - p2004r   Учебные задачки надо самому решать вообщето. А та...   12.12.2014 - 13:34
|- - Insaniens   Цитата(p2004r @ 12.12.2014 - 13:34) ...   13.12.2014 - 15:55
|- - p2004r   Преамбула: Ну вот редкий из "медистов" с...   13.12.2014 - 17:33
|- - Insaniens   Цитата(p2004r @ 13.12.2014 - 17:33) ...   13.12.2014 - 20:19
|- - p2004r   Цитата(Insaniens @ 13.12.2014 - 20:1...   13.12.2014 - 20:50
|- - Insaniens   Цитата(p2004r @ 13.12.2014 - 20:50) ...   14.12.2014 - 14:30
- - DoctorStat   Цитата(Insaniens @ 12.12.2014 - 12:2...   13.12.2014 - 13:31
|- - Insaniens   Цитата(DoctorStat @ 13.12.2014 - 13...   13.12.2014 - 15:56
|- - DoctorStat   Цитата(DoctorStat @ 13.12.2014 - 13...   4.01.2015 - 21:30
|- - p2004r   без разделения fixed и random эффекта этот расчет ...   4.01.2015 - 22:39
- - Диагностик   Цитата(Insaniens @ 12.12.2014 - 17:2...   27.12.2014 - 07:07
|- - 100$   Цитата(Диагностик @ 27.12.2014 - 08...   27.12.2014 - 13:27
||- - p2004r   Цитата(100$ @ 27.12.2014 - 13:2...   27.12.2014 - 16:19
|- - Диагностик   Цитата(Диагностик @ 27.12.2014 - 12...   1.01.2015 - 06:39
|- - 100$   Цитата(Диагностик @ 1.01.2015 - 07:3...   1.01.2015 - 18:08
|- - nokh   Цитата(Диагностик @ 1.01.2015 - 08:3...   6.01.2015 - 00:39
|- - p2004r   Тут принципиально получается может быть несколько ...   7.01.2015 - 19:30
- - DrgLena   ?Статисты? ?сатанисты продолжают издеваться над м...   27.12.2014 - 18:08
- - 100$   А вообще очень забавно: 1. Критерий типа омега-кв...   27.12.2014 - 19:44
|- - Диагностик   Цитата(100$ @ 28.12.2014 - 00:4...   28.12.2014 - 13:55
|- - 100$   Цитата(Диагностик @ 28.12.2014 - 14...   28.12.2014 - 18:59
- - 100$   Ну никак не дает мне модератор культурно пообщатьс...   29.12.2014 - 10:12
|- - p2004r   Цитата(100$ @ 29.12.2014 - 10:1...   29.12.2014 - 12:37
|- - 100$   Цитата(p2004r @ 29.12.2014 - 13:37) ...   29.12.2014 - 13:41
||- - p2004r   Цитата(100$ @ 29.12.2014 - 13:4...   29.12.2014 - 17:27
||- - 100$   Цитата(p2004r @ 29.12.2014 - 18:27) ...   29.12.2014 - 17:40
|- - Диагностик   Цитата(p2004r @ 29.12.2014 - 17:37) ...   29.12.2014 - 15:51
- - Олег Кравец   5 суток. Всем необходимо остыть и расслабиться. И ...   2.01.2015 - 18:35
- - p2004r   Цитата(Insaniens @ 12.12.2014 - 12:2...   3.01.2015 - 13:34
- - 100$   Здесь надо для каждого пациента несколько раз изме...   3.01.2015 - 22:11
|- - Диагностик   Цитата(100$ @ 4.01.2015 - 03:11...   4.01.2015 - 10:59
- - DrgLena   Цитата(Диагностик @ 4.01.2015 - 10:5...   4.01.2015 - 12:20
|- - p2004r   Цитата(DrgLena @ 4.01.2015 - 12:20) ...   5.01.2015 - 19:14
- - Диагностик   Показано конечно и очень коряво, но почему Вы не у...   4.01.2015 - 12:39
- - DrgLena   Цитата(Диагностик @ 4.01.2015 - 12:3...   4.01.2015 - 14:00
- - Диагностик   Из их сообщений этого не следует.   4.01.2015 - 15:44
- - DrgLena   Да уж, действительно забавно, не знаком ( по ссылк...   4.01.2015 - 18:05
- - 100$   Я исхожу из той вероятностно-статистической модели...   4.01.2015 - 21:55
|- - p2004r   Цитата(100$ @ 4.01.2015 - 21:55...   4.01.2015 - 22:34
|- - 100$   ЦитатаR -- _выборочный_ показатель Золотые слова....   5.01.2015 - 01:53
|- - p2004r   Цитата(100$ @ 5.01.2015 - 01:53...   5.01.2015 - 12:19
||- - 100$   Цитата(p2004r @ 5.01.2015 - 13:19) 1...   5.01.2015 - 14:07
|- - Диагностик   Цитата(100$ @ 5.01.2015 - 06:53...   6.01.2015 - 15:25
- - p2004r   Господин Модератор! Мне самовыпилиться с этог...   5.01.2015 - 21:20
- - anserovtv   В некоторых пакетах программ (XLStat) имеется про...   6.01.2015 - 11:26
- - Insaniens   Уважаемые участники форума! Благодарю за выска...   6.01.2015 - 15:09
|- - Диагностик   Цитата(Insaniens @ 6.01.2015 - 20:09...   6.01.2015 - 15:39
- - Insaniens   Ага Именно коэффициент корреляции я и использовал...   6.01.2015 - 16:37
|- - 100$   Insaniens, не в службу, а в дружбу сообщите, пож-с...   6.01.2015 - 20:22
|- - Insaniens   Цитата(100$ @ 6.01.2015 - 20:22...   7.01.2015 - 10:53
|- - Диагностик   Цитата(Insaniens @ 7.01.2015 - 15:53...   7.01.2015 - 15:14
|- - Insaniens   Цитата(Диагностик @ 7.01.2015 - 15:1...   14.01.2015 - 10:28
|- - nokh   Цитата(Insaniens @ 14.01.2015 - 12:2...   19.01.2015 - 15:42
- - Диагностик   Я сейчас решаю подобную задачу, поэтому и зантерес...   8.01.2015 - 12:02
- - Insaniens   Цитата(nokh @ 19.01.2015 - 15:42) Во...   20.01.2015 - 16:59


Добавить ответ в эту темуОткрыть тему