给出平均时间图的准确值 (seqmtplot)

【问题标题】：Give the exact values of the mean time plots (seqmtplot)给出平均时间图的准确值 (seqmtplot)
【发布时间】：2014-11-10 03:46:05
【问题描述】：

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

这里有一个示例脚本，它使用 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?

提前致谢，约翰内斯

【问题讨论】：

标签： r plot sequence traminer

【解决方案1】：

可以使用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)

【讨论】：