用于获取二元逻辑回归的对数似然的 R 代码

问题描述

我在 R 中使用 glm 函数开发了一个二项逻辑回归。我需要三个输出 对数似然（无系数） 对数似然（仅常量） 对数似然（最佳）

我需要哪些函数或包才能获得这些输出？

解决方法

假设我们有一个拟合模型 m。

完整模型的对数似然（即在 MLE）：logLik(m)

仅拦截模型的对数似然：logLik(update(m,. ~ 1))

尽管如果我们足够仔细地考虑 deviance() 和 $null.deviance 组件（这些是针对 饱和模型 )