# HG changeset patch # User testtool # Date 1506948665 14400 # Node ID b4d7a110dd32ed895e956aab31eb427377986b5b # Parent d95c8e3b24e7aed0e60e36bdad919097d1e27078 Uploaded diff -r d95c8e3b24e7 -r b4d7a110dd32 survival/survival.R --- a/survival/survival.R Mon Oct 02 08:44:25 2017 -0400 +++ b/survival/survival.R Mon Oct 02 08:51:05 2017 -0400 @@ -4,12 +4,12 @@ input = args[1] type = args[2] error = args[3] -conf.type = args[4] -conf.lower = args[5] +conf_type = args[4] +conf_lower = args[5] fit = args[6] png = args[7] -tofit <- survival::survfit(Surv(obsT, status) ~ X1, type = type, error = error, conf.type = conf.type , conf.lower = conf.lower, data = read.csv(input)) +tofit <- survival::survfit(Surv(obsT, status) ~ X1, type = type, error = error, conf.type = conf_type , conf.lower = conf_lower, data = read.csv(input)) 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")])