погрешность на ROC-кривой методом бутстрепа |
Здравствуйте, гость ( Вход | Регистрация )
погрешность на ROC-кривой методом бутстрепа |
17.03.2015 - 20:54
Сообщение
#1
|
|
Группа: Пользователи Сообщений: 3 Регистрация: 17.03.2015 Пользователь №: 27109 |
Добрый день.
Недавно столкнулся с статьей, где авторы сделали РОК кривую и отметили на ней погрешности каждой точки методом бутстреп, может кто то сталкивался с подобной задачей и может помочь реализовать этот метод в какой-либо программе (предпочтительно SPSS). Спасибо за ваше время и силы! |
|
17.03.2015 - 20:56
Сообщение
#2
|
|
Группа: Пользователи Сообщений: 1091 Регистрация: 26.08.2010 Пользователь №: 22699 |
Добрый день. Недавно столкнулся с статьей, где авторы сделали РОК кривую и отметили на ней погрешности каждой точки методом бутстреп, может кто то сталкивался с подобной задачей и может помочь реализовать этот метод в какой-либо программе (предпочтительно SPSS). Спасибо за ваше время и силы! Такие графики можно построить в pROC. На форуме я по моему постил такие графики. |
|
17.03.2015 - 21:02
Сообщение
#3
|
|
Группа: Пользователи Сообщений: 3 Регистрация: 17.03.2015 Пользователь №: 27109 |
|
|
17.03.2015 - 21:16
Сообщение
#4
|
|
Группа: Пользователи Сообщений: 1091 Регистрация: 26.08.2010 Пользователь №: 22699 |
Спасибо! Но я пытался загуглить что это такое - безрезультатно... можно чуточку поподробней что такое pROC (как я понимаю это отельная программа)? Это пакет R http://cran.r-project.org/web/packages/pROC/pROC.pdf R есть как расширение для текущих версий практически всех коммерческих статпакетов. Но можно и просто скачать R и использовать отдельно. |
|
17.03.2015 - 21:55
Сообщение
#5
|
|
Группа: Пользователи Сообщений: 3 Регистрация: 17.03.2015 Пользователь №: 27109 |
Это пакет R http://cran.r-project.org/web/packages/pROC/pROC.pdf R есть как расширение для текущих версий практически всех коммерческих статпакетов. Но можно и просто скачать R и использовать отдельно. Спасибо, буду пытаться разобраться как этим воспользоваться! |
|
17.03.2015 - 22:06
Сообщение
#6
|
|
Группа: Пользователи Сообщений: 219 Регистрация: 4.06.2013 Из: Тверь Пользователь №: 24927 |
В SPSS это сделать нельзя ни в самом пакете ни с помощью синтаксиса или скрипта.
Готового ничего нет, нужно писать специально. Сообщение отредактировал anserovtv - 17.03.2015 - 22:15 |
|
17.03.2015 - 22:10
Сообщение
#7
|
|
Группа: Пользователи Сообщений: 1091 Регистрация: 26.08.2010 Пользователь №: 22699 |
Спасибо, буду пытаться разобраться как этим воспользоваться! Используйте функции Код ?ci.se? Compute the CI of sensitivities at given specificities ?ci.sp? Compute the CI of specificities at given sensitivities ?ci.thresholds? Compute the CI of specificity and sensitivity of thresholds ?ci.coords? Compute the CI of arbitrary coordinates или можно прямо из расчета roc Код # Alternatively, you can get the CI directly from roc():
rocobj <- roc(aSAH$outcome, aSAH$s100b, ci=TRUE, of="se", boot.n=100) rocobj$ci # Plotting the CI plot(rocobj) plot(rocobj$ci) |
|
17.03.2015 - 22:12
Сообщение
#8
|
|
Группа: Пользователи Сообщений: 1091 Регистрация: 26.08.2010 Пользователь №: 22699 |
В SPSS это сделать нельзя ни в самом пакете ни с помощью скрипта или макроса. http://www.ibm.com/developerworks/library/ba-call-r-spss/ |
|
17.03.2015 - 22:56
Сообщение
#9
|
|
Группа: Пользователи Сообщений: 219 Регистрация: 4.06.2013 Из: Тверь Пользователь №: 24927 |
Согласен, но я этим никогда не занимался, пока не владею.
Думаю, что это не очень просто. Асимптотический 95% ДИ для площади под ROC-кривой вычисляется в обычной процедуре SPSS. Бутсрепа в этой процедуре нет. Сообщение отредактировал anserovtv - 17.03.2015 - 23:12 |
|
18.03.2015 - 12:43
Сообщение
#10
|
|
Группа: Пользователи Сообщений: 1091 Регистрация: 26.08.2010 Пользователь №: 22699 |
нет там ничего сложного (лишь бы library() была реализована
|
|