Какое значение приведено в сводке модели Кокша в R? Например,
Rsquare= 0.186 (max possible= 0.991 )
Я по глупости включил его в рукопись в качестве значения и рецензент вскочил на него, сказав, что он не знал об аналоге статистики из классической линейной регрессии, разрабатываемой для модели Кокса, и, если она была, пожалуйста, предоставить ссылку. Любая помощь будет отличной!
Ответы:
Используя
getS3method("summary","coxph")
вы можете посмотреть, как он рассчитывается.Соответствующие строки кода следующие:
Здесь
cox$loglik
это «вектор длины 2, содержащий логарифмическое правдоподобие с начальными значениями и с конечными значениями коэффициентов» (см.?coxph.object
) Иcox$n
«число наблюдений, использованных в подгонке».источник
Деление на числа наблюдений в резюме неправильно, это должно быть числом без цензуры событий; см. O'Quigley et al. (2005) Объясненная случайность в моделях пропорциональных опасностей Статистика в медицине с. 479-489.n
coxph
источник