Seqedplot с использованием TraMineR: как построить выбранную кривую опасности на одном графике?

Question

Я использую 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)

Я хотел бы только построить выбранные кривые в соответствии с переменной sex используя ignore опция:

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

Как я могу построить выбранные мужские и женские кривые на одном графике, а не на двух отдельных?

заранее спасибо

Answer 1

Одним из способов может быть привязка пола к типу события:

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

так как нет возможности использовать groups параметр для разделения внутри графика. Так что попробуйте что-нибудь в этом духе и полностью опустите опцию групп