[R] Reparamaterization for threshold values for ordinal regression

[treeyp]

I'm writing a function for ordinal regression for an ad hoc situation.  To 
estimate threshold parameters (4 levels), I applied reparameterization 
method to avoid using linear constraint suggested  by C. Croux. This is 
working fine, but I have a problem finding varaince for the  estimated  
parameters because it repamaterize as follows : 



r1 = b1, 

r2 = b1+b2^2 

r3=b1+b2^2+b3^2. 



It ensures ordering the constraints , but I don't know how to get variances 
for r2 and r3 . Is ther e any other method? 



Kyong
        [[alternative HTML version deleted]]


______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.