对于使用软件包 relsurv 观察到的生存率，我确实有一个关于计算相对生存率的问题。
发布生殖示例：
为此，我们需要软件包relsurv library("relsurv")
d1 <- subset(colrec, site == "colon" & diag >= as.date("1Jan1994") &diag <= as.date("31Dec1995"))
d1$d.int <- 1
d2 <- subset(colrec, site == "colon" & diag >= as.date("1Jan1999") &diag <= as.date("31Dec2000"))
d2$d.int <- 2
d <- rbind(d1, d2)
ind <- which(d$time > 365.241 * 10)
d$time[ind] <- 365.241 * 10
d$stat[ind] <- 0
slopop[c("50", "70"), c("1990", "2000"), ]
从HMD下载文件http://www.mortality.org/。国家斯洛文尼亚 mltper_1x1.txt和fltper_1x1.txt
的年度生命周期表文件 slotab <- transrate.hmd(male = "mltper_1x1.txt",
+ female = "fltper_1x1.txt")
通过以下方式导入您自己的数据：
rmap = list(age = agey * 365.241, sex = sex, year = diag)
然后最后：
fit_rsr <- rs.surv(Surv(time, stat) ~ d.int,data = d, ratetable = slopop, method = "ederer1",add.times = c(5, 10) * 365.241,rmap = list(age = age, sex = sex, year = diag))
但是在最后一个代码中，我得到了这个错误： 评估中的错误（predvars，data，env）：找不到对象'year'。
这到底是什么意思？我该如何调整？
预先感谢，如果我需要添加一些内容，请告诉我。
Kweet