Question

我正在使用TraMineR，并且对R和TraMineR还是陌生的。

实际上，我使用TraMineR和R中的簇库对生命历程数据集进行了分类。（使用本指南：http://traminer.unige.ch/preview-typology.shtml）

现在，我将不同的案例分类为不同的类型（所有4种类型）。

我想对某种类型进行更深入的分析，但是我需要知道哪些案例（我有案例编号）属于哪种类型。＃是否有可能将案例分类的特定类型作为新变量写入数据集本身？还有另一种方法吗？

Answer 1

在参考指南的示例中，使用最佳匹配距离并根据过渡概率基于替换成本，按如下方式获得类型

library(TraMineR)
data(mvad)
mvad.seq <- seqdef(mvad, 17:86)
dist.om1 <- seqdist(mvad.seq, method = "OM", indel = 1, sm = "TRATE")

library(cluster)
clusterward1 <- agnes(dist.om1, diss = TRUE, method = "ward")
cl1.4 <- cutree(clusterward1, k = 4)

cl.4是一个向量，其序列的簇成员资格与mvad数据集相对应。（将它转换成一个因子可能很方便。）因此，您可以简单地将此变量作为附加列添加到数据集中。例如，如果我们想将此新列命名为Type

mvad$Type <- cl1.4
tail(mvad[,c("id","Type")])  ## id and Type of the last 6 sequences

##     id Type
## 707 707    3
## 708 708    3
## 709 709    4
## 710 710    2
## 711 711    1
## 712 712    4

可以将类型学写入数据集吗？

1 个答案: