我正在使用mhurdle包来估计截断的正常障碍模型。函数mhurdle包含一个“权重”参数,该参数应该以与lm相同的方式工作(根据mhurdle的帮助)。但是,当我使用这个参数时,我得到的结果与没有权重的情况完全相同。
我想知道为什么会这样,而且,作为替代方案,如果可以在估计模型之前修改变量以包括权重(包使用最大似然估计)。
(我在一个论坛上读过另一个导致同样问题的功能因为权重从未在估计中实际使用过。这里可能是相同但我不知道如何检查而我没有能够再次找到这个问题)
这就是我写函数的方式:
depvar <- mhurdle (y ~ indvar1 + indvar2 + indvar3 |
indvar1 + indvar2 + indvar3 |
0,
data = mydata,
na.action = na.omit,
dist = "tn",
weights = mydata$weights)
有任何想法/建议吗?
由于
答案 0 :(得分:0)
这是一个CRAN包。如果包裹没有像宣传的那样工作,请与作者联系。 https://cran.r-project.org/web/packages/mhurdle/index.html