[R] Parametric survival models with left truncated, right censored data

Beth Ireland Beth.Ireland at hymans.co.uk
Wed Jan 23 16:46:31 CET 2008


Dear All,

I would like to fit some parametric survival models using left
truncated, right censored data in R. However I am having problems
finding a function to fit parametric survival models which can handle
left truncated data.

I have tested both the survreg function in package survival: 

fit1 <- survreg(Surv(start, stop, status) ~ X + Y + Z, data=data1)

and the psm function in package Design: 

fit2 <- psm(Surv(start, stop, status) ~ X + Y + Z, data=data1)


But neither function appears to work with left truncated data. The error
message "Invalid survival type" is received for both functions when left
truncated data is specified.

The function weibreg in the eha package fits Weibull survival models and
works with left truncated data. This function is useful, but I need to
fit parametric models other than Weibull.

Any suggestions of functions that can fit parametric survival models
other than Weibull on left truncated, right censored data would be
greatly appreciated.

Kind Regards,

Beth 
--------------------
Hymans Robertson LLP is a limited liability partnership registered in
England and Wales with registered number OC310282. A list of members 
of Hymans Robertson LLP is available for inspection at One London
Wall, London, EC2Y 5EA, the firm's registered office. Hymans
Robertson LLP is authorised and regulated by the Financial Services
Authority. Certain regulated services are provided by Hymans
Robertson Financial Services LLP, a subsidiary of Hymans Robertson
LLP. Hymans Robertson Financial Services LLP is a limited liability
partnership registered in England and Wales with registered number
OC310836. A list of members of Hymans Robertson Financial Services
LLP is available for inspection at One London Wall, London, EC2Y 5EA,
the firm's registered office.  Hymans Robertson Financial Services
LLP is authorised and regulated by the Financial Services Authority.
This e-mail and any attachments are confidential. If it ...{{dropped:9}}



More information about the R-help mailing list