R 检验各组之间生存率没有差异的假设

我试图验证以下假设：863例肾移植患者死亡时间的数据（肾脏）显示，四组患者的存活率没有差异。患者可按种族和性别分为四组

我正在使用代码：

library(KMsurv)
data(kidney)
kidney
head(kidney)
str(kidney)

group<- rep(NA, dim(kidney)[1])
group[ kidney$gender==1 & kidney$race==1 ] <- 1
group[ kidney$gender==1 & kidney$race==2 ] <- 2
group[ kidney$gender==2 & kidney$race==1 ] <- 3
group[ kidney$gender==2 & kidney$race==2 ] <- 4

library(survival)
survdiff(Surv(obstime, death)~group, data=kidney)

survdiff(Surv(obstime, death)~race, data=kidney, subset=gender==1)

survdiff(Surv(obstime, death)~race, data=kidney, subset=gender==2)

survdiff(Surv(obstime, death)~race+strata(gender), data=kidney)

库（KMsurv）
数据（肾脏）
肾
头（肾）
str（肾脏）
组根据
？KMsurv:：RENDY
或
dim（RENDY）
，加载
数据（RENDY）
的RENDY数据框有119行和3列。如果您定义了自己的
renty
变量，请确保在加载肾脏数据后调用
data（renty）
不会覆盖该变量。您不需要创建
group
，您可以执行
survdiff（Surv（obstime，death）~gender+race，data=renty）
，但如果您是，我建议将其作为变量添加到数据框中，而不是单独的对象
getwd()
kidney<- read.csv("kidneytransplant.csv")
head(Kidney)
str(Kidney)