Question

以下列模型为例

library('lme4')
foo <- lmer(cty ~ hwy + (1|model), data=mpg, REML=F)

现在我们可以检索模型的LogLikelihood

sum.foo <- summary(foo)
LL <- sum.foo[["logLik"]]
LL
'log Lik.' -343 (df=4)
typeof(LL)
[1] "double"

问题：如何将所有这些信息存储在双重对象中？似乎至少有两个字符串，即“log Lik”。和“（df = 4）”。有没有办法从这个对象中检索df的值？

Answer 1

LL实际上只是一个数字（带属性）。打印LL时看到的行由print.logLik函数创建，如下所示：

getAnywhere(print.logLik)
A single object matching ‘print.logLik’ was found
It was found in the following places
  registered S3 method for print from namespace stats
  namespace:stats
with value

function (x, digits = getOption("digits"), ...) 
{
    cat("'log Lik.' ", paste(format(c(x), digits = digits), collapse = ", "), 
        " (df=", format(attr(x, "df")), ")\n", sep = "")
    invisible(x)
}

这是运行LL时调用的内容，但LL实际上只有一个数字（长度为1的数字向量）。 cat函数用于打印您在控制台上看到的'log Lik.' -343 (df=4)。

为了获得你可以做的df值（从上面的函数中也可以看出）：

format(attr(LL[['logLik']], "df"))

查看示例（来自lmer的文档）：

fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
a <- summary(fm1)# (with its own print method)

> a[['logLik']]
'log Lik.' -871.8141 (df=6)
> format(attr(a[['logLik']], "df"))
[1] "6"

根据@BenBolker的评论format中的提及只会将其转换为字符。

> attr(a[['logLik']], "df")
[1] 6

可能更好。

lme4 :: lmer摘要对象包含带有字符串的双重对象

1 个答案: