[R] ESTIMATES SURVIVAL ANALYSIS WEIBULL

CHIRIBOGA Xavier xavier.chiriboga at unine.ch
Fri Apr 10 19:23:05 CEST 2015


Dear members,



I did a survival analysis, Weibul fits better for my data: this is the "head" of my data



 start stop state soil volatile hours replicate
1     0   48     1   ca       nc    48         3
2     0   48     1   ca       nc    48         4
3     0   48     1   ca       nc    48         5
4     0   48     1   ca       nc    48         8
5     0   60     1   ca        c    60         1
6     0   60     1   ca        c    60         1





I have 4 soils: ca, cs, cl and sand



and I have 2 volatiles: nc and c



However, when I get the stimates my soil "ca" does not appear and the volatile "c" neither



Call:
survreg(formula = Surv(hours, state) ~ soil * volatile, data = surgal,
    dist = "weibull")
                      Value Std. Error      z         p
(Intercept)          4.3605     0.0417 104.51  0.00e+00
soilcl              -0.0784     0.0590  -1.33  1.84e-01
soilcs              -0.3865     0.0567  -6.81  9.59e-12
soilsand            -0.1913     0.0579  -3.30  9.61e-04
volatilenc          -0.1600     0.0579  -2.76  5.74e-03
soilcl:volatilenc    0.1320     0.0823   1.60  1.09e-01
soilcs:volatilenc    0.2679     0.0812   3.30  9.64e-04
soilsand:volatilenc  0.0821     0.0809   1.01  3.10e-01
Log(scale)          -1.4788     0.0442 -33.46 1.73e-245

Scale= 0.228

Weibull distribution
Loglik(model)= -1036.3   Loglik(intercept only)= -1065
        Chisq= 57.41 on 7 degrees of freedom, p= 5e-10
Number of Newton-Raphson Iterations: 5
n= 257



Anyone knows why? or what does it mean? how should I interpret the values in the table?



Thanks for your help,



Xavier



More information about the R-help mailing list