Как рассчитать с-индекс - PullRequest
0 голосов
/ 29 мая 2018

Существует ли какой-либо код Stata для вычисления c-index (также известного как c-статистика и статистика соответствия) для чистого серого?

Я нашел пакет R, который делает это, но я не знаю, как это сделать в Stata.

1 Ответ

0 голосов
/ 29 мая 2018

Я думаю, вам нужно следующее:

webuse lbw, clear

logit low age lwt i.race smoke ptl ht ui
lroc

Logistic model for low

number of observations =      189
area under ROC curve   =   0.7462

В последней строке указывается c-statistic.

Введите help lroc в командной строке Stata для получения дополнительной информации.


РЕДАКТИРОВАТЬ:

Я только что заметил, что вы упомянули "Fine-Грей ", поэтому я предполагаю, что вы имеете в виду регрессию анализа выживания. команда , предоставленная сообществом somersd делает то, что вы хотите в Stata.

Последнюю версию можно получить здесь:

net describe somersd, from(http://www.rogernewsonresources.org.uk/stata12)

После установки пакетавведите help somersd для получения дополнительной информации и синтаксиса.

Пример использования файла справки автора иллюстрирует его использование:

use http://www.stata-press.com/data/r9/drugtr, clear

generate youth = 100 - age
generate byte censind = 1 - died

somersd studytime drug youth, tr(c) cenind(censind)

Somers' D with variable: studytime
Transformation: Harrell's c
Valid observations: 48

Symmetric 95% CI for Harrell's c
------------------------------------------------------------------------------
             |              Jackknife
   studytime |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        drug |   .7275986   .0367931    19.78   0.000     .6554855    .7997117
       youth |   .6415771   .0528314    12.14   0.000     .5380295    .7451246
------------------------------------------------------------------------------
...