[R] Extract p-value from aftreg object

Philipp Rappold philipp.rappold at gmail.com
Thu Feb 18 17:07:13 CET 2010


Dear all,

does anyone know how I can extract specific p-values for covariates 
from an aftreg object? After fitting a model with aftreg I can find 
all different variables by using str(), but there's no place where 
p-values are kept. The odd thing is that print() displays them 
correctly.

EXAMPLE:

 > testdata
   start stop censor groupvar       var1      var2
1     0    1      0        1 0.91663902 0.0847912
2     1    2      0        1 0.60470753 0.6487798
3     2    3      0        1 0.09599891 0.2195178
4     3    4      1        1 0.86384189 0.6667897
5     0    1      0        2 0.07747445 0.8782836
6     1    2      0        2 0.44608030 0.2218685
7     2    3      1        2 0.77317152 0.3813840


 > fit1 <- aftreg(Surv(start, stop, censor)~var1, data=testdata)


 > fit1
Call:
aftreg(formula = Surv(start, stop, censor) ~ var1, data = testdata)

Covariate          W.mean      Coef Exp(Coef)  se(Coef)    Wald p
var1                0.540     0.150     1.162     0.770     0.845

log(scale)                    1.358     3.890     0.260     0.000
log(shape)                    2.015     7.502     0.781     0.010

Events                    2
Total time at risk             7
Max. log. likelihood      -1.4026
LR test statistic         0.05
Degrees of freedom        1
Overall p-value           0.816436
	WALD P IS DISPLAYED CORRECTLY.


Any help is highly appreciated, I'm going nuts here ;)

Thanks
Philipp



More information about the R-help mailing list