问题标签 [marginal-effects]

我正在使用 R 来复制一项研究，并获得与作者报告的几乎相同的结果。然而，在某一时刻，我计算的边际效应似乎小得不切实际。如果您能看看我的推理和下面的代码，看看我是否在某个地方或另一个地方弄错了，我将不胜感激。

我的样本包含 24535 个观测值，因变量“x028bin”是取值 0 和 1 的二元变量，此外还有 10 个解释变量。其中九个自变量具有数字级别，自变量“f025grouped”是由不同宗教教派组成的因素。

我想运行一个概率回归，包括宗教教派的假人，然后计算边际效应。为此，我首先消除缺失值并使用因变量和自变量之间的交叉表来验证没有小单元格或 0 单元格。然后我运行运行良好的概率模型，我也获得了合理的结果：

但是，当从概率系数和比例因子计算所有变量的均值时，我获得的边际效应太小（例如 2.6042e-78）。代码如下所示：

很抱歉，由于我的数据集太大，我无法为您提供工作示例。任何评论将不胜感激。非常感谢。

最好的，

托拜厄斯

我正在使用 logit 模型估计随机效应glmer，并且我想报告自变量的边际效应。对于glm模型，包mfx有助于计算边际效应。对象是否有任何包或功能glmer？

谢谢你的帮助。

下面给出了一个可重现的例子

我在 Cross Validated 上问过这个问题，但我认为我可能得不到帮助，因为这更像是一个编程问题，而不是对统计数据的理论/解释。

我正在尝试mlogit在 R 中使用该包，并且一直在关注小插图，试图弄清楚如何获得我的数据的边际效应。提供的示例使用连续变量，但我想知道如何使用分类解释变量来做到这一点。

我的风险值作为协变量是连续的，但我也有年龄、阶级和性别作为协变量。我想看看“女性”或“年轻女性”在风险方面的边际效应。我该怎么做？

帮助文档说：

我不确定如何操纵z数据框来获得女性或年轻女性的平均风险，然后才能计算边际效应。我会分开做吗？我是否以某种方式将数据框按年龄类别划分（比如我只有 2 个年龄类别：年轻人和老年人），这样我就有 1 个年轻人数据框和一个单独的新数据框用于老年人，然后计算平均风险？

我希望从我自己的数据中得到的是能够解释产生我的后代类别的可能性的大小。举个例子，我想说的是，如果风险增加 1 个单位，那么年长的雌性产生 2 个后代的可能性会增加 10%。由于风险增加 1 个单位，年轻女性生育 2 个后代的可能性增加 15%。

我不确定如何手动计算边际效应，因此对如何获得一个包来为我做这件事感到困惑。我也一直在nnet图书馆或图书馆尝试VGAM，但这些似乎都没有提供很大的帮助。

R中的效果包可以用于绘制不包括截距的线性模型的（边际）效果吗？

这是我为此类模型尝试的一些 R 代码：

它产生的错误如下：

如果有人对如何克服此错误有任何想法，请告诉我。

谢谢，

伊莎贝拉

我是 R 新手，我还不完全理解它的计算逻辑......

在以前的帖子的帮助下，我也无法克服我的问题。

我有一个包含 11 个变量的大约 600 个观察值的数据集。我已经成功地运行了多项式模型，但是我无法实现边际效应，因为我的mean()命令正在获取 NA：

数据集：

用于转换数据的命令：

均值矩阵的命令：

我得到的输出：

您能否就此命令的 NA 原因提出建议？

如果它不能为主要是 NAs 的 DEBT 提供正确的输出，我会理解，但为什么它不计算 RSVS、CHINN 和其他人的平均值？

为什么要计算 RGDP 的平均值？两个变量都有：

如何克服这个问题？

更新：

感谢 David Arenburg 的评论，我设法计算出平均值：

但是，在计算边际效应时会出现另一个错误：

在复制平均值 5 次以获得 6 行的矩阵 zz 后，有 6 个备选方案：

effects() 命令给出一个输出：

但是我无法对这个输出得出任何结论——我不知道我可以将这些边际效应分配给哪个变量。

当我尝试一个一个运行时，我得到另一个错误：

模型命令：

我正在使用的软件包：

提前谢谢你， Zyta

我通常使用 mfx 包和 logitmfx 函数生成 logit 模型边际效应。但是，我正在使用的当前调查具有权重（由于在某些人群中过采样，这对样本中 DV 的比例有很大影响），并且 logitmfx 似乎没有任何方法包含权重。

我已经为模型安装了 svyglm，如下所示：

如何从这些结果中产生边际效应？

数据：数据

代码：

问题：

当我运行 mfxboot 函数时，我收到以下错误消息：

bootvals[i, ] <- pdf1 * coef(x1) 中的错误：要替换的项目数不是替换长度的倍数

知道为什么会这样吗？以及如何解决这个问题的任何建议？

谢谢。

在 R 中成功运行 mlogit 模型后，我在尝试获得边际效应时收到错误消息：

我什至尝试按照另一篇文章中的说明更改源代码中的第 16 行，但仍然出现相同的错误。任何帮助，将不胜感激。我所有的变量都不是替代特定的。我有 4 种选择。

以下是网络上发布的类似问题的链接，但我很难使用那里发布的解决方案。 mlogit 在 R 中的边际效应 在此 先感谢

我使用 mlgit() 包在 R 中运行了嵌套 Logit 模型。我现在正在尝试测量边际效应/弹性并继续遇到错误。在这里，我通过修改包作者的小插图重新创建了错误：

我收到以下错误：

当我没有嵌套模型（如常规多项式 Logit）时，这可以正常工作，并且在之前的一些 stackoverflow 问题中已经涵盖了这一点，但是在重新预测更改的数据帧的步骤中特别发生了一些奇怪的事情（在这种情况下意味着框架z）。

我会注意到这里的解决方案：mlogit 在 R 中的边际效应对我没有帮助。

有没有人对如何使用面板校正标准误差在 R 中制作边际效应图有任何建议？

为了估计 R 中的面板校正标准误差，我使用 plm 和 lmtest 包。

首先我用 plm 估计回归模型：

然后我使用 lmtest 包通过 coeftest 传递这个 plm 对象：

plm2 然后包含系数列表和面板校正标准误差。我可以使用什么包和代码然后使用 plm2 创建边际效应图？