在逻辑回归的glm函数中，似然函数存储在哪里？ 有家人吗

Question

我目前正在尝试使用optim ， nlm等函数自行运行逻辑回归。但是，与运行glm相比，我在某种程度上为自定义模型获取了不同的系数值。

我怀疑我的似然函数可能有所不同，但是找不到logit似然函数在glm存储的位置。 有人知道吗？

Answer 1

似然函数未存储在glm对象中。 在R 3.1中，仅存在以下字段：

> m <- glm(y ~ x, data=df, family=binomial(link=logit))
> names(m)
 [1] "coefficients"      "residuals"         "fitted.values"    
 [4] "effects"           "R"                 "rank"             
 [7] "qr"                "family"            "linear.predictors"
[10] "deviance"          "aic"               "null.deviance"    
[13] "iter"              "weights"           "prior.weights"    
[16] "df.residual"       "df.null"           "y"                
[19] "converged"         "boundary"          "model"            
[22] "call"              "formula"           "terms"            
[25] "data"              "offset"            "control"          
[28] "method"            "contrasts"         "xlevels"

如果在R提示符下键入glm.fit ，您将看到优化的过程。 特别是，没有单独的似然函数可以用参数空间中的不同点进行评估。 请参阅史蒂夫评论中的链接以了解需要执行的操作。