产生特定变量的状态频率和熵指数

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)

希望这有帮助。