ROC-анализ, помогите разобраться Опции

[docent]

Доброго времени суток всем профи. Помогите разобраться "чайнику" в ROC-анализе. ROC-кривую строю с помощью пакета SPSS.
У меня много вопросов. Но самый первый - может ли в качестве тестовой переменной выступать переменная, которая имеет интервальную норму?
Т.е. у меня в качестве такой переменной выступает иммуноглобулин А (норма 1,6-5,0). Если можно, то как интерпретировать данные? Как найти пороги отсечения (их же в таком случае должно быть два)?
Привожу исходные данные (0 - здоровый, 1 - больной). Заранее всем спасибо

0 1,4
0 2,6
0 4,1
0 2,6
0 1,4
1 2,8
1 2,8
1 2,8
0 2,8
0 1,0
1 2,1
1 4,4
0 0,9
0 0,8
1 0,9
0 1,2
0 2,2
1 0,9
1 2,8
1 6,7
1 3,0
0 0,6
1 10,1
0 1,8
1 0,8
1 5,0
0 0,8
1 2,8
1 2,8
0 5,0
1 2,8
0 1,8
1 2,8
1 10,7
1 2,3
1 3,4
1 4,8
1 7,9
1 4,6
0 2,3
1 3,1
1 4,3
1 3,8
1 1,4
1 4,1
0 12,4
0 2,2
0 0,9
1 3,1
1 8,2
0 3,1
0 1,4
0 2,7
0 2,5
0 2,2
0 0,9
1 4,9
1 3,2
1 1,8
1 2,2
1 2,3
1 2,0
1 4,3
1 4,8
1 3,5
1 12,1
1 3,5
0 1,8
1 2,1
1 1,4
0 1,7
1 8,2
1 4,5
0 3,1
1 2,1
1 1,9
1 1,7
0 1,2
1 2,3
1 3,5
1 1,2
1 1,1
0 0,8
0 0,9
0 1,4
1 0,8
1 1,2
0 1,5
1 1,4
0 1,1
1 1,3
1 0,8
1 0,6
0 0,9
0 1,1
0 0,5

[p2004r]

Доверительный интервал для AUC

Код

> library(pROC)
> ci(roc(data$V1, data$V2),  method="bootstrap")
Loading required package: tcltk
95% CI: 0.6202-0.8203 (2000 stratified bootstrap replicates)

Принятие оптимального решения в условии равенства частот и "стоимости" ошибок первого и второго рода.

Код

> coords(roc(data$V1, data$V2), "best", ret=c("threshold", "specificity", "sensitivity", "1-npv"))
  threshold specificity sensitivity       1-npv
  2.7500000   0.8421053   0.5862069   0.4285714

Ну и результат применения

Код

> table(data$V1,data$V2>2.75)
  
    FALSE TRUE
  0    32    6
  1    24   34

Сообщение отредактировал p2004r - 16.07.2014 - 01:14

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