Как с помощью TraMineR определить количество переходов в отдельных последовательностях для каждого типа перехода?

В TraMineR функция seqtransn позволяет вычислить количество переходов в отдельных последовательностях.

Мне было интересно: с TraMineR, есть ли способ определить количество переходов в отдельных последовательностях, но для каждого типа перехода?

Например, если у меня есть последовательности состояний с тремя различными состояниями, у меня будет 6 возможных переходов. И я хотел бы подсчитать для каждого человека количество появлений каждого из 6 возможных переходов.

Заранее спасибо за вашу помощь,

Лучший,


person Scoub'    schedule 03.04.2018    source источник


Ответы (1)


Я покажу вам здесь первые 3 последовательности данных biofam, как вы можете получить количество различных возможных переходов для каждой последовательности:

data(biofam)
biofam.seq <- seqdef(biofam[1:3,10:25])

## First, make each sequence an element of a list    
biofam.seq.list <- setNames(split(biofam.seq, seq(nrow(biofam.seq))), rownames(biofam.seq))

## and apply seqtrate to each element of the list
tr.count <- lapply(biofam.seq.list, seqtrate, count=TRUE)
tr.count

## $`1167`
##        [-> 0] [-> 1] [-> 3] [-> 6]
## [0 ->]      8      0      1      0
## [1 ->]      0      0      0      0
## [3 ->]      0      0      0      1
## [6 ->]      0      0      0      5
##
## $`514`
##        [-> 0] [-> 1] [-> 3] [-> 6]
## [0 ->]      0      1      0      0
## [1 ->]      0      9      1      0
## [3 ->]      0      0      0      1
## [6 ->]      0      0      0      3
##
## $`1013`
##        [-> 0] [-> 1] [-> 3] [-> 6]
## [0 ->]      6      1      0      0
## [1 ->]      0      4      1      0
## [3 ->]      0      0      0      1
## [6 ->]      0      0      0      2
person Gilbert    schedule 03.04.2018
comment
Привет. Хорошее решение. Спасибо за помощь! Лучший. - person Scoub'; 05.04.2018
comment
Спасибо за этот комментарий, но посетите страницу справки Что мне делать, когда кто-то ответит на мой вопрос?. - person Gilbert; 06.04.2018