Формат SPS репрезентативной последовательности из seqrep

Кто-нибудь знает, как извлечь репрезентативные последовательности из вывода seqrep в формате SPS? Это было бы очень полезно для зрителей сюжета из seqrplot.


person Yolande    schedule 12.03.2018    source источник


Ответы (1)


Вы можете распечатать вывод с аргументом format = 'SPS'. Я иллюстрирую biofam данными

data(biofam)
biofam.lab <- c("Parent", "Left", "Married", "Left+Marr",
                "Child", "Left+Child", "Left+Marr+Child", "Divorced")
biofam.seq <- seqdef(biofam, 10:25, labels=biofam.lab)

## Computing the distance matrix
costs <- seqcost(biofam.seq, method="INDELSLOG")
biofam.om <- seqdist(biofam.seq, method="OM", sm=costs$sm, indel=costs$indel)

## Representative set using the neighborhood density criterion
biofam.rep <- seqrep(biofam.seq, diss=biofam.om, criterion="density")
print(biofam.rep, format="SPS")

##     Sequence         
## [1] (0,16)           
## [2] (0,9)-(3,1)-(6,6)
## [3] (0,6)-(1,10)     

И если вы хотите получить репрезентативные последовательности в форме SPS, вы можете использовать seqformat и seqconc:

biofam.rep.sps <- seqconc(seqformat(biofam.rep, to='SPS'))

Результатом является матрица с одним столбцом, в которой каждая репрезентативная последовательность хранится в виде строки символов.

person Gilbert    schedule 13.03.2018
comment
Можно ли иметь представителя по кластеру, используя группу = фактор кластера? - person Yolande; 14.03.2018
comment
Это другой вопрос, и вам лучше всего создать для него отдельную ветку. Но вы можете получить представителей по кластеру, используя функцию seqrep.grp для TraMineRextras. - person Gilbert; 14.03.2018
comment
@ Гилберт. Благодарю вас! - person Yolande; 15.03.2018