Untitled

 avatar
unknown
plain_text
5 months ago
1.5 kB
7
Indexable

library(readxl)
D1 <- read_excel("D:\\Academic\\4H\\402\\Data-1.xlsx")

library(survival)
km.model <- survfit(Surv(Time,Died) ~ 1,
                    data = D1,
                    type = "kaplan-meier")
km.model
summary(km.model)

plot(km.model)
plot(km.model, conf.int = T, xlab = "Time(month)", ylab = "%Alive",
     main = "kM-Model", las =1, mark.time = T)

km.model2 <- survfit(Surv(Time, Died) ~ Over_40,
                     data = D1,
                     type = "kaplan-meier")
km.model2
summary(km.model2)

plot(km.model2)
plot(km.model2, conf.int = F, xlab = "Time(month)", ylab = "%Alive",
     main = "kM-Model", las = 1, mark.time = T)

survdiff(Surv(Time, Died) ~ Over_40,
         data = D1)

#2


library(readr)
D2 = read_csv("D:\\Academic\\4H\\402\\Data-2.csv")

D2$censored[D2$Turnover == 1] <- 1
D2$censored[D2$Turnover == 0 | D2$Turnover == 2] <- 0

library(survival)
km.fit1 <- survfit(Surv(LOS, censored)~1,
                   data = D2,
                   type = "kaplan-meier")

km.fit1
summary(km.fit1)
summary(km.fit1, times = c(30, 60, 90*(1:30)))

plot(km.fit1)
plot(km.fit1, conf.int = T, main = "KM")


km.fit2 <- survfit(Surv(LOS, censored) ~ Gender,
                   data = D2,
                   type = "kaplan-meier")

km.fit2
print(km.fit2)
summary(km.fit2)
summary(km.fit2, times = c(30, 60, 90*(1:30)))


plot(km.fit2)
plot(km.fit2, conf.int = T, main = "KM")

coxph <- coxph(Surv(LOS, censored) ~ Gender, data = D2 )
summary(coxph)


Editor is loading...
Leave a Comment