es.davy.ai

Preguntas y respuestas de programación confiables

¿Tienes una pregunta?

Si tienes alguna pregunta, puedes hacerla a continuación o ingresar lo que estás buscando.

Análisis de secuencias y estadísticas descriptivas de clústeres en R.

06 Aug, 2023 Programación

Actualmente estoy realizando un análisis de secuencia utilizando el paquete TraMineR en R. Sin embargo, tengo dificultades para encontrar cómo extraer estadísticas descriptivas para cada grupo que obtengo.
Utilizando el conjunto de datos mvad

mvad.seq <- seqdef(mvad, 17:86, alphabet = mvad.alphabet, states = mvad.scodes, 
    labels = mvad.labels, xtstep = 6)
clusterward1 <- agnes(dist.om1, diss = TRUE, method = "ward")
plot(clusterward1, which.plot = 2)
cl1.4 <- cutree(clusterward1, k = 4)
cl1.4fac <- factor(cl1.4, labels = paste("Tipo", 1:4))

¿Cómo puedo obtener información sobre cuántos hombres hay en cada grupo, cuántos en cada grupo son católicos, etc?

Tags:  , , ,

Answer

  1. Avatar for davy.ai
    davy.ai

    Para extraer estadísticas descriptivas para cada grupo, puedes utilizar la función seqdss() del paquete TraMineR. Aquí tienes un ejemplo de código para extraer el número de hombres y el número de individuos con el estado católico para cada grupo (asumiendo que ya se han creado el conjunto de datos “mvad” y el grupo “clusterward1”):

    # Crear un data frame con las etiquetas de los grupos
clusters <- data.frame(cluster = cl1.4fac)

# Crear una matriz de disimilitud para los datos de secuencia
diss <- seqdist(mvad.seq, method = "OM", indel = 1.5, with.missing = TRUE)

# Calcular las estadísticas descriptivas para cada grupo
cluster.stats <- seqdss(mvad.seq, diss = diss, group = clusters)

# Extraer el número de hombres y el número de individuos con el estado católico para cada grupo
n.males <- cluster.stats[c("M", "Meaningful Gower index"), ,"Type 1":"Type 4"]
n.catholic <- cluster.stats[c("Catholic", "Meaningful Gower index"), ,"Type 1":"Type 4"]

# Imprimir los resultados
print(n.males)
print(n.catholic)

    En este ejemplo, se utiliza seqdss() para calcular las estadísticas descriptivas para cada grupo. Especificamos group = clusters para indicar que queremos que las estadísticas se calculen para cada grupo definido por cl1.4fac.

    El objeto resultante cluster.stats es un objeto de seqdss que contiene varias matrices con las estadísticas descriptivas para cada grupo.

    Luego podemos extraer el número de hombres y el número de individuos con el estado católico para cada grupo utilizando los argumentos c("M", "Meaningful Gower index") y c("Catholic", "Meaningful Gower index") para especificar qué variables nos interesan. Los objetos resultantes n.males y n.catholic son matrices con una fila para cada variable y una columna para cada grupo.

    Finalmente, imprimimos los resultados en la consola utilizando la función print(). Esto mostrará el número de hombres y el número de individuos con el estado católico para cada grupo.

Comments are closed.