traminer - 为特定变量生成状态频率和熵指数

Question

我能够从我的数据集 ( DISDARAE) 中为不同的变量 ( GENDER, RACE) 生成单独的图，例如

seqIplot(DISDATAE.seq, border = NA, group = DISDATAE$GENDER, sortv = "from.start")
seqIplot(DISDATAE.seq, border = NA, group = DISDATAE$RACE, sortv = "from.start")

如何为每个变量生成单独的状态频率和熵表？我对整个数据集使用了这种语法：seqstatd(DISDATAE.seq[, 1:4])，但无法为单独的变量创建一个

Answer 1

只需使用by. 我说明了使用mvad随TraMineR

library(TraMineR)
data(mvad)
# creating the state sequence object
mvad.seq <- seqdef(mvad[, 15:86])

## Distributions and cross-sectional entropies by sex
by(mvad.seq, mvad$male, seqstatd)

希望这可以帮助。