[R] heteroskedastic ordered probit

From: Rendas-Baum, Regina <Regina_Rendas-Baum_at_brown.edu>
Date: Sat 08 Jul 2006 - 13:12:37 EST

Dear all,  

is there existing code to fit a ordinal regression model with probit link where the variance can be modeled? For example, the polr function from the MASS package fits:  

probitP(Y <= k | x) = zeta_k - eta

Instead of,  

probitP(Y <= k | x) = (zeta_k - eta)/scale, where scale = exp(b*z)  

which is what I need.  

I've seen these models go by the name of heteroskedastic ordered probit, location-scale glms, generalized nonlinear models, ... . I know the package gnlm (Jim Lindsey) has a function to fit this model using the logit link, so I'm wondering if this is already done in R using probit instead.  



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Sat Jul 08 13:20:34 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sat 08 Jul 2006 - 20:15:24 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.