Суммировать атрибуты последовательностей в одном объекте последовательности?

Я использую TraMineR для анализа наборов последовательностей. Каждый связанный набор последовательностей может содержать 100 рабочих процессов из одного проекта за один период времени. Используя TraMineR, я могу легко вычислить описательную статистику для каждой последовательности, однако меня больше интересует описательная статистика самого объекта последовательности, включающая в себя все более мелкие последовательности.

Например, чтобы получить частоты состояний, я запускаю:

seqstatd(sequences.sts)

Однако это дает мне частоты состояний для каждой последовательности в моем объекте последовательности. Я хочу получить доступ к частотам состояний во всех последовательностях внутри моего объекта последовательности. Как я могу это сделать?


person histelheim    schedule 04.11.2013    source источник


Ответы (1)


Я не уверен, что понимаю ваш вопрос, поскольку seqstatd() возвращает частоты поперечного сечения в каждой последующей позиции, а НЕ частоты состояний для каждой последовательности. Последний возвращается seqistatd().

Предполагая, что вы ссылаетесь на результат seqistatd(), вы получите среднее время, проведенное в каждом состоянии с seqmeant(sequence.sts).

Для других сводок вы можете использовать функцию apply. Например, вы получаете дисперсию времени, проведенного в каждом состоянии с

tab <- seqistatd(mvad.seq)
vart <- apply(tab,2,var)
head(vart)

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

person Gilbert    schedule 04.11.2013