Хрупкие оценки в объекте coxph
Если кто-то использует obj=coxph(... + frailty(id)), то объект также возвращает (log) оценки слабости для каждого человека, которые можно извлечь с помощью obj$frail.
Кто-нибудь знает, как эти оценки получаются? Это эмпирические байесовские оценки?
Спасибо!
Theodor
1 ответ
Распределение по умолчанию для frailty можно увидеть в ?frailty страница должна быть "гамма". Если вы посмотрите на функцию уязвимости (которая не скрыта), то увидите, что она просто вставляет имя распределения в "уязвимость". и использует get() чтобы получить правильную функцию. Так посмотри на frailty.gamma (также не скрытый), чтобы найти ответы на свой вопрос. Снова оглядываясь на страницу справки, вы видите, что я должен был разобраться со всем этим, не глядя на код, так как он находится вверху страницы. Но есть много путей к знаниям с R. (Они являются ML, а не "эмпирическими байесовскими" оценками.)
Страница справки подсказывает мне, что автор (Терно) ожидает, что вы проконсультируетесь с Терно и Грамбшем для получения дополнительной информации, которая не очевидна при чтении кода. Если вы серьезно работаете с моделями выживания в R, это очень полезная книга. Это очень ясно и полезно для понимания основ пакета "выживания".