Mercurial > repos > testtool > survival

annotate survival/survival.R @ 2:ce6df2bd21cd draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author testtool
date Fri, 29 Sep 2017 18:41:35 -0400
parents 5b70a33bf327
children 744cf0f4d1b1
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1
5b70a33bf327 Uploaded
testtool
parents: 0
diff changeset
1 require(survival, quietly = TRUE)
0
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
2 args <- commandArgs(trailingOnly = TRUE)
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
3
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
4 input =args[1]
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
5 fit = args[2]
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
6 png = args[3]
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
7
2
ce6df2bd21cd Uploaded
testtool
parents: 1
diff changeset
8 tofit <- survival::survfit(Surv(obsT, status) ~ X1, type="kaplan-meier", conf.type="log", data=read.csv(input))
0
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
9
2
ce6df2bd21cd Uploaded
testtool
parents: 1
diff changeset
10 res <- summary(tofit,times=c(1,2,3))
0
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
11 save.df <- as.data.frame(res[c("strata", "time", "n.risk", "n.event", "surv", "std.err", "lower", "upper")])
2
ce6df2bd21cd Uploaded
testtool
parents: 1
diff changeset
12 write.csv(save.df, fit)
0
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
13
2
ce6df2bd21cd Uploaded
testtool
parents: 1
diff changeset
14 png(png)
ce6df2bd21cd Uploaded
testtool
parents: 1
diff changeset
15 plot(tofit)
0
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
16 dev.off()
6cc0c5eb71b2 Uploaded
testtool
parents:
diff changeset
17