Произведите частоты состояния и индекс энтропии для конкретной переменной

Я могу генерировать отдельные графики из моего набора данных (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]), но не смог создать его для отдельных переменных.


person user3355411    schedule 28.05.2015    source источник


Ответы (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)

Надеюсь это поможет.

person Gilbert    schedule 18.06.2015