0

我正在使用 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)

有什么想法/建议吗?

谢谢

4

1 回答 1

0

这是一个 CRAN 包。如果包没有像宣传的那样工作,请联系作者。https://cran.r-project.org/web/packages/mhurdle/index.html

于 2015-07-24T19:28:19.627 回答