r - 检验组间生存率没有差异的假设

Question

我正在尝试检验数据（肾脏）中四组之间的存活率没有差异的假设，该数据报告了 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)

但是我没有得到863例肾移植患者的数据，我只得到了119个观察值。我使用的代码有错误吗？我无法获得任何结果，因为我无法获得全部数据。

我试过这段代码来获取数据，但我仍然遇到同样的问题：

getwd()
kidney<- read.csv("kidneytransplant.csv")
head(Kidney)
str(Kidney)

标签： rsurvival