Помощь - Поиск - Пользователи - Календарь
Полная версия этой страницы: погрешность на ROC-кривой методом бутстрепа
Форум врачей-аспирантов > Разделы форума > Медицинская статистика
KFateev
Добрый день.

Недавно столкнулся с статьей, где авторы сделали РОК кривую и отметили на ней погрешности каждой точки методом бутстреп, может кто то сталкивался с подобной задачей и может помочь реализовать этот метод в какой-либо программе (предпочтительно SPSS).




Спасибо за ваше время и силы!
p2004r
Цитата(KFateev @ 17.03.2015 - 20:54) *
Добрый день.

Недавно столкнулся с статьей, где авторы сделали РОК кривую и отметили на ней погрешности каждой точки методом бутстреп, может кто то сталкивался с подобной задачей и может помочь реализовать этот метод в какой-либо программе (предпочтительно SPSS).




Спасибо за ваше время и силы!


Такие графики можно построить в pROC. На форуме я по моему постил такие графики.
KFateev
Цитата(p2004r @ 17.03.2015 - 21:56) *
Такие графики можно построить в pROC. На форуме я по моему постил такие графики.



Спасибо! Но я пытался загуглить что это такое - безрезультатно... можно чуточку поподробней что такое pROC (как я понимаю это отельная программа)?
p2004r
Цитата(KFateev @ 17.03.2015 - 21:02) *
Спасибо! Но я пытался загуглить что это такое - безрезультатно... можно чуточку поподробней что такое pROC (как я понимаю это отельная программа)?


Это пакет R http://cran.r-project.org/web/packages/pROC/pROC.pdf

R есть как расширение для текущих версий практически всех коммерческих статпакетов. Но можно и просто скачать R и использовать отдельно.
KFateev
Цитата(p2004r @ 17.03.2015 - 22:16) *
Это пакет R http://cran.r-project.org/web/packages/pROC/pROC.pdf

R есть как расширение для текущих версий практически всех коммерческих статпакетов. Но можно и просто скачать R и использовать отдельно.



Спасибо, буду пытаться разобраться как этим воспользоваться!
anserovtv
В SPSS это сделать нельзя ни в самом пакете ни с помощью синтаксиса или скрипта.
Готового ничего нет, нужно писать специально.
p2004r
Цитата(KFateev @ 17.03.2015 - 21:55) *
Спасибо, буду пытаться разобраться как этим воспользоваться!


Используйте функции
Код
       ?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)
p2004r
Цитата(anserovtv @ 17.03.2015 - 22:06) *
В SPSS это сделать нельзя ни в самом пакете ни с помощью скрипта или макроса.


http://www.ibm.com/developerworks/library/ba-call-r-spss/
anserovtv
Согласен, но я этим никогда не занимался, пока не владею.
Думаю, что это не очень просто.
Асимптотический 95% ДИ для площади под ROC-кривой вычисляется в обычной процедуре SPSS.
Бутсрепа в этой процедуре нет.
p2004r
нет там ничего сложного (лишь бы library() была реализована smile.gif
Для просмотра полной версии этой страницы, пожалуйста, пройдите по ссылке.
Форум IP.Board © 2001-2025 IPS, Inc.