Использование пакета msSurv в R - PullRequest
3 голосов
/ 30 апреля 2019

Я пытаюсь использовать msSurv для задачи моделирования с несколькими состояниями, которая рассматривает переход отдельных людей на разные стадии. Частично это создает объект дерева, в котором, как мне кажется, я совершаю ошибку, но не могу понять, что это такое. Я включу здесь минимально работоспособный пример.

Nodes <- c("1", "2", "3", "4", "5", "6")
Edges <- list("1" = list(edges = c("2", "3", "4", "5", "6")),
"2" = list(edges = c("1", "3", "4", "5", "6")),
"3" = list(edges = c("1", "2", "4", "5", "6")),
"4" = list(edges = c("1", "2", "3", "5", "6")),
"5" = list(edges = c("3", "4", "6")),
"6" = list(edges = NULL))
treeobj <- new("graphNEL", nodes = Nodes, edgeL = Edges, edgemode = "directed")
fit3 <- msSurv(df, treeobj, bs = TRUE, LT = TRUE)

Я получаю следующую ошибку.

No states eligible for exit distribution calculation.
Entry distributions calculated for states 6 .
Error in bs.IA[, , j, b] : subscript out of bounds

Данный набор данных можно найти здесь .

Любая помощь искренне приветствуется.

1 Ответ

0 голосов
/ 09 мая 2019

Возможно, я неправильно понимаю, но ваша 6 группа не имеет 1-6 в качестве ребра, поэтому программа возвращает ошибку, потому что, по сути, вы говорите, что 6 не связано с вычислением.Что касается решения, я считаю, что у 6 должны быть ребра, так как в этой строке, возможно, должны быть ребра: "6" = list(edges = NULL))

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...