Williams Scott
2006-Feb-26 03:02 UTC
[R] frailty in coxph or repeated measures in cph (Design)
I am trying to build a model to aid a clinical decision. Certain patients have a blood marker measured at each visit - a rise of this may indicate recurrence of the cancer after treatment (endpoint is "clinical recurrence", censored). In a proportion (up to 30%), this rise is a false positive - hence I wish to correlate factors at the time of the rising test to clinical recurrence, preferably expressed using nomogram (Design library). Many patients have more than one rise and sometimes even 5 rises doesn't mean inevitable failure. The aim is to identify the true positives early, to offer further therapy. The data form presently disregards the potential multiplicity of measurements per individual, treating each as independent. For example: id <- c(a,a,a,b,b,b) # patient id risenumber <- c(1,2,3,1,2,3) # number of each rise per patient clinfail <- c(1,1,1,0,0,0) # censored indicator of clinical failure status (the endpoint) clinfailtime <- c(5,4,3,10,9,8) # time from rise to clinical failure endpoint riseval <- c(10,20,30,1,2,3) # value of test at the time of rise timesinceRx <- c(1,2,3,1,2,3) # years since treatment Rxtype <- c(c,c,c,d,d,d) # type of treatment ...plus other variables at time of rise plus pre-treatment variables In generic terms, analysis would be: fit <- cph(Surv(clinfailtime, clinfail) ~ riseval +Rxtype...) nomogram(fit) I could easily convert these data to counting process notation and use time-dependent covariates aligned to treatment date. I do think the question is different however - not "what predicts clinical failure from the time of treatment" - rather, "what predicts the risk of clinical failure beyond the time of rise". To complicate things, the different treatments will be likely to have non-proportional hazards, with Rxtype=='d' likely to have large numbers of false positive rises at 1-2 years due to differing biology. The obvious lack of independence between the successive rises in an individual is the problem - my feeling is to use the current data format with: fit <- coxph(Surv(clinfailtime, clinfail) ~ riseval +Rxtype +frailty(id) The frailty option is not available in cph (Design), meaning an illustrative nomogram would be more difficult to produce. Any thoughts on what might be the best solution to this would be most welcome. Scott Williams MD Radiation Oncology Peter MacCallum Cancer Centre Melbourne, Australia
Reasonably Related Threads
- Nomogram after Cox Random Effect (frailty) model
- Approaches of Frailty estimation: coxme vs coxph(...frailty(id, dist='gauss'))
- coxph models with frailty
- frailty models in survreg() -- survival package (PR#2933)
- Frailty by strata interactions in coxph (or coxme)?