Здравствуйте, гость ( Вход | Регистрация )
20.11.2009 - 17:31
Сообщение
#1
|
|
|
Группа: Пользователи Сообщений: 120 Регистрация: 27.08.2009 Пользователь №: 6284 |
День добрый!
Есть таблица Число случаев Год показательБ1 показательБ2 показательБ3 Всего(сумма Б1+Б2+Б3) 2000 2001 2002 .... на пересечении стоят целые числа - количество случаев, выявленных за год. Задача, показать, что есть изменения во временном тренде, т.е. с течением времени. Как вариант - показать, что( к примеру) Б1 растет в удельном количестве по сравнению с Б2, который падает, а Б3 не изменился за эти годы. Есть мысль насчет регресии, но....по ней вопросы. Но если у кого есть опыт, как анализировать эти данные - поделитесь, пожалуйста. Сообщение отредактировал Green - 20.11.2009 - 17:33 ![]() Это не кованализ :)
|
|
|
![]() |
![]() |
![]() |
27.11.2009 - 21:46
Сообщение
#2
|
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
Регрессионные модели для А я сделала и линейную и Пуассона, могу выложить, сравнить оценки, остатки т.т.д. Но не совсем понятно, как с трендом по Neumann.
Для А=средний квадрат (n-1) последовательных разностей =((66-60)^2+(60-64)^2+(64-33)^2+(33-53)^2+(53-38)^2)/5=327,6 Дисперсия=192,3 Отношение =1,7 в таблице критич. знач для n=6 и 5% уровня 0,89 и нулевая гипотеза не отклоняется. Для болезни Б отношение =192/241=0,82 что ниже критического и принимается гипотеза о наличие тренда. Green, что у меня не так? Не понятна для меня идея DoctorStat относительно использования LR для этих данных, поскольку оценить вероятность получить заболевание можно только если есть альтернатива ( не заболел), а в этих данных все больные. |
|
|
![]() |
![]() |
27.11.2009 - 22:25
Сообщение
#3
|
|
|
Группа: Пользователи Сообщений: 381 Регистрация: 18.08.2008 Из: Москва Златоглавая Пользователь №: 5224 |
Не понятна для меня идея DoctorStat относительно использования LR для этих данных, поскольку оценить вероятность получить заболевание можно только если есть альтернатива ( не заболел), а в этих данных все больные. Применение метода отношения правдоподобия для проверки наличия линейного временнОго тренда в вероятности заболевания. Предположим, что распределение наступления болезни подчиняется биномиальному закону. У нас есть результаты испытаний (известные величины) в виде объема выборки n и числа случаев заболеваний m по нескольким последовательным годам. Есть две неизвестные величины, которые нужно рассчитать, вероятность заболевания для одного человека в год p и ее годовой тренд deltap. Основной вопрос исследования: тренд deltap равен нулю или нет (заболеваемость меняется или нет)? Для ответа на этот вопрос нужно проверить статистическую гипотезу методом отношения правдоподобия. Сначала найдем максимальную вероятность L1 (правдоподобие без ограничений) получения наших данных, не налагая ограничений на p и deltap с помощью перемножения отдельных биномиальных вероятностей и дифференцирования произведения по p и deltap (это сделать непросто!). Аналогично найдем максимальную вероятность L0 (правдоподобие с ограничениями) в предположении, что тренда нет: deltap=0. Далее, вычисляем величину Lambda=ln[L1/L0]. Величина Lambda асимптотически (когда наблюдений много) распределена по закону хи-квадрат с (количество неизвестных-1)=1 степенью свободы. Следовательно, выбирая уровень значимости alfa (обычно alfa =0.05), если Lambda больше критического значения хи-квадрат, то тренд присутствует и наоборот.
![]() Просто включи мозги => http://doctorstat.narod.ru
|
|
|
![]() |
![]() |
Green Годовые тренды 20.11.2009 - 17:31
nokh DoctorStat не так давно выкладывал очень хорошую к... 20.11.2009 - 19:22
Green nokh,
Во-первых, спасибо за оперативность.
Во-вто... 20.11.2009 - 23:34
DoctorStat Цитата(Green @ 20.11.2009 - 23:34) К... 21.11.2009 - 11:02
nokh Цитата(Green @ 21.11.2009 - 02:34) .... 23.11.2009 - 07:18
Green DoctorStat, спасибо.
Именно так и думала сделать... 21.11.2009 - 18:21
DrgLena А можно ли подробнее, в каких статистических пакет... 21.11.2009 - 22:26
Green DrgLena,
Речь шла не о Нойманне ("Доверитель... 22.11.2009 - 10:33
DrgLena Критерий Нойманна для тренда из Закса. Я про него,... 22.11.2009 - 12:59
Green nokh, спасибо еще раз!
1. Еще раз напишу пост... 23.11.2009 - 19:46
DoctorStat Цитата(Green @ 23.11.2009 - 19:46) З... 23.11.2009 - 22:25
nokh Цитата(Green @ 23.11.2009 - 22:46) Е... 25.11.2009 - 07:08
DrgLena Green, желательно привести реальные данные, тогда ... 24.11.2009 - 01:35
Green Спасибо DrgLena за все комментарии
1. Фон Нейман ... 24.11.2009 - 14:05
DrgLena Вы хотите грамотно посчитать и оформить результат.... 24.11.2009 - 21:39
Green > Вы многократно повторяете, что у вас есть на ... 25.11.2009 - 00:00
DrgLena Новый поворот мысли, относительные риск посчитать.... 25.11.2009 - 00:50
Green Что значит писать поздно вечером:)
1.В предыдущем... 25.11.2009 - 02:45
DoctorStat Цитата(Green @ 25.11.2009 - 02:45) Я... 25.11.2009 - 11:41
nokh Цитата(DoctorStat @ 25.11.2009 - 14... 25.11.2009 - 11:53
DoctorStat Цитата(nokh @ 25.11.2009 - 11:53) Ши... 25.11.2009 - 14:51
nokh Спасибо, Кендалл и Стюарт есть, читал/листал когда... 25.11.2009 - 15:50
DrgLena Уважаемые коллеги-волонтеры (злопыхателей нет), по... 25.11.2009 - 12:23
DrgLena Цитата(Green @ 25.11.2009 - 03:45) 3... 25.11.2009 - 12:55
Green nokh, спасибо!
Ткнули в нос. Эта работа свалил... 25.11.2009 - 16:07
DrgLena Конечно, ехать, снимается не корректно употребленн... 25.11.2009 - 18:39
nokh ОК, но быстро не сделаю: хочу попробовать ещё пуас... 25.11.2009 - 21:48
nokh Цитата(DrgLena @ 28.11.2009 - 00:46)... 27.11.2009 - 23:03
DrgLena Меня, все же, интересует прогноз, поэтому регресси... 27.11.2009 - 22:17
nokh Цитата(DrgLena @ 28.11.2009 - 01:17)... 27.11.2009 - 23:10
DrgLena Пуассонову регрессию я провела в двух программах и... 28.11.2009 - 00:03
nokh Коэффициенты потому и отличаются, что я анализиров... 28.11.2009 - 00:27
DrgLena Nokh, спасибо, стало понятней, но все же, на бытов... 28.11.2009 - 00:44
nokh Так Вы - врач, я - эколог, откуда браться арифмети... 28.11.2009 - 01:06
DrgLena Это среди математиков, я врач, а среди врачей - ма... 28.11.2009 - 01:24
Green nokh, LrgLena, DoctorStat - всем спасибо!
nok... 30.11.2009 - 12:10
Green DrgLena,
число событий в единицу времени (Пуассон... 30.11.2009 - 13:01
Green nokh,
1. почему вы считали доли для JP?
Что страшн... 30.11.2009 - 14:13
nokh Цитата(Green @ 30.11.2009 - 14:10) В... 30.11.2009 - 17:28
Green nokh,
в пакете JP я выставила модель "using ... 30.11.2009 - 17:57
Green nokh, я еще поспрашиваю.
В дополнение к предыдуще... 1.12.2009 - 21:20
nokh Цитата(Green @ 1.12.2009 - 23:20) ..... 2.12.2009 - 18:43
Green понятно, спасибо!
поскольку занималась анализ... 2.12.2009 - 21:47
Игорь Почему для решения задачи выделения годового тренд... 3.12.2009 - 08:31
Green Только мое мнение:
ССА требует значительно больше... 3.12.2009 - 17:28
Green Игорь, я еще добавлю рассуждение такого плана:
При... 4.12.2009 - 10:24![]() ![]() |