2

我想为我定义的每个不同状态提取平均时间图的确切值，以进一步分析我的集群。平均时间图不允许我轻松读取这些值。有没有办法提取这些值？

这里有一个示例脚本，它使用 TraMineR 中可用的 mvad 数据集，R 中可用。

library(TraMineR) 
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)
## seqmtPlot
seqmtplot(mvad.seq, group=mvad$gcse5eq, withlegend=TRUE, border=NA, xtstep=3, sortv="from.start") #Example of seqmtplot. How can I extract the mean time values from these plots?

在此先感谢，约翰内斯

4

1 回答 1

2

可以使用该seqmeant函数获得该值。

seqmeant(mvad.seq)

但是，此函数没有group参数。您需要为每个组分别计算此信息。对于“否”组，您可以使用

seqmeant(mvad.seq[mvad$gcse5eq=="no", ])

如果要为每个组计算它，可以使用以下by函数：

by(mvad.seq, mvad$gcse5eq, seqmeant)

希望这可以帮助！

PS：您使用的某些参数seqmtplot在这种情况下没有意义，您应该使用

seqmtplot(mvad.seq, group=mvad$gcse5eq, withlegend=TRUE)
于 2014-11-07T08:13:40.820 回答