Дать точные значения средних временных графиков (seqmtplot)

Question

Я хотел бы получить точные значения моих временных графиков времени для каждого различного состояния, которое я определил для дальнейшего анализа моих кластеров. Средние графики времени не позволяют мне легко прочитать эти значения . Есть ли способ извлечь эти значения?

Вот вам пример сценария с использованием набора данных 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? 

Спасибо заранее, Johannes

Answer 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)