R() without object failed for some forms of interval censoring

git-svn-id: svn://svn.r-forge.r-project.org/svnroot/ctm@2611 edb9625f-4e0d-4859-8d74-9fd3b1da38cb

Author: thothorn

pkg/mlt/R/R.R

@@ -219,6 +219,8 @@ R.numeric <- function(object = NA, cleft = NA, cright = NA,
 
     ### treat extremely small intervals as `exact' observations
     d <- cright - cleft
+    if (length(object) == 1L && all(is.na(object)))
+        object <- rep(NA, length(d))
     if (any(!is.na(d) | is.finite(d))) {
         if (any(d < 0, na.rm = TRUE)) stop("cleft > cright")
         i <- (d < tol)

pkg/mlt/tests/bugfixes.R

@@ -134,3 +134,9 @@ m <- mlt(ctm(response = Bernstein_basis(numeric_var("y", bounds = c(0, 1)), orde
 mi <- mlt(ctm(response = Bernstein_basis(numeric_var("y", bounds = c(0, 1)), order = 1)),
               data = d[-ic,,drop = FALSE], theta = c(-1, 1))
 all.equal(logLik(m), logLik(mi))
+
+### R with interval censoring
+lwr <- as.double(1:3)
+upr <- as.double(2:4)
+upr[2] <- 2
+R(cleft = lwr, cright = upr)

pkg/mlt/tests/bugfixes.Rout.save

@@ -1,6 +1,6 @@
 
-R version 4.4.2 (2024-10-31) -- "Pile of Leaves"
-Copyright (C) 2024 The R Foundation for Statistical Computing
+R version 4.4.3 (2025-02-28) -- "Trophy Case"
+Copyright (C) 2025 The R Foundation for Statistical Computing
 Platform: x86_64-pc-linux-gnu
 
 R is free software and comes with ABSOLUTELY NO WARRANTY.
@@ -170,6 +170,13 @@ Coefficients: 30.951 38.885 42.523 -0.01508 -0.058906 -0.12136 -0.054061 -0.4769
 > all.equal(logLik(m), logLik(mi))
 [1] TRUE
 > 
+> ### R with interval censoring
+> lwr <- as.double(1:3)
+> upr <- as.double(2:4)
+> upr[2] <- 2
+> R(cleft = lwr, cright = upr)
+[1] ( 1,  2]  2       ( 3,  4]
+> 
 > proc.time()
    user  system elapsed 
-  0.933   0.068   0.997 
+  0.816   0.068   0.879