search for: unvar

...poor job asking the question. So, here is the general question: how can I run a run a multivariate (more than one non-independent, response variables) ordered probit regression model? I've had success doing this in the univariate case using the vglm() function in the VGAM package. For example: unvar.prob<-vglm(y~x,cumulative(link="probit",parallel=FALSE,r everse=TRUE)) This would fit a unrestricted cumulative probit model (or a "thresholds of change" model) with the ordered response variable treated as a latent trait modeled with a standard normal distribution. What I...