问题标签 [r-zelig]

我想编写一个涉及向 Zelig 提供参数值的函数。但是，在我这样做之前，我需要弄清楚如何将值传递给 Zelig，即使没有函数。这是我想要的结果：

请注意race已成功设置为“其他”并educate已正确设置为 9 的输出。但是，如果我尝试提供带有这些参数的命名向量，它将不起作用：

输出显示race并educate保持默认值（可能是原始数据集中这些变量的模式和均值）。我也尝试过提供命名列表。我该如何克服呢？

此外，我似乎根本无法修改这个data论点。文档说我应该能够通过单独命名参数（使用...参数）或通过提供一个新的 data.frame 来进行平均（使用data参数）来设置参数。这里也没有运气：

我正在使用 R 3.6.3、Zelig 版本 5.1.7、macOS Mohave 10.14.5。谢谢！

我想估计一个 SUR（看似不相关的回归）模型。

我尝试使用systemfit和它的包装器Zelig。但我无法理解如何指定要预测的因素（即添加固定效应）和聚类标准误差，就像我们在felm().

此外，如果我只是将固定效应变量添加到回归方程中，则会出现以下错误：

非常感谢你的帮助！

我正在从我的数据中添加一个数据样本：

（由于这是一个样本数据集，因此，上述两个回归会给出我上面提到的错误，但在我的实际数据上运行良好。）

但是，我感兴趣的是使用 SUR 估计上述公式，其中X_var8固定X_var9效应和标准误差聚集在X_var8水平上。

如果我们使用felm()，规范是

但是，由于我的标准误差与方程式相关，因此我需要使用 SUR。

任何帮助将非常感激。谢谢你！

我正在尝试运行由另一位研究人员编写的 R 代码，但在基于 zelig 函数建模的后验模拟返回风险比时，它总是给我错误。复制的材料可以在这里下载。更具体地说，我无法运行“_replic.ind2006.AJPS”代码，该代码使用“individualdataset2006”数据（它都在 zip 存档中）。

以下代码是从第 74 行到第 92 行编写的。

这部分从使用匹配程序开始，然后是通过 logit 函数估计模型。那么问题从哪里开始呢？它在运行循环时精确开始

其中length(inc)表示分类变量的级别数，这意味着，在实践中，1:4。

代码无法运行

它应该存储来自估计的风险比率。它给出了以下错误：

rr[i] <- summary(sim(fit1, x = bf0, x1 = bf1))$qi.stats$rr[1] 中的错误：替换长度为零

我首先认为这是一个循环问题，如果我“手动”存储风险比，但运行

独自返回NULL。事实上，如果我只返回没有 的摘要$qi.stats$rr[1]，则输出中没有风险比，它只给我期望值、预测值和一阶差异。

运行时也报“replacement has zero length”的问题

但不同的是我可以通过 得到期望值get_qi，所以现在这不是一个大问题，虽然它很麻烦。

有人对此有任何线索吗？

先感谢您。