Befehl, um Schwellenwerte von geordneten Probit zu verwenden

Question

Ich schätze ein bestelltes Probit (für diejenigen, die nur probit kennen, fügte ich auch eine sehr kurze Erklärung in den umseitigen Hyperlink unten hinzu). Meine abhängige Variable ist jedoch ein Prozentsatz, der in acht Prozent-Gruppen eingeteilt wurde. Ich weiß z.B. diese Kategorie 1 bedeutet 0 Prozent, diese Kategorie 2 bedeutet 0<y<5 usw. Folglich kenne ich alle Schwellenwerte alpha und könnte sie in my likelihood function verwenden (siehe Gleichung (2) in diesem umseitigen Link). Kennt jemand einen Befehl dafür in R oder Stata oder existiert ein solcher Befehl überhaupt?

Answer 1

Ich denke, man es so mit oprobit oder oglm in Stata tun kann:

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 verwendet Pr (y = j | x) = Pr (cut_ {j-1} < x'b + u < = cut_ {j}), also passt das nicht genau zu Ihren Bins. Möglicherweise müssen Sie c (mindouble) vom Schnittpunkt addieren/subtrahieren, um zu erhalten, was Sie wollen.

Dies ist off-topic hier, aber ich denke, die ordinale Ansatz nicht mit der Beobachtung, dass Ihre "latente" Variable hat eine begrenzte Reichweite mit einer unbekannten Skala.

Ich würde intreg als Robustheitsprüfung versuchen. Es behandelt nicht das Bereichsproblem, aber das Skalierungsproblem ist in dieser Einstellung kein Problem.