R:GLM中饱和模型的对数似然

时间:2018-03-05 06:45:13

标签: r glm log-likelihood

LL = loglikelihood

剩余偏差 = 2(LL(饱和模型) - LL(建议模型)

但是,当我使用glm函数时,似乎

Residual Deviance = -2LL(Proposed Model)

例如,

mydata <- read.csv("https://stats.idre.ucla.edu/stat/data/binary.csv")
mydata$rank <- factor(mydata$rank)
mylogit <- glm(admit ~ gre + gpa + rank, data = mydata, family = "binomial")
summary(mylogit)
###
Residual deviance: 458.52  on 394  degrees of freedom
AIC: 470.52
#Residual deviance
-2*logLik(mylogit)
##'log Lik.' 458.5175 (df=6)
#AIC
-2*logLik(mylogit)+2*(5+1)
##470.5175

LL(饱和模型)在哪里?如何在R中获得它的值?

谢谢。

1 个答案:

答案 0 :(得分:0)

我得到了答案:它只发生在饱和模型的对数似然为0时,对于离散模型意味着饱和模型下观测数据的概率为1.二进制数据几乎就是这是真的情况(因为个别拟合概率变为零或一)。HHere以获取详细信息。