使用 TraMineR 的 seqedplot：如何在一张图中绘制选定的危险曲线？答案

【问题标题】：seqedplot using TraMineR: how to plot selected hazard curve in one single graph?使用 TraMineR 的 seqedplot：如何在一张图中绘制选定的危险曲线？
【发布时间】：2014-04-24 18:48:20
【问题描述】：

我正在使用 TraMineRextras，特别是我对 seqedplot 函数感兴趣。

这里有一个例子：

data(actcal.tse)
actacl.tse <- actcal.tse[1:500,]
iseq <- unique(actcal.tse$id)
nseq <- length(iseq)
data(actcal)
actcal <- actcal[rownames(actcal) %in% iseq,]
actcal.seqe <- seqecreate(actcal.tse)
seqelength(actcal.seqe) <- rep(12, nseq)

而不是通过以下方式绘制所有曲线：

seqedplot(actcal.seqe, type="hazard", breaks=6, group=actcal$sex, lwd=3)

我只会使用ignore 选项根据变量sex 绘制选定的曲线：

seqedplot(actcal.seqe, type="survival", group=actcal$sex, lwd=3, ignore= c("NoActivity", "Start", "FullTime", "Stop", "LowPartTime", "Increase", "Decrease"))

如何在一个图表中而不是在两个单独的图表中绘制男性和女性选择的曲线？

提前谢谢你

【问题讨论】：

标签： r traminer

【解决方案1】：

一种方法可能是将性别绑定到事件类型：

actcal.tse$event <- sprintf("%s-%s", actcal.tse$event, actcal$sex[match(actcal.tse$id, rownames(actcal))])

因为似乎没有使用groups 参数在图表中拆分的选项。所以尝试这些方法，并完全省略组选项

【讨论】：

谢谢，我会按照您的建议尝试一下。
这是个好主意。然而，如果不调整每种性别的风险人群，则由此产生的生存概率将是错误的。所以恐怕你需要忘记使用seqedplot按性别叠加曲线。
谢谢吉尔伯特，我会在找到其他解决方案后立即更新帖子