Bug fix for aft with left truncation

ellessenne · Aug 9, 2024 · 26ed06a · 26ed06a
1 parent 4c6ddf5
commit 26ed06a
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 1 deletion.
diff --git a/R/aft.R b/R/aft.R
@@ -411,7 +411,7 @@ aft <- function(formula, data, smooth.formula = NULL, df = 3,
         warning("Design matrix for the acceleration factor is not full rank")
     X <- X[, X.index, drop=FALSE]
     XD0 <- X0 <- XD <- matrix(0,1,ncol(X))
-    X_list = list()
+    X_list <- X_list0 <- list()
     gauss = gauss.quad(control$nNodes)
     if (delayed && all(time0==0)) delayed <- FALSE # CAREFUL HERE: delayed redefined
     if (tvc.integrated) {

diff --git a/inst/working_code.R b/inst/working_code.R
@@ -18,6 +18,15 @@
 ##   require(bbmle)
 ## }
 
+## Bug report: aft with left truncation
+library(rstpm2)
+brcancer$start <- 0
+fit <- aft(Surv(start, rectime,censrec==1)~hormon,data=brcancer,df=4) # now okay:)
+brcancer$start <- 1
+fit <- aft(Surv(start, rectime,censrec==1)~hormon,data=brcancer,df=4) # now okay:)
+
+length(rstpm2:::lhs(Surv(start, rectime, censrec == 1) ~ hormon))
+
 ## test for predict(..., type="lpmatrixD")
 library(rstpm2)
 fit = aft(Surv(rectime,censrec==1)~hormon,data=brcancer)