使用 r 中的预测概率查找对数似然

【问题标题】：Find log-likelihood using the predicted probabilities in r使用 r 中的预测概率查找对数似然
【发布时间】：2021-06-26 17:13:01
【问题描述】：

我建立了一个glm 模型并用它来预测测试数据的概率。

model = glm(y ~ x, data=dt, family=binomial(link='logit'))
pred = predict(model, newdata=test.dt, type='response')

如何找到预测概率的检验对数似然？

【问题讨论】：

【解决方案1】：

类似

dbinom(dt$y, prob=pred, size=1, log=TRUE)

应该给你对数似然（假设这些是二元/伯努利响应）。您也可以直接将对数似然计算为dt$y*log(pred) + (1-dt$y)*log(1-pred)（我认为）

如果您想要组合对数似然，只需 sum() 上述值 ... 或 logLik(model)

【讨论】：