Я пытаюсь использовать пакет mstate и функцию coxph с помощью команды strata примерно для 10 000 наблюдений. Но каждый раз, когда я запускаю команду coxph, сеанс прерывается. Я действительно нуждаюсь в помощи. Вот код:
tmat <- transMat(x = list(c(2, 3, 4, 8), c(5,6,8), c(5,7,8), c(6,7,8),c(8),c(8),c(8),c()),
names = c("health", "IA", "IS", "DM", "IA+IS","IA+DM","IS+DM","death"))
cmdata <- read.csv("save_m_df.csv")
cmdata_30 <- subset(cmdata,age_5g=="5")
mscm_30 <- msprep(time = c(NA, "IA_date", "IS_date","DM_date", "IA_IS_date","IA_DM_date","IS_DM_date","death_date"),
status = c(NA, "IA", "IS", "DM","IA_IS","IA_DM","IS_DM","is_death"),
data = cmdata_30 , trans = tmat,id=cmdata_30$THJID)
mscm_30[, c("Tstart", "Tstop", "time")] <- mscm_30[, c("Tstart","Tstop", "time")]/365.25
c0_30 <- coxph(Surv(Tstart, Tstop, status) ~ strata(trans), data = mscm_30,method = "breslow")