Re: [R] aftreg vs survreg loglogistic aft model (different intercept term)

2010-11-29 Thread Terry Therneau
Survreg maximizes the log-likelihood to a relative tolerance of 1e-9 (?survreg.control). The printout shows -379503.5, to see the rest of the digits you need something like: fit - survreg( print(fit$loglik, digits=9) Aftreg printed even less digits; you would have to do the

Re: [R] aftreg vs survreg loglogistic aft model (different intercept term)

2010-11-26 Thread Immanuel
I add an example , all the variables are mutually excluding dummy variables, notice the different intercept: 5.627 vs 5.545: survreg: Value Std. Error zp (Intercept) 5.6270.00887 634.3 0.00e+00 Var1.recR2 -0.1080.01026 -10.5 1.00e-25 Var1.recR3 -0.490

[R] aftreg vs survreg loglogistic aft model (different intercept term)

2010-11-25 Thread Immanuel
Hi, I'm estimating a loglogistic aft (accelerated failure time) model, just a simple plain vanilla one (without time dependent covariates), I'm comparing the results that I obtain between aftreg (eha package) and survreg(surv package). If I don't use any covariate the results are identical , if I