Cómo producir una secuencia de trama de frecuencia para un solo grupo dentro de una solución de clúster

Question

No podía encontrar una respuesta suficiente a mi problema, tal vez alguien pueda ayudar aquí?(Soy un principiante a R)

Hago análisis de la secuencia, el espacio de estado es n = 10 y el espacio de tiempo es t = 168 (meses).Dibujé una secuencia de trama de frecuencia para una solución de clúster con 8 grupos.Sin embargo, la trama no es muy abierto a interpretaciones, porque las parcelas son demasiado formalizado, o demasiado pequeño resp.(ver gráfico a continuación)

Hice los siguientes procedimientos hasta ahora (muy cerca de las instrucciones en el TraMineR-Ayuda-documento):

dist.om1 <- seqdist(neu.seq, method = "OM", indel = 1, sm = submat)
clusterward1 <- agnes(dist.om1, diss = TRUE, method = "ward")
cluster8 <- cutree(clusterward1, k = 8)
cluster8 <- factor(cluster8, labels = c("Typ 1", "Typ 2", "Typ 3", "Typ 4", "Typ 5", "Typ 6", "Typ 7", "Typ 8"))
seqfplot(neu.seq, group = cluster8, pbarw = T, withlegend = T)

Traté de volver a configurar los márgenes, pero el resultado era siempre el mismo diagrama (el adjunto de la parcela fue hecho con la configuración predeterminada).Eso pensaba yo, en su lugar, tal vez yo podría dibujar la secuencia de trama de frecuencia para un solo grupo dentro de mi 8-cluster-solución.(en Stata-código, me gustaría escribir algo como para un único índice de secuencia de la trama sqindexplot if cluster8 == 4) Sin embargo, no sé cómo se hace en R.Si alguien tiene una idea de cómo obtener una bonita secuencia de la trama de frecuencia, estaría muy agradecido!Gracias!Oliver

Answer 1

Con 8 grupos a los que usted puede necesitar para reducir el tamaño de fuente de las etiquetas de los ejes mediante el cex.plot argumento.Por ejemplo:

seqfplot(neu.seq, group = cluster8, withlegend = T, cex.plot=.5)

También puede obtener un mejor aspecto de las parcelas con el border=NA argumento que suprime el borde negro alrededor de las barras que representan a cada secuencia de patrones.

Alternativamente, si usted está utilizando dispositivos gráficos tales como pdf, png o jpeg para crear su parcela de archivos, trate de jugar con los parámetros width y height una de las funciones.El más grande de la height valor, menor es el texto que se ve fuera.

Para obtener sólo clúster 4, uso

seqfplot(neu.seq[cluster8=="Typ 4",], withlegend = T)

(Ver también Cómo identificar las secuencias dentro de cada grupo? )

Y si quieres combinar las parcelas por sí mismo, usando por ejemplo par(mfrow=c(.,.)) tienes que desactivar la corrección automática de la leyenda, e insertar la leyenda manualmente con seqlegend, por ejemplo,

par(mfrow=c(2,2))
seqfplot(neu.seq[cluster8=="Typ 4",], withlegend = F)
seqfplot(neu.seq[cluster8=="Typ 5",], withlegend = F)
seqfplot(neu.seq[cluster8=="Typ 6",], withlegend = F)
seqlegend(neu.seq)
dev.off()

Espero que esto ayude.