2

クラスターをさらに分析するために定義したさまざまな状態ごとに、平均時間プロットの正確な値を抽出したいと考えています。平均時間プロットでは、これらの値を簡単に読み取ることができません。これらの値を抽出する方法はありますか?

ここに、R で利用可能な TraMineR で利用可能な mvad データセットを使用したサンプル スクリプトがあります。

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 に答える