require(survival)
set.seed(1)
S <- Surv(runif(100))
x <- runif(100)
dd <- datadist(x); options(datadist='dd')
f <- psm(S ~ x, dist="exponential")
summary(f) # effects on log(T) scale
f.ph <- pphsm(f)
if (FALSE) summary(f.ph) # effects on hazard ratio scale
options(datadist=NULL)
Run the code above in your browser using DataLab