Mercurial > repos > testtool > survival

changeset 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 e31c47bf46a4
files survival/survival.R
diffstat 1 files changed, 5 insertions(+), 7 deletions(-) [+]
line wrap: on
line diff
--- a/survival/survival.R	Fri Sep 29 18:32:23 2017 -0400
+++ b/survival/survival.R	Fri Sep 29 18:41:35 2017 -0400
@@ -5,15 +5,13 @@
 fit = args[2]
 png = args[3]
 
-#fit a Kaplan-Meier and plot it
-input<-("test-data/input.csv")
-fit <- survival::survfit(Surv(obsT, status) ~ X1, type="kaplan-meier", conf.type="log", data=read.csv(input))
+tofit <- survival::survfit(Surv(obsT, status) ~ X1, type="kaplan-meier", conf.type="log", data=read.csv(input))
 
-res <- summary(fit,times=c(1,2,3))
+res <- summary(tofit,times=c(1,2,3))
 save.df <- as.data.frame(res[c("strata", "time", "n.risk", "n.event", "surv", "std.err", "lower", "upper")])
-write.csv(save.df, file = "survfit.csv")
+write.csv(save.df, fit)
 
-png("fit.png")
-plot(fit)
+png(png)
+plot(tofit)
 dev.off()