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
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
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
3 matches
Mail list logo