[R] multivariate ordered probit regression---use standard bivariate normal distribution?

2012-04-03 Thread Trey Batey
mulative(link = "probit", parallel = FALSE, reverse = TRUE), data = refdata) Coefficients: (Intercept):1 (Intercept):2 (Intercept):3 age:1 age:2 -2.07719235 -2.43422370 -2.99123098 0.07319632 0.05133132 age:3 0.03797696 Degrees of Freedom: 2238 Total; 2232 Residual Residual Deviance: 1583.4

[R] multivariate ordinal probit regression vglm()

2012-03-21 Thread Trey Batey
cs.google.com/open?id=0B5zZGW2utJN0TEctcW1oblFRcTJrNDVLOVBmRWRaQQ). Thanks in advance. --Trey *** Trey Batey---Anthropology Instructor Mt. Hood Community College 26000 SE Stark St. Gresham, OR 97030 trey.ba...@mhcc.edu or ekt.ba...@gmail.com ## unrestricted cumula

[R] plotting survival curves (multiple curves on single graph)

2011-07-05 Thread Trey Batey
Hello. This is a follow-up to a question I posted last week. With some previous suggestions from the R-help community, I have been able to plot survival (, hazard, and density) curves using published data for Siler hazard parameters from a number of ethnographic populations. Can the function belo

[R] plotting survival curves with model parameters

2011-06-28 Thread Trey Batey
Hello. I am trying to write an R function to plot the survival function (and associated hazard and density) for a Siler competing hazards model. This model is similar to the Gompertz-Makeham, with the addition of a juvenile component that includes two parameters---one that describes the initial in