问题
there is quite a lot of information (internet and textbooks) on how to do survival analysis in R with the survival package. But I don't find any information on how to do this when you have left censored data.
Problem background:
I have a self constructed data set with published survival data. Usually the event time and the date of the last follow-up (right censoring) is given. There is however one study that only states that the event happened before day 360. So I left censored this data.
What I want to do:
I want to analyse the complete data set with left truncation, events, and right truncation. I want to plot the Kaplan-Meier curve by gender and then
- do a log-rank test
- do a Cox regression
What I need:
I am able to create a Surv object with type = interval2. But this does neither allow to calculate survdiff, nor coxph of the survival package.
The intcox package was removed from CRAN and I don't find what I search in the icenReg or interval packages.
Can anyone please give me a hind how to solve my problem or where to find practical information on this? I am already spending days on this one.
Many thanks!
回答1:
You can fit a Cox-PH model with both right and left censoring in icenReg by using the ic_sp function. You can fit this using the standard Surv response variable, i.e.
fit <- ic_sp(Surv(L, R, type = 'interval2') ~ treatment, data = myData)
or a little more succinctly with
fit <- ic_sp(cbind(L, R) ~ treatment, data = myData)
Log-rank tests are not available in icenReg, but can be found in the interval package.
来源：https://stackoverflow.com/questions/42069425/interval-censored-data-cox-proportional-hazard-and-surival-difference-in-r