goranbrostrom/eha

Small bug in the print.phreg function?

Opened this issue · 0 comments

First I would like to say that I really appreciate your book, both because it gives me an unusual perspective on the subject and because the synergy with the R package has great educational value.

During my exploration I noted that this code,

library(eha)
fit <-phreg(Surv(enter,exit,event)~1,data=mort)
fit
Call:
phreg(formula = Surv(enter, exit, event) ~ 1, data = mort)

Covariate          W.mean      Coef Exp(Coef)  se(Coef)    Wald p
log(scale)                    3.785               0.066     0.000 
log(shape)                    0.332               0.058     0.000 

Events                    276 
Total time at risk         17038 
Max. log. likelihood      -1399.3 

I derived from pag. 240, didn't actually show the Exp(Coef), although it was announced in the header.

I thought it might also be a wanted behaviour, in that it doesn't refer to a coefficient of a term included in the right side of the equation. However, the book on pag. 241 reports the R output showing the exponentiated value too, albeit in a transposed format, so I believe it could not be working as intended.

            log(scale) log(shape)
Coef             3.785      0.332
se(Coef)         0.066      0.058
exp(Coef)       44.015      1.393

Thank you for sharing your knowledge and your experience in the field.

Paolo

> R.version
               _                                
platform       x86_64-w64-mingw32               
arch           x86_64                           
os             mingw32                          
crt            ucrt                             
system         x86_64, mingw32                  
status                                          
major          4                                
minor          2.2                              
year           2022                             
month          10                               
day            31                               
svn rev        83211                            
language       R                                
version.string R version 4.2.2 (2022-10-31 ucrt)
nickname       Innocent and Trusting            

> packageVersion("eha")
[1] ‘2.10.1’