我已经在R中使用glm函数开发了一个二项式逻辑回归。我需要三个输出,它们是对数似然性(无系数(对数似然性(仅限常数(对数似然(最优(
我需要什么功能或程序包才能获得这些输出?
假设我们有一个合适的型号m
。
完整模型(即MLE(的对数似然:logLik(m)
仅拦截模型的对数似然:logLik(update(m, . ~ 1))
尽管如果我们足够仔细地考虑deviance()
和$null.deviance
分量(这些分量是关于饱和模型定义的(