cutp returns an empty named list

[lbraglia]

Hi Christopher

is there something wrong here?

df <- structure(list(time = c(69, 2493, 2469, 113, 2099, 328, 183, 
2289, 2302, 2110, 2206, 1059, 2273, 880, 404, 2196, 2039, 324, 
1811, 853, 2160, 2133, 758, 2036, 197, 90, 2093, 2083, 1939, 
2041, 2057, 381, 945, 2081, 2182, 2009, 2064, 1275, 2045, 353, 
2032, 1961, 761, 1917, 1436, 1138, 1757, NA, 1096, 1832, 1708, 
1794, 1917, 1905, 247, 1892, 75, 121, 1766, 1964, 849, 1796, 
1826, 1827, 1709, 338, 1814, 1708, 1639, 1616, 1476, 1569, 1482, 
667, 1642, 1464, 1607, 1690, 1592, 1562, 154, 1681, 1713, 1490, 
1689, 462, 1559, 1520, 736, 1459, 1441, 95, 616, 1722, 1580, 
1555, 1617, 1407, 1433, 553, 40, 1310, 1231, 1245, 1321, 1298, 
1272, 1255, 225, 1127, 1283, 1289, 1279, 128, 166, 1324, 1211, 
601, 1310, 86, 1098, 1152, 993, 101, 1058, 1205, 1058, NA, 1155, 
1253, 398, 1208, 482, 1095, 1117, 1167, 1070, 1143, 998, 504, 
1020, 1068, 917, 364, 1154, 278, 343, 354, 413, 196, 460, 106, 
158, 9, 880, 384, 271, 995, 473, 1005, 321, 612, 207, 44, 790, 
1034, 869, 900, 823, 945, 291, 417, 901, 554, 249, 763, 22, 753, 
193, 754), status = c(1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 
0L, 1L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 
0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 
1L, 0L, 1L, 0L, 0L, NA, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 
0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 
0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 
0L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 
0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 
0L, 1L, 0L, 0L, 0L, NA, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 
0L, 1L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 
0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 
1L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L), marker = c(-1.8, -2.5, 
-2.8, -2.5, -0.2, -0.5, -1, -0.4, -2.7, -1.6, -1.8, -6.2, -2.5, 
-0.7, -0.5, -0.6, -2.8, -3.2, -0.7, -4, -0.3, -4, -2.6, -1.4, 
-1, -2.5, -3.9, -0.3, -1, -1, -1.2, -1, -0.7, -3, -0.9, -0.5, 
-0.1, -1.1, -1.1, -2.8, -0.7, -0.2, -2, -1.5, -3.5, -7, -3.2, 
-1, -1, -17.9, -1.9, -1.1, -0.8, -2, -1.4, -1.1, -0.5, -2.5, 
-2.2, -0.1, -0.6, -0.2, -1.5, -0.3, -6, -3.8, -4.2, -1.7, -3.9, 
-2.5, -5.5, -2.6, -4.3, -3, -0.7, -0.5, -0.3, -1, -0.5, -1.7, 
-5, -3.1, -1.4, -0.4, -0.1, -3, -0.7, -0.3, -3.7, -1, -3.5, -10.5, 
-4.2, -2.2, -0.4, -4.5, -1, -2.8, -0.6, -4.7, -1, -4.9, -0.3, 
-0.6, -6.6, -3.6, -1.5, -0.9, -1.2, -1.4, -2.1, -1, -1.6, -3.2, 
-1, -3.4, -0.6, -0.3, -2.5, -1, -3.2, -3.8, -2.5, -1.5, -2.8, 
-0.4, -4, -0.6, -2, -0.5, -1.9, -2.8, -3.6, -3.8, -0.4, -2.4, 
-1.4, -1.8, -4.6, -3, -10.1, -3.9, -2, -1.3, -1.5, -0.2, -1, 
-0.3, -3, -0.7, -1.5, -1, -0.5, -2, -3.5, -0.7, -0.6, -0.3, -0.8, 
-3.9, -0.6, -2.9, -0.8, -0.8, -4, -1, -1.1, -2.6, -0.7, -2.7, 
-2.5, -18.8, -0.4, -7, -7.1, -4, -7, -4.2, -1, -0.5)), class = "data.frame", 
row.names = c(NA, -180L))

model <- survival::coxph(Surv(time = time, event = status) ~ marker, data = df)
survMisc::cutp(model)

to me it returns an empty named list. If it's not a bug, how to interpret it (not found in the man)?

thanks a lot for the function and the package, best, Luca

[lbraglia]

survMisc_0.4.tar.gz returns,

 > survMisc::cutp(model, var = 'marker')
    CP      Q      p
1 -1.5 0.6969 0.7164

so i guess it's a 0.5 issue