我正在使用 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)
有什么想法/建议吗?
谢谢