在R中使用msSurv软件包
我试图用msSurv来解决一个多状态建模问题,这个问题着眼于个体向不同阶段的过渡。其中一部分是创建一个树对象,我认为这是一个错误,但我不明白它是什么。我将在这里包括最低可行的示例在R中使用msSurv软件包,r,statistics,survival-analysis,R,Statistics,Survival Analysis,我试图用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"
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)
可以找到有问题的数据集
非常感谢您的帮助。我可能有误解,但您的6组没有1-6作为优势,因此程序返回错误,因为实际上您说6与计算无关。关于解决方案,我认为6应该有边,因为在这一行中可能需要有边:
“6”=list(edges=NULL))
此外,您应该包括用于此问题的所有包,并且不要发布到数据集的链接。使用dput(数据)
-复制并粘贴输出到您的问题。
No states eligible for exit distribution calculation.
Entry distributions calculated for states 6 .
Error in bs.IA[, , j, b] : subscript out of bounds