інтерпретація оцінок засмічення логістичної регресії


21

Чи міг би хтось порадити мені, як інтерпретувати оцінки з логістичної регресії, використовуючи посилання забивання?

Я встановив таку модель у lme4:

glm(cbind(dead, live) ~ time + factor(temp) * biomass,
    data=mussel, family=binomial(link=cloglog))

Наприклад, оцінка часу - 0,015. Чи правильно сказати, що шанси смертності за одиницю часу множать на exp (0,015) = 1,015113 (~ 1,5% збільшення за одиницю часу).
Іншими словами, чи є оцінки, отримані в засміченні, виражені в коефіцієнтах журналу, як це стосується логістичної логістичної регресії?


Відредагуйте код, щоб дотримуватися Rправил синтаксису. Ти не можеш (після '
Френк Харрелл

Відредагуйте оригінальну публікацію та видаліть коментар.
Френк Харрелл

Відповіді:


30

Функція зв’язку доповнення лог-журнал не є логістичною регресією - термін "логістичний" означає посилання logit. Це все ще біноміальна регресія.

оцінка часу 0,015. Чи правильно сказати, що шанси смертності на одиницю часу множать на exp (0,015) = 1,015113 (~ 1,5% збільшення за одиницю часу)

Ні, тому що він не моделює з точки зору журналів. Це те, що ви мали б із посиланням logit; якщо ви хочете, щоб модель працювала з точки зору логіки, скористайтеся logit-посиланням.

Про це говорить функція зв’язку допоміжного журналу-журналу

η(х)=журнал(-журнал(1-πх))=хβ

де .πх=П(Y=1|Х=х)

Отже, - не коефіцієнт шансів; дійсно exp ( η ) = - log ( 1 - π x ) .досвід(η)досвід(η)=-журнал(1-πх)

досвід(-досвід(η))=(1-πх)1-досвід(-досвід(η))=πхх

х


Як Бен обережно натякав у своєму питанні в коментарях:

чи правда сказати, що ймовірність смертності за одиницю часу (тобто небезпеку) збільшується на 1,5%?

Параметри в додатковій моделі журналу журналів мають чітку інтерпретацію з точки зору небезпеки. У нас це є:

еη(х)=-журнал(1-πх)=-журнал(Sх)S

(Отже, виживання журналу зменшиться приблизно на 1,5% за одиницю часу у прикладі.)

год(х)=-ггхжурнал(Sх)=ггхеη(х) , тому справді здається, що у прикладі, наведеному у питанні, ймовірність смертності * за одиницю часу збільшується приблизно на 1,5%

П(Y=1)


7
чи правда сказати, що ймовірність смертності за одиницю часу (тобто небезпеку) збільшується на 1,5%?
Бен Болкер
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.