TraMineR中seqcost时变替换成本使用报错求助

问题原因 #

当用seqcost(..., time.varying = T)生成时变成本时，返回的lifcost$indel是随时间变化的向量，但seqdist的method="OM"（最优匹配）默认仅接受单个数值的插入/删除成本，无法直接解析向量形式的时变indel，因此触发Error: 'indel' has an unknown value。

解决方案 #

要保留时变成本特性，需改用seqdist的CUSTOM方法，该方法专门支持时变或状态依赖的成本设置：

library(TraMineR)
data(biofam)

biofam.seq <- seqdef(biofam[501:600,10:25])

# 生成时变插入删除成本和替换矩阵
lifcost <- seqcost(biofam.seq, method="INDELSLOG", time.varying = T)

# 使用CUSTOM方法计算时变最优匹配距离
biofam.om <- seqdist(biofam.seq, method="CUSTOM", indel=lifcost$indel, sm=lifcost$sm)

备选方案（不推荐，丢失时变特性） #

如果不需要保留时变特性，可以将时变indel取均值作为单个固定值传入OM方法：

# 计算时变indel的平均值
avg_indel <- mean(lifcost$indel)
# 使用固定indel值计算OM距离
biofam.om <- seqdist(biofam.seq, method="OM", indel=avg_indel, sm=lifcost$sm)

内容的提问来源于stack exchange，提问作者giac