Skip to content

Commit d4274f2

Browse files
committed
Fix wrong assumption that time > 0 in ppc_km_overlay()
1 parent d6e2135 commit d4274f2

File tree

1 file changed

+6
-7
lines changed

1 file changed

+6
-7
lines changed

R/ppc-censoring.R

Lines changed: 6 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -91,12 +91,6 @@ ppc_km_overlay <- function(
9191

9292
data <- ppc_data(y, yrep, group = status_y)
9393

94-
if (!is.null(left_truncation_y)) {
95-
data$left_trunc <- left_truncation_y[data$y_id]
96-
} else {
97-
data$left_trunc <- 0
98-
}
99-
10094
# Modify the status indicator:
10195
# * For the observed data ("y"), convert the status indicator back to
10296
# a numeric.
@@ -109,7 +103,12 @@ ppc_km_overlay <- function(
109103
as.numeric(as.character(.data$group)),
110104
1))
111105

112-
sf_form <- survival::Surv(time = data$left_trunc, time2 = data$value, event = data$group) ~ rep_label
106+
if (is.null(left_truncation_y)) {
107+
sf_form <- survival::Surv(time = data$value, event = data$group) ~ rep_label
108+
} else {
109+
sf_form <- survival::Surv(time = left_truncation_y[data$y_id], time2 = data$value, event = data$group) ~ rep_label
110+
}
111+
113112
if (!is.null(add_group)) {
114113
data <- dplyr::inner_join(data,
115114
tibble::tibble(y_id = seq_along(y),

0 commit comments

Comments
 (0)