погрешность на ROC-кривой методом бутстрепа Опции

[KFateev]

Добрый день.

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




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

Эскизы прикрепленных изображений

 

[p2004r]

Добрый день.

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




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

Такие графики можно построить в pROC. На форуме я по моему постил такие графики.

[KFateev]

Такие графики можно построить в pROC. На форуме я по моему постил такие графики.

Спасибо! Но я пытался загуглить что это такое - безрезультатно... можно чуточку поподробней что такое pROC (как я понимаю это отельная программа)?

[p2004r]

Спасибо! Но я пытался загуглить что это такое - безрезультатно... можно чуточку поподробней что такое pROC (как я понимаю это отельная программа)?

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

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

[KFateev]

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

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

Спасибо, буду пытаться разобраться как этим воспользоваться!

[p2004r]

Спасибо, буду пытаться разобраться как этим воспользоваться!

Используйте функции

Код

       ?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)