#r #regression #cox-regression #goodness-of-fit
#r #регрессия #cox-регрессия #соответствие
Вопрос:
Я создал следующую модель:
data("nwtco")
subcoh <- nwtco$in.subcohort
selccoh <- with(nwtco, rel==1|subcoh==1)
ccoh.data <- nwtco[selccoh,]
ccoh.data$subcohort <- subcoh[selccoh]
ccoh.data$age <- ccoh.data$age/12 # Age in years
fit.ccSP <- cch(Surv(edrel, rel) ~ stage histol age,
data =ccoh.data,subcoh = ~subcohort, id=~seqno, cohort.size=4028, method="LinYing")
Однако я просто не могу понять, как получить информацию о том, насколько хорошо модель соответствует данным.
Кто-нибудь знает, как это сделать?
Комментарии:
1. Похоже, это не конкретный вопрос программирования, который подходит для переполнения стека. Если у вас есть общие вопросы об оценке соответствия статистической модели, вам следует задать такие вопросы в разделе Перекрестная проверка . У вас больше шансов получить там лучшие ответы.
2. Спасибо за информацию. Проблема, однако, заключается в том, как конкретно получить информацию из R, то есть какой код использовать, поскольку cch-функция усложняет ситуацию.
3. Но какую информацию вы пытаетесь получить? Как вы хотите проверить правильность подгонки? Вы должны указать, какое точное значение вы пытаетесь вычислить в вопросе.
4. Я понимаю. В идеале, я бы, например, хотел иметь возможность выводить остатки Кокса-Снелла или Мартингейла.
5. Возможно
riskRegression
, пакет мог бы вам помочь, но у меня есть вопрос. Вы используете модель пропорционального риска Кокса?