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.
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’
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,
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.
Thank you for sharing your knowledge and your experience in the field.
Paolo