计算TraMiner中每个序列的转换率

2 回答

• 0

  来自 TraMineRextras 的 TSE_to_STS 和来自 TraMineR 的 seqtrate 的函数用于一组序列，不适用于单个序列 . 这是因为在内部它们使用的函数用于不适用于向量的表 .

  解决方法是添加具有虚拟事件的虚拟序列，并从结果的概率转换矩阵中删除虚拟事件 .

  test.events <- c("A","B","C","X")
  test.stm <- seqe2stm(test.events, dropList=list("A"=test.events[-1],
       B=test.events[-2], C=test.events[-3], X=test.events[-4]))
  test.tse <- data.frame(id = c(99,1,1,1), time = c(0,1,2,3), 
       event = c("X","A","B","C"))
  test.sts <- TSE_to_STS(test.tse, id = "id", timestamp = "time", 
       event = "event", stm=test.stm, tmin=1, tmax=4, firstState="None")
  test.seqdef <- seqdef(test.sts,informat = "STS")
  test.trate <- seqtrate(test.seqdef)
  test.trate[-nrow(test.trate),-ncol(test.trate)]
  

  希望这可以帮助 .

  回复于 2024-05-19T12:11:00+08:00
• 1

  基于吉尔伯特的解释，这是我修改过的代码 . 它创建一个具有不同id（= 99）的相同序列 . 利用两个序列的转换率相同，转换矩阵与用一个序列计算的相同 . 它无需创建虚拟事件即可工作 .

  test.events <- c("A","B","C")
  test.stm <- seqe2stm(test.events, dropList=list("A"=test.events[-1], B=test.events[-2], C=test.events[-3]))
  test.tse <- data.frame(id = c(1,1,1), time = c(1,2,3), event = c("A","B","C"))
  test.tse.bis <- test.tse
  test.tse.bis[,1] <- 99
  test.tse <- rbind(test.tse,test.tse.bis)
  test.sts <- TSE_to_STS(test.tse, id = "id", timestamp = "time", event = "event", stm=test.stm, tmin=1, tmax=4, firstState="None")
  test.seqdef <- seqdef(test.sts,informat = "STS")
  seqtrate(test.seqdef)
  

  回复于 2024-05-19T12:11:00+08:00