从 plm() 中的随机效应模型访问 theta

Question

我正在运行建立一个随机效应模型，想获取 theta 对象（准去均值方差）并将其存储到一个对象中。环顾四周，不清楚它存储在 plm() 对象中的内容，或者当我运行 summary() 函数时它来自哪里。

这是一个代码示例：

library(plm)
data("Produc", package = "plm")
zz <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp, 
    data = Produc, index = c("state", "year"),
    model = "random")

summary(zz)

...
Effects:
                   var  std.dev share
idiosyncratic 0.001454 0.038137 0.175
individual    0.006838 0.082691 0.825
theta: 0.8888
...

我想我可以通过 theta <- zz$theta 访问它，但是 returns NULL。

任何更熟悉这个包如何工作的人如果能提供建议，我们将不胜感激！

Answer 1

模型对象的ercomp元素包含各种分量误差的估计。您可能是想从 zz$ercomp$theta 获取 theta；但是，您原来的 post?

中没有随机效应模型

从 plm() 中的随机效应模型访问 theta

Access theta from random effect model in plm()

r

plm