不平衡面板数据集上的指数分数回归模型：处理缺失值

Question

我正在使用不平衡的面板数据集。我的因变量是自然界的。因此，我应用了由（Ramallho等人，2015）开发的面板指数分数回归模型。他们给出了以下命令：

frm(y, X, linkfrac = "logit", var.type = "cluster", var.cluster = id, intercept = F)

我的查询与如何解决面板数据集中的缺失值有关。

首先，我将零值替换为缺失值，但是我发现这不是解决缺失值的正确方法。因此，我将其替换为na并运行以下命令。

frm(y, X, linkfrac = "logit", var.type = "cluster", var.cluster = id, intercept = F)(na = TRUE)

我得到了结果，但是有这个错误：

"Error: attempt to apply non-function"

此外，当我置零时，我的结果与理论非常一致。 “ na”使系数的幅值非常小。请指导我解决这个问题。