Я пытаюсь проверить гипотезу о том, что нет различий в выживаемости между четырьмя группами в данных (почки), сообщенных о случаях смерти 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)