Procesamientos basicos. Analisis de correspondencia
Beatriz Valdez
02 julio 2015
Algunos ejemplos de análisis de correspondencia. Dado que la data no está totalmente limpia, el ejemplo es sólo para mostrar el código y sobre su interpretación.
Subimos la data que deseamos explorar. Los paquetes necesarios los subiremos en forma progresiva, en la medida en que los necesitemos.
read.csv("~/nvH.csv", header=TRUE)
require(memisc)
require(FactoMineR)Preparamos la data:
activos <-data.set(cComunal = nvH$CONSEJO.COMUNAL, nevera =nvH$NEVERA,
lavadora = nvH$NEVERA, comedor =nvH$COMEDOR,
muebles = nvH$MUEBLES,
utencilios = nvH$UTENCILIOS.D.COCINA,
ventilador = nvH$VENTILADOR, computador = nvH$COMPUTADOR,
tv = nvH$TV
)
activosDf <-as.data.frame(activos)Invocamos el análisis de correspondencia múltiple
mca1<-MCA(activosDf, graph = FALSE)
plot(mca1)
Imposible de entender. Veamos algunas alternativas:
plot(mca1,invisible=c("ind","quali.sup","quanti.sup"),cex=0.7)
Posiblemente debemos hacer algo con las no respuestas, en aras de poder mejorar la visualización. Queda claro sí que la primera dimensión tiene que ver con la tenencia exclusiva de utencilios del hogar, o la tenencia compartida; en tanto que la dimensión 2 se caracteriza por quienes no respondieron sobre la modalidad de tenencia y los que sí respondieron.
Del lado de los que tienen activos, observamos cuatro subgrupos: - los que comparten los activos, - los que no respondieron,
- los que tienen activos propios. - Los que no tienen activos.
Cerca de los que no responden o comparte activos, encontramos al Consejo comunal Nueva Esperanza. En este grupo, se comparten utencilios, lavadora, nevera, muebles, comedor, televisor y computador, o se dejó de responder sobre la tenencia.
El Consejo comunal La Pica está cercano al grupo de los que predominar la no tenencia de activos. El consejo comunal Ali Primera se encuentra cerca del grupo en los que predomina la tenencia no compartida de activos
Para facilitar entender las dimensiones, podemos invocar:
dimdesc(mca1)## $`Dim 1`
## $`Dim 1`$quali
## R2 p.value
## nevera 0.81749858 0.000000e+00
## lavadora 0.81749858 0.000000e+00
## utencilios 0.75423614 0.000000e+00
## muebles 0.61775415 1.355526e-240
## comedor 0.60209091 1.580421e-230
## tv 0.52582584 1.407342e-186
## ventilador 0.37855331 8.202218e-119
## computador 0.16612341 3.033039e-45
## cComunal 0.02901919 7.406223e-07
##
## $`Dim 1`$category
## Estimate p.value
## utencilios_Compartido 1.109631102 3.976474e-317
## lavadora_Compartida 0.986357450 0.000000e+00
## nevera_Compartida 0.986357450 0.000000e+00
## muebles_Compartido 1.143580591 1.892919e-215
## comedor_Compartido 1.091913091 5.360187e-188
## tv_Compartido 1.468715960 4.605045e-178
## ventilador_Compartido 1.371454207 2.965731e-111
## computador_Compartida 0.954849014 3.457204e-35
## comedor_NR 0.382470094 6.901885e-14
## lavadora_NR 0.803805471 3.465230e-13
## nevera_NR 0.803805471 3.465230e-13
## computador_NR 0.102413042 1.610635e-08
## muebles_NR 0.189708088 1.882366e-07
## utencilios_NR 0.307157415 2.103648e-07
## ventilador_NR 0.006985899 2.092459e-06
## Con la Unión la Nueva Esperanza 0.158793149 1.448027e-03
## utencilios_No tiene -0.652476008 4.592828e-02
## La Pica -0.218154535 2.182954e-02
## tv_NR -0.166172187 2.520590e-04
## ventilador_No tiene -0.696213427 7.485169e-05
## computador_No tiene -0.428280893 5.376832e-06
## lavadora_No tiene -0.889761362 4.736250e-06
## nevera_No tiene -0.889761362 4.736250e-06
## computador_Propia -0.628981164 1.968899e-06
## Ali Primera -0.134588568 1.290866e-06
## comedor_No tiene -0.650187167 5.519974e-09
## ventilador_Propio -0.682226679 2.118135e-11
## muebles_No tiene -0.603387107 6.080386e-13
## muebles_Propio -0.729901572 2.757654e-28
## comedor_Propio -0.824196018 8.931378e-32
## tv_Propio -0.673025720 8.083733e-68
## lavadora_Propia -0.900401560 3.702628e-121
## nevera_Propia -0.900401560 3.702628e-121
## utencilios_Propio -0.764312509 1.119911e-166
##
##
## $`Dim 2`
## $`Dim 2`$quali
## R2 p.value
## nevera 0.63274603 1.266195e-250
## lavadora 0.63274603 1.266195e-250
## muebles 0.55277549 3.045240e-201
## comedor 0.55171036 1.201781e-200
## ventilador 0.52745814 1.924017e-187
## utencilios 0.51882305 6.643589e-183
## computador 0.44581700 1.653033e-147
## tv 0.30731566 1.222709e-91
## cComunal 0.04734827 1.982001e-11
##
## $`Dim 2`$category
## Estimate p.value
## lavadora_NR 3.81281106 3.055198e-223
## nevera_NR 3.81281106 3.055198e-223
## muebles_NR 2.24897414 2.575746e-187
## comedor_NR 2.25125035 1.899634e-184
## ventilador_NR 2.26287114 2.367700e-183
## utencilios_NR 2.74076732 2.040013e-163
## computador_NR 1.47859774 3.632568e-148
## tv_NR 1.38779498 2.006676e-88
## Con la Unión la Nueva Esperanza 0.28413647 1.051990e-13
## tv_No tiene -0.44509087 8.796790e-03
## Nuevo Horizonte III -0.07526275 4.870944e-03
## Ali Primera -0.08433264 4.462510e-03
## computador_Compartida -0.64917502 3.360471e-03
## ventilador_Propio -0.61037751 1.277113e-03
## ventilador_No tiene -0.65998362 7.932980e-04
## ventilador_Compartido -0.99251002 1.731942e-06
## comedor_No tiene -0.68090657 4.744725e-07
## tv_Compartido -0.64206509 3.914981e-07
## tv_Propio -0.30063902 7.898577e-08
## muebles_No tiene -0.70740961 6.971939e-08
## muebles_Compartido -0.97580207 4.139352e-11
## comedor_Compartido -1.01724084 7.638652e-13
## utencilios_Compartido -1.16412411 5.669541e-17
## lavadora_Compartida -1.53204276 4.124027e-17
## nevera_Compartida -1.53204276 4.124027e-17
## computador_No tiene -0.46142655 4.793833e-29
##
##
## $`Dim 3`
## $`Dim 3`$quali
## R2 p.value
## nevera 0.7282375 0.000000e+00
## lavadora 0.7282375 0.000000e+00
## comedor 0.2902854 1.467890e-85
## muebles 0.2808646 2.930129e-82
## utencilios 0.2189478 1.342009e-61
## tv 0.1851034 5.377175e-51
## ventilador 0.1496260 2.359944e-40
## computador 0.1155458 1.488512e-30
## cComunal 0.0649354 5.763806e-16
##
## $`Dim 3`$category
## Estimate p.value
## lavadora_No tiene 0.88588692 2.147345e-303
## nevera_No tiene 0.88588692 2.147345e-303
## comedor_No tiene 0.32961605 3.960745e-81
## muebles_No tiene 0.30408971 2.025398e-75
## utencilios_No tiene 0.63454450 4.543050e-61
## tv_No tiene 0.63726203 3.588762e-51
## ventilador_No tiene 0.30728269 1.698787e-41
## computador_No tiene 0.27793984 2.692389e-27
## Ali Primera 0.07425228 2.070781e-07
## La Pica 0.26280226 9.781365e-05
## utencilios_NR 0.01601114 1.602533e-02
## lavadora_NR -0.00107350 2.514948e-02
## nevera_NR -0.00107350 2.514948e-02
## muebles_Compartido -0.10256163 1.639214e-02
## ventilador_Compartido -0.19440792 1.173617e-02
## comedor_Compartido -0.11254050 4.010708e-03
## computador_Compartida -0.18714799 1.186629e-05
## Nuevo Horizonte III -0.24908234 1.034816e-14
## computador_Propia -0.23285696 2.404396e-25
## utencilios_Propio -0.36008712 1.057951e-25
## tv_Propio -0.30500494 4.250135e-26
## ventilador_Propio -0.15622070 5.351555e-32
## muebles_Propio -0.32371330 1.035580e-68
## comedor_Propio -0.30358393 4.883084e-70
## lavadora_Propia -0.55984777 1.965209e-139
## nevera_Propia -0.55984777 1.965209e-139Estos indicadores no dan información adicional sobre el tipo de activos que determina cada dimensión. Asi, observamos que los activos más relevantes tanto para la dimensión 1 como para la dimesuón 2 son nevera ylavadora. Esta información nos permite entender el siguiente gráfico
plotellipses(mca1,keepvar=1:8)
Veamos el resultado empleando el método “Burt”
mca12<-MCA(activosDf, method ="Burt", graph = FALSE)
dimdesc(mca12)## $`Dim 1`
## $`Dim 1`$quali
## R2 p.value
## nevera 0.81749858 0.000000e+00
## lavadora 0.81749858 0.000000e+00
## utencilios 0.75423614 0.000000e+00
## muebles 0.61775415 1.355526e-240
## comedor 0.60209091 1.580421e-230
## tv 0.52582584 1.407342e-186
## ventilador 0.37855331 8.202218e-119
## computador 0.16612341 3.033039e-45
## cComunal 0.02901919 7.406223e-07
##
## $`Dim 1`$category
## Estimate p.value
## utencilios_Compartido 1.109631102 3.976474e-317
## lavadora_Compartida 0.986357450 0.000000e+00
## nevera_Compartida 0.986357450 0.000000e+00
## muebles_Compartido 1.143580591 1.892919e-215
## comedor_Compartido 1.091913091 5.360187e-188
## tv_Compartido 1.468715960 4.605045e-178
## ventilador_Compartido 1.371454207 2.965731e-111
## computador_Compartida 0.954849014 3.457204e-35
## comedor_NR 0.382470094 6.901885e-14
## lavadora_NR 0.803805471 3.465230e-13
## nevera_NR 0.803805471 3.465230e-13
## computador_NR 0.102413042 1.610635e-08
## muebles_NR 0.189708088 1.882366e-07
## utencilios_NR 0.307157415 2.103648e-07
## ventilador_NR 0.006985899 2.092459e-06
## Con la Unión la Nueva Esperanza 0.158793149 1.448027e-03
## utencilios_No tiene -0.652476008 4.592828e-02
## La Pica -0.218154535 2.182954e-02
## tv_NR -0.166172187 2.520590e-04
## ventilador_No tiene -0.696213427 7.485169e-05
## computador_No tiene -0.428280893 5.376832e-06
## lavadora_No tiene -0.889761362 4.736250e-06
## nevera_No tiene -0.889761362 4.736250e-06
## computador_Propia -0.628981164 1.968899e-06
## Ali Primera -0.134588568 1.290866e-06
## comedor_No tiene -0.650187167 5.519974e-09
## ventilador_Propio -0.682226679 2.118135e-11
## muebles_No tiene -0.603387107 6.080386e-13
## muebles_Propio -0.729901572 2.757654e-28
## comedor_Propio -0.824196018 8.931378e-32
## tv_Propio -0.673025720 8.083733e-68
## lavadora_Propia -0.900401560 3.702628e-121
## nevera_Propia -0.900401560 3.702628e-121
## utencilios_Propio -0.764312509 1.119911e-166
##
##
## $`Dim 2`
## $`Dim 2`$quali
## R2 p.value
## nevera 0.63274603 1.266195e-250
## lavadora 0.63274603 1.266195e-250
## muebles 0.55277549 3.045240e-201
## comedor 0.55171036 1.201781e-200
## ventilador 0.52745814 1.924017e-187
## utencilios 0.51882305 6.643589e-183
## computador 0.44581700 1.653033e-147
## tv 0.30731566 1.222709e-91
## cComunal 0.04734827 1.982001e-11
##
## $`Dim 2`$category
## Estimate p.value
## lavadora_NR 3.81281106 3.055198e-223
## nevera_NR 3.81281106 3.055198e-223
## muebles_NR 2.24897414 2.575746e-187
## comedor_NR 2.25125035 1.899634e-184
## ventilador_NR 2.26287114 2.367700e-183
## utencilios_NR 2.74076732 2.040013e-163
## computador_NR 1.47859774 3.632568e-148
## tv_NR 1.38779498 2.006676e-88
## Con la Unión la Nueva Esperanza 0.28413647 1.051990e-13
## tv_No tiene -0.44509087 8.796790e-03
## Nuevo Horizonte III -0.07526275 4.870944e-03
## Ali Primera -0.08433264 4.462510e-03
## computador_Compartida -0.64917502 3.360471e-03
## ventilador_Propio -0.61037751 1.277113e-03
## ventilador_No tiene -0.65998362 7.932980e-04
## ventilador_Compartido -0.99251002 1.731942e-06
## comedor_No tiene -0.68090657 4.744725e-07
## tv_Compartido -0.64206509 3.914981e-07
## tv_Propio -0.30063902 7.898577e-08
## muebles_No tiene -0.70740961 6.971939e-08
## muebles_Compartido -0.97580207 4.139352e-11
## comedor_Compartido -1.01724084 7.638652e-13
## utencilios_Compartido -1.16412411 5.669541e-17
## lavadora_Compartida -1.53204276 4.124027e-17
## nevera_Compartida -1.53204276 4.124027e-17
## computador_No tiene -0.46142655 4.793833e-29
##
##
## $`Dim 3`
## $`Dim 3`$quali
## R2 p.value
## nevera 0.7282375 0.000000e+00
## lavadora 0.7282375 0.000000e+00
## comedor 0.2902854 1.467890e-85
## muebles 0.2808646 2.930129e-82
## utencilios 0.2189478 1.342009e-61
## tv 0.1851034 5.377175e-51
## ventilador 0.1496260 2.359944e-40
## computador 0.1155458 1.488512e-30
## cComunal 0.0649354 5.763806e-16
##
## $`Dim 3`$category
## Estimate p.value
## lavadora_No tiene 0.88588692 2.147345e-303
## nevera_No tiene 0.88588692 2.147345e-303
## comedor_No tiene 0.32961605 3.960745e-81
## muebles_No tiene 0.30408971 2.025398e-75
## utencilios_No tiene 0.63454450 4.543050e-61
## tv_No tiene 0.63726203 3.588762e-51
## ventilador_No tiene 0.30728269 1.698787e-41
## computador_No tiene 0.27793984 2.692389e-27
## Ali Primera 0.07425228 2.070781e-07
## La Pica 0.26280226 9.781365e-05
## utencilios_NR 0.01601114 1.602533e-02
## lavadora_NR -0.00107350 2.514948e-02
## nevera_NR -0.00107350 2.514948e-02
## muebles_Compartido -0.10256163 1.639214e-02
## ventilador_Compartido -0.19440792 1.173617e-02
## comedor_Compartido -0.11254050 4.010708e-03
## computador_Compartida -0.18714799 1.186629e-05
## Nuevo Horizonte III -0.24908234 1.034816e-14
## computador_Propia -0.23285696 2.404396e-25
## utencilios_Propio -0.36008712 1.057951e-25
## tv_Propio -0.30500494 4.250135e-26
## ventilador_Propio -0.15622070 5.351555e-32
## muebles_Propio -0.32371330 1.035580e-68
## comedor_Propio -0.30358393 4.883084e-70
## lavadora_Propia -0.55984777 1.965209e-139
## nevera_Propia -0.55984777 1.965209e-139No se observan grandes diferencias entre los dos métodos