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

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。我用TraMineR 说明使用mvad 数据传送

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)

希望这会有所帮助。