Displaying 2 results from an estimated 2 matches for "auric4".
Did you mean:
auric
2012 Mar 21
0
multivariate ordinal probit regression vglm()
...edictor (age) and
one or more response variables (the stages of skeletal morphology).
For example, the pubic symphysis and auriclar surface (two joint
surfaces of the pelvis) may be observed in three and four stages,
respectively (see sample dataframe "refdata" below).
age pube3 auric4
1 32 3 2
2 42 3 2
3 27 2 1
4 39 2 1
5 85 3 4
I've had some success in fitting the ordinal probit model using both
polr(method="probit") in the MASS library and vglm() in the VGAM
library, b...
2012 Apr 04
0
multivariate ordered probit regression---use standard bivariate normal distribution?
...data = refdata)
Coefficients:
(Intercept):1 (Intercept):2 age:1 age:2
-1.65895567 -2.14755951 0.06688242 0.04055919
Degrees of Freedom: 1492 Total; 1488 Residual
Residual Deviance: 1188.909
Log-likelihood: -594.4543
##########################################################
unvar.prob2<-vglm(auric4~age,data=refdata,cumulative(link="probit",parallel=FALSE,reverse=TRUE))
> unvar.prob2
Call:
vglm(formula = auric4 ~ age, family = cumulative(link = "probit",
parallel = FALSE, reverse = TRUE), data = refdata)
Coefficients:
(Intercept):1 (Intercept):2 (Intercept):3 age:1 age:...