[R] Extract p-value from aftreg object

David Winsemius dwinsemius at comcast.net
Thu Feb 18 19:16:03 CET 2010


On Feb 18, 2010, at 12:33 PM, David Winsemius wrote:

> I don't seem to have the unnamed package loaded that has aftreg(),  
> but in general you ought to be able to get what you want by looking  
> not just at the aftreg-object but also at the print(aftreg)-object  
> using str().

Trying that approach after loading package eha shows it will not work,  
and looking at the print method shows that no object is returned.

-- 
David.
>
>
> On Feb 18, 2010, at 11:07 AM, Philipp Rappold wrote:
>
>> 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
> -- 
> David Winsemius, MD
> Heritage Laboratories
> West Hartford, CT
>
> ______________________________________________
> R-help at 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.

David Winsemius, MD
Heritage Laboratories
West Hartford, CT



More information about the R-help mailing list