順序付きプロビットのしきい値を使用するコマンド

Question

私は順序付きプロビットを推測しています（プロビットを知っている人には、以下のオーバーリーフハイパーリンクにも非常に短い説明が追加されています）。しかし、私の従属変数は、8パーセンテージグループに分類されたパーセンテージです。私が知っている意味カテゴリ1は0％を意味し、カテゴリ2は0<y<5などを意味する。結果として、すべての閾値はalphaであり、これをmy likelihood functionで使用することができる（このオーバーリーフリンクの式（2）参照）。誰かがRやStataでこれについてのコマンドを知っていますか、あるいはそのようなコマンドが存在していますか？

Answer 1

私はあなたがStataの中oprobitまたはoglmと、このようにそれを行うことができると思います。

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} < X'B +をu < = cut_ {j}）、これはあなたのビンに正確に収まらないかもしれません。カットポイントからc（mindouble）を追加/減算して、必要なものを得る必要があるかもしれません。

ここでは話題にはなりませんが、「latent」変数の範囲が未知のスケールであるという観測からは、序数的アプローチが正方形にならないと思います。

私はロバストネスチェックとしてintregを試してみます。範囲の問題は扱いませんが、スケールの問題はその設定で問題ではありません。