diff --git a/coxph.R b/coxph.R index 3ccf4fb..5393f31 100644 --- a/coxph.R +++ b/coxph.R @@ -44,7 +44,7 @@ sma_cox_res <- function(m) { # sma_coxph_univ <- function() { univ_formulas <- sapply(covariates, function(x) as.formula( - paste('Surv(time_at_risk, time_at_risk_status) ~ ', + paste('Surv(time_at_risk, ifelse(time_at_risk_status == 1, 1, 0)) ~ ', x))) univ_models <- lapply(univ_formulas, function(x) {coxph(x, data = secmal)})