![]() |
Здравствуйте, гость ( Вход | Регистрация )
![]() |
![]()
Сообщение
#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 |
|
![]() |
![]() |
![]() |
![]()
Сообщение
#2
|
|
![]() Группа: Пользователи Сообщений: 1091 Регистрация: 26.08.2010 Пользователь №: 22699 ![]() |
Учебные задачки надо самому решать вообщето.
А так вот например --- рандомизацией получаем перцентили необходимые для отбрасывания "нулевой гипотезы" о отсутствии связи. Критерий не достигает реальное значение --- значит связь не может быть доказана. Отстутвие связи доказать значительно сложнее, так как она может быть очень незначительная и потребуется значительно большая группа (её размер вы думаю посчитаете сами?) Код > quantile(replicate(10000, cor(gembn$Гемоглобин, sample(gembn$Гемоглобин_лаб))), probs=c(0.025,0.5,0.975)) 2.5% 50% 97.5% -0.55839229 0.01049748 0.50118705 > quantile(replicate(100000, cor(gembn$Гемоглобин, sample(gembn$Гемоглобин_лаб))), probs=c(0.025,0.5,0.975)) 2.5% 50% 97.5% -0.56194921 0.01313884 0.49301626 > quantile(replicate(100000, cor(gembn$Гемоглобин, sample(gembn$Гемоглобин_лаб))), probs=c(0.025,0.5,0.975)) 2.5% 50% 97.5% -0.56014658 0.01265719 0.49350855 > quantile(replicate(1000000, cor(gembn$Гемоглобин, sample(gembn$Гемоглобин_лаб))), probs=c(0.025,0.5,0.975)) 2.5% 50% 97.5% -0.55927770 0.01301754 0.49402244 Ладно, вот так можно определеить объем выборки для "доказательсва значимости" такого небольшого по ожидаемой величине эффекта (естетственно что и сам эффект будет меняться и ещё как ![]() Код > quantile(replicate(10000, cor(c(gembn$Гемоглобин, gembn$Гемоглобин), sample(c(gembn$Гемоглобин_лаб, gembn$Гемоглобин_лаб)))), probs=c(0.025,0.5,0.975))
2.5% 50% 97.5% -0.3871983241 -0.0001393141 0.3637828166 > quantile(replicate(10000, cor(c(gembn$Гемоглобин, gembn$Гемоглобин), sample(c(gembn$Гемоглобин_лаб, gembn$Гемоглобин_лаб)))), probs=c(0.025,0.5,0.975)) 2.5% 50% 97.5% -0.384000417 0.006484305 0.362028738 > quantile(replicate(10000, cor(c(gembn$Гемоглобин, gembn$Гемоглобин, gembn$Гемоглобин), sample(c(gembn$Гемоглобин_лаб, gembn$Гемоглобин_лаб, gembn$Гемоглобин_лаб)))), probs=c(0.025,0.5,0.975)) 2.5% 50% 97.5% -0.308052837 0.004971756 0.300314873 ![]() |
|
![]() |
![]() |
![]() ![]() |