命令使用有序概率的阈值

时间:2016-07-22 15:30:08

标签: r stata

我估计一个有序的概率(对于那些只知道概率的人我也在下面的背面超链接中添加了一个非常简短的解释)。但是,我的因变量是一个百分比,已分为八个百分比组。意思我知道,例如类别1表示0%,类别2表示0<y<5等。因此,我知道所有阈值alpha并且可以在my likelihood function中使用它们(参见本背景中的等式(2) -链接)。有人在R或Stata中知道这个命令还是存在这样的命令?

1 个答案:

答案 0 :(得分:0)

我认为你可以在Stata中使用oprobitoglm这样做:

webuse nhanes2f
constraint 1 [cut1]_cons=-3
constraint 2 [cut2]_cons=-2 
constraint 3 [cut3]_cons=-1 
constraint 4 [cut4]_cons=0
oprobit health female black age c.age#c.age, constraint(1 2 3 4)
oglm health female black age c.age#c.age, link(probit) constraint(1 2 3 4)

Stata使用Pr(y = j | x)= Pr(cut_ {j-1}&lt; x&#39; b + u&lt; = cut_ {j}),因此这可能不适合你的箱子。您可能需要从切割点添加/减去c(mindouble)以获得您想要的结果。

这在这里是偏离主题的,但我认为顺序方法与观察你的&#34;潜在的&#34;变量的范围有限且规模未知。

我会尝试intreg作为健壮性检查。它不涉及范围问题,但规模问题在该设置中不是问题。