ROC-кривая по результатам логистической регрессии Опции

[nokh]

По результатам множественной логистической регрессии программа MedCalc наряду её параметрами и отношениями шансов выдаёт также площадь под характеристической кривой. Когда признак один - всё ясно. Но
(1) Не вполне понятен механизм расчёта чувствительности и специфичности когда в уравнении, например, 5-10 параметров. Неужто алгоритм перебирает все возможные варианты, а кривую строит по наиболее отдалённым от диагонали точкам?
(2) программа MedCalc саму кривую не строит, может кто знает в каком пакете по результатам множественной логистической регрессии строится график ROC-кривой или хотя бы рассчитываются точки для её построения.

[nokh]

В.П. иллюстрирует лекцию своими реальными данными, почему, собственно кривая не может уходить под диагональ? Теоретически, и практически может.

Хм. Пока не посчитал Ваш пример думал не может. Выходит ошибался, беру свои слова про Акелу обратно. Но пример неудачный, я за последние года полтора строил и помогал строить ROC-кривые больше десятка раз и ни разу они не уходили под диагональ.

Теперь относительно того, где считать этот пример. В Attestat логистическая регрессия по 8 предикторам выдается без свободного члена.
В MedCalc полное совпадение с результатами на сайте.

Для того чтобы выдавался свободный член, нужно набить ещё один столбец, заполненный одинаковым числом, и включить его в область анализа. Интуитивно это совсем непонятно, но описано в помощи к AtteStat. Когда тестировал этот модуль таким способом было полное совпадение с MedCalc, а десятичных знаков даже больше.

[Pyrosmani]

Для того чтобы выдавался свободный член, нужно набить ещё один столбец, заполненный одинаковым числом, и включить его в область анализа. Интуитивно это совсем непонятно, но описано в помощи к AtteStat. Когда тестировал этот модуль таким способом было полное совпадение с MedCalc, а десятичных знаков даже больше.

Еще раз убедился в правоте слов Игоря, о том что у нас не принято утруждать себя прочтением manual`ов(я и сам это редко делаю). Действительно при введении дополнительного регрессора в область анализа получаются совсем другие результаты. Кроме указания на необходимость введения дополнит. регрессора в AtteStat есть и наглядный пример как это делать.

Сообщение отредактировал Pyrosmani - 29.05.2010 - 10:12