如何使用 R 包 MSTATE 计算每个非吸收状态的平均逗留时间

Question

我正在进行生存分析，但似乎不知道该怎么做。

以下是来自 MSTATE 教程的代码块，用于简单的 Cox 回归。
如何计算每种非吸收状态下的平均逗留时间？

代码：

library(mstate)
data(ebmt3)
tmat <- trans.illdeath(names=c("Tx","PR","RelDeath"))
ebmt3$prtime <- ebmt3$prtime/365.25
ebmt3$rfstime <- ebmt3$rfstime/365.25
covs <- c("dissub", "age", "drmatch", "tcd", "prtime")
msbmt <- msprep(time = c(NA, "prtime", "rfstime"), status = c(NA, "prstat", "rfsstat"), data = ebmt3, trans = tmat, keep = covs)
expcovs <- expand.covs(msbmt, covs[2:3], append = FALSE)
msbmt <- expand.covs(msbmt, covs, append = TRUE, longnames = FALSE)
c1 <- coxph(Surv(Tstart, Tstop, status) ~ dissub1.1 + dissub2.1 +
 age1.1 + age2.1 + drmatch.1 + tcd.1 + dissub1.2 + dissub2.2 +
 age1.2 + age2.2 + drmatch.2 + tcd.2 + dissub1.3 + dissub2.3 +
 age1.3 + age2.3 + drmatch.3 + tcd.3 + strata(trans), data = msbmt,
 method = "breslow")

newd <- data.frame(dissub = rep(0, 3), age = rep(0, 3), drmatch = rep(0,
 3), tcd = rep(0, 3), trans = 1:3)
 newd$dissub <- factor(newd$dissub, levels = 0:2, labels = levels(ebmt3$dissub))
 newd$age <- factor(newd$age, levels = 0:2, labels = levels(ebmt3$age))
 newd$drmatch <- factor(newd$drmatch, levels = 0:1, labels = levels(ebmt3$drmatch))
 newd$tcd <- factor(newd$tcd, levels = 0:1, labels = levels(ebmt3$tcd))
 attr(newd, "trans") <- tmat
 class(newd) <- c("msdata", "data.frame")
 newd <- expand.covs(newd, covs[1:4], longnames = FALSE)
 newd$strata = 1:3
 newd
 msf1 <- msfit(c1, newdata = newd, trans = tmat)

谢谢！

Answer 1

我认为您正在寻找 mstate 中的 ELOS 函数 - 它代表一个州的预期停留时间 - 要完成您的示例，您需要使用 probtrans 计算转换概率，然后您可以计算 ELOS每个州。

pt <- probtrans(msf1,predt=0)
# ELOS until last observed time point
ELOS(pt)

如何使用 R 包 MSTATE 计算每个非吸收状态的平均逗留时间

How to calculate mean sojourn time in each nonabsorbing state using R package MSTATE

r

markov-chains

survival-analysis