summary（）函数的方法

Question

我在summary()函数的输出上使用mle(stats4)，其输出属于类mle。我想了解summary()如何估算mle(stats4)返回的系数的标准偏差，但我在summary.mle打印的列表中看不到methods(summary)，为什么可以＆＃39;}我发现summary.mle()功能？

（我认为正确的功能是summary.mlm()，但我不确定并且不知道原因mlm，而不是mle）

Answer 1

如

所示

>library(stats4)
>showMethods("summary")
Function: summary (package base)
object="ANY"
object="mle"

以summary方式解释S4。我不知道如何直接检查R中的代码，因此我直接为您搜索stats4的来源。

在stats4/R/mle.R中，有：

setMethod("summary", "mle", function(object, ...){
    cmat <- cbind(Estimate = object@coef,
                  `Std. Error` = sqrt(diag(object@vcov)))
    m2logL <- 2*object@min
    new("summary.mle", call = object@call, coef = cmat, m2logL = m2logL)
})

因此它会创建一个S4对象summary.mle。我猜你现在可以自己追踪代码了。