问题描述
有人知道我们是否(或何时?)可以在保证金软件包中使用“ at”规范来表示变量?
根据边距vignette(日期为2018-05-22),这并不总是可能的。我引用下面的插图:
“ Stata margins命令的一个非常不错的功能是它可以优雅地处理因子变量。很难在R中模拟此功能,但是margins()函数会尽力而为。 …。 margins()识别因子,并分别显示因子每个级别的边际效应。 (注意:这可能不适用于某些“ at”规范。)”
确实,当我尝试完全依赖于作为因素的自变量(例如下面的代码中的“ LOC”)进行逻辑回归时,会出现错误。
边距(me.model.glm,at = list(LOC = 0:1))
attributes(.Data)1的长度必须与矢量[0]相同
解决方法
试试
margins(me.model.glm,at = list(LOC = c("0","1"))
在类似的情况下,这对我有用。请注意,如果我只为因子变量输入一个值,例如
margins(me.model.glm,at = list(LOC = "0"))
我收到一个错误。我必须为因子变量的每个级别提供值