Hello, I have a question regarding the use of an offset term with survreg(), in the Survival library. In particular, I am trying to figure out on what scale the offset term should be.
Here's a simple example with no censoring and no coefficients: --------- y = rlnorm(1000, meanlog = 10, sdlog = 2) delta = rep(1, 1000) int = rep(1, 1000) survreg(Surv(y,delta)~offset(10*int), dist = "lognormal") Call: survreg(formula = Surv(y, delta) ~ offset(10 * int), dist = "lognormal") Coefficients: (Intercept) -0.1317181 Scale= 1.985558 Loglik(model)= -11973.1 Loglik(intercept only)= -11973.1 n= 1000 --------- I would have expected the intercept estimate to be 10, since I specified offset = 10, unless the offset term should be specified on another scale. Thank you, Ming ______________________________________________ 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.
