用有限元估计法计算id的个数

Question

我正在估算一系列国家和年份的固定效应面板模型。我想计算一下每项估计所包括的国家数目。如果我使用plm函数运行回归，我可以很容易地做到这一点，如下所示。但是，出于其他原因，我不想使用plm运行我的回归。是否有一种方法来计算在一个犯罪估计中使用的id的数量？

下面给出了一个示例代码：

library(plm); library(lfe)
data("Grunfeld", package = "plm")
reg_plm <- plm(inv ~ value + capital,   data = Grunfeld, model = "within", effect = "twoways")
reg_felm <- felm(inv ~ value + capital  | 
                  firm + year, data=Grunfeld)

pdim(reg_plm)[["panel.names"]]$id.names 

请注意，我的数据是一个不平衡的面板，pdim不适用于felm对象。

Answer 1

使用str(reg_felm) --就像在这种情况下，搜索对象时通常使用的那样--显示出固定的效果存储在fe中。

我们可以显示levels

levels(reg_felm$fe$firm)
# [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10"

他们的号码在哪里

length(levels(reg_felm$fe$firm))
# [1] 10

或者使用table来查看它是否平衡。

reg_felm$fe$firm |> table()
#  1  2  3  4  5  6  7  8  9 10 
# 20 20 20 20 20 20 20 20 20 20 

Answer 2

使用nlevels可能更容易

nlevels(reg_felm$fe$firm)