annotate survival/survival.R @ 0:6cc0c5eb71b2 draft

Uploaded
author testtool
date Fri, 29 Sep 2017 18:11:27 -0400
parents
children 5b70a33bf327
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
1 args <- commandArgs(trailingOnly = TRUE)
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
2
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
3 input =args[1]
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
4 fit = args[2]
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
5 png = args[3]
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
6
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
7 #fit a Kaplan-Meier and plot it
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
8 input<-("test-data/input.csv")
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
9 fit <- survfit(Surv(obsT, status) ~ X1, type="kaplan-meier", conf.type="log", data=read.csv(input))
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
10
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
11 res <- summary(fit,times=c(1,2,3))
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
12 save.df <- as.data.frame(res[c("strata", "time", "n.risk", "n.event", "surv", "std.err", "lower", "upper")])
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
13 write.csv(save.df, file = "survfit.csv")
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
14
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
15 png("fit.png")
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
16 plot(fit)
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
17 dev.off()
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
18