flexsurvreg `[.default`(Y, , "time") 中的错误：下标越界

Question

鉴于

Surv

函数的以下数据，并尝试使用

flexsurvreg

拟合模型，我收到一条错误消息：

[.default

(Y, , "time") 错误：下标越界

我不确定为什么会出现此错误。我必须与我的

Surv

对象

Y

中的数据有关，但我无法弄清楚它是什么，或者为什么会发生。

library(flexsurv)

Y <- structure(c(15.9474257914955, 92.8646984681524, 60.7937542591016, 
39.0043051102707, 22.2735472359904, 278.426492884738, 23.248652819273, 
4.17711769344757, 2.24005517260132, 274.773385211189, 21.3520961126228, 
18.2518223950669, 125.312959877672, 85.3612084472231, 27.2734023160286, 
36.3620638309367, 403.428793492735, 152.710509835273, 3.10392989813245, 
127.28582946954, 63.2099882741064, 44.5655964009001, 2.78577709845033, 
10.1090428767665, 102.198023729669, 343.099100197125, 17.2156423414208, 
29.4150785094056, 32.4848806916993, 3.6542577596261, 6.95172956945501, 
5.87113165766953, 56.8024939247263, 264.697690687394, 71.7590478868313, 
14.6108577317073, 20.8013409027592, 70.8457539651151, 25.791752475344, 
3.26064956790472, 57.7113479672798, 235.938105241419, 194.317991436403, 
42.4436470057455, 257.834274067571, 102.290239804092, 10.100812340795, 
34.5474927005247, 10.0337594860639, 383.575147185537, 391.74058499839, 
9.47984328773105, 7.35426482323835, 67.7660799235721, 10.5793123247722, 
389.128247057677, 58.5388316464529, 143.045014956796, 23.8206669979057, 
6.13273314050339, 40.376519704907, 1.66133588374672, 100.825597177946, 
11.5298575402338, 403.428793492735, 278.651209575017, 21.1528979520127, 
125.560296291348, 16.9371425236829, 8.58662139224698, 39.1666176547387, 
10.6516800452639, 54.6887503526312, 49.9744555673712, 8.29554464332079, 
7.00761139711887, 28.7522035577908, 168.782824659952, 1.31802200518933, 
25.852352792303, 102.306773590937, 250.435169558085, 14.4728494359721, 
32.5760886362187, 45.4089016526169, 11.0372942507258, 68.5019135755667, 
52.2303579024536, 403.428793492735, 403.428793492735, 141.801217464064, 
7.40242943834269, 29.100557117208, 8.96866215119505, 12.4354900525681, 
8.16703087179552, 58.1586028095045, 11.4000475582873, 25.4583340915452, 
6.16028494710057, 403.428793492735, 84.4824268536514, 50.7414776998092, 
18.8735292056578, 403.428793492735, 19.5885585047774, 11.4281913868452, 
160.288238443443, 11.4617504201645, 18.2467841931974, 6.20805769909054, 
5.66932843252971, 109.730417416886, 286.568917910547, 47.0260093039365, 
23.5178655461873, 403.428793492735, 40.8042534601449, 19.3234743656241, 
17.820853712635, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1
), .Dim = c(120L, 2L), .Dimnames = list(NULL, c("", "status")), type = "right", inputAttributes = list(
    time = list(dim = c(120L, 1L), dimnames = list(c("1", "2", 
    "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", 
    "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", 
    "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", 
    "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", 
    "44", "45", "46", "47", "48", "49", "50", "51", "52", "53", 
    "54", "55", "56", "57", "58", "59", "60", "61", "62", "63", 
    "64", "65", "66", "67", "68", "69", "70", "71", "72", "73", 
    "74", "75", "76", "77", "78", "79", "80", "81", "82", "83", 
    "84", "85", "86", "87", "88", "89", "90", "91", "92", "93", 
    "94", "95", "96", "97", "98", "99", "100", "101", "102", 
    "103", "104", "105", "106", "107", "108", "109", "110", "111", 
    "112", "113", "114", "115", "116", "117", "118", "119", "120"
    ), NULL))), class = "Surv")

flexsurvreg(
  Y ~ 1,
  dist = "lnorm"
  )

Answer 1

Surv

对象看起来不正确。

尝试显式命名

time

维度并重新运行回归：

attr(Y,'dimnames')[[2]][[1]] <- 'time'