Hi,
I am trying to simulate a regression on survival data under a few
conditions:
1. Under different error distributions
2. Have the error term be dependent on the covariates

But I'm not sure how to specify either conditions. I am using the Design
package to perform the survival analysis using the survreg, bj, coxph
functions.   Any help is greatly appreciated.

This is what I have so far:
survtime  <- 10*rexp(500)  #distribution of survival time
cens <- ifelse(survtime > 10, 0, 1) #indicator for censored/observed
survtime  <- pmin(survtime, 10) #new survival time values with censored
info
age <- rnorm(200, 40, 10)  #age variable
race <- factor(sample(c('a','b'),500,TRUE))  #categorical variable
test <- bj(Surv(survtime, cens) ~ rcs(age,5) + race)

	[[alternative HTML version deleted]]

