使用TraMineR中的组从seqrep中提取单个代表性序列图-Java 学习之路

我正在尝试使用 TraMineR 中的 group= 选项重现 seqrep 所做的图表，因为我想要一个图像文件，其中包含每个组的代表性序列，以便在演示文稿中使用 . 似乎我在语法中缺少一个重要的部分，因为序列对象的子集产生了与 TraMineR 中的板载 group 选项非常不同的结果 . 正如您在下面的图片中看到的那样，由 seqrep 识别的代表性序列强烈不同 . 你知道我的语法中哪些部分缺失或错误了吗？

library(TraMineR)
library(WeightedCluster)

data(mvad)
mvad.alphabet <- c("employment", "FE", "HE", "joblessness", "school",
                   "training")
mvad.labels <- c("Employment", "Further Education", "Higher Education",
                 "Joblessness", "School", "Training")
mvad.scodes <- c("EM", "FE", "HE", "JL", "SC", "TR")

## Define sequence objects
mvad.seq <- seqdef(mvad[, 17:86], alphabet = mvad.alphabet, 
                   states = mvad.scodes, labels = mvad.labels, 
                   weights = mvad$weight, xtstep = 6)

## Defining cost matrices for OM - custom and transition based
subm.trate <- seqsubm(mvad.seq, method="TRATE")

## Computing dissimilarity matrices
mvad.ham <- seqdist(mvad.seq, method="HAM", sm=subm.trate)

## Clustering OM dissimilarities
wardCluster <- hclust(as.dist(mvad.ham), method = "ward",
                      members = mvad$weight)
clust4 <- cutree(wardCluster, k = 4)

## Plotting representative sequences
seqrplot(mvad.seq, criterion="dist", group=clust4, dist.matrix=mvad.ham, 
         withlegend=FALSE, border=NA, weighted=TRUE, stats=FALSE)
seqrplot(mvad.seq[clust4==2,], criterion="dist", dist.matrix=mvad.ham, 
         withlegend=FALSE, border=NA, weighted=TRUE, stats=FALSE)

Four groups of representative sequences: see second group

Second group of representative sequences

1 回答

2
在对状态序列对象进行子集化时，应相应地对不相似矩阵进行子集化
```
seqrplot(mvad.seq[clust4==2,], criterion="dist", 
         dist.matrix=mvad.ham[clust4==2, clust4==2], 
         withlegend=FALSE, border=NA, weighted=TRUE, stats=FALSE)
```
如果这样做， seqrep 将使用错误的距离来确定代表性序列 .
回复于 2024-05-03T09:13:19+08:00

使用TraMineR中的组从seqrep中提取单个代表性序列图

1 回答

相关问题