[R-es] AYUDA CON ERROR CON LA LIBRERIA PCA

[Camilo Calle]

Cordial saludo adjunto la base de datos y el script:

eu60<- read.csv(file.choose(), header=T, sep=";", dec=".", row.names=1)
eu60.pca <- PCA(eu60, quali.sup=19)
eu60.data <- cbind.data.frame(eu60[,19], eu60.pca$ind$coord)
eu60.ellipse <- coord.ellipse(eu60.data, bary=TRUE)
plot.PCA(ellipse=eu60.ellipse, cex=0.8)



El 4 de diciembre de 2013 17:35, Camilo Calle <camicalle en gmail.com>escribió:

> Muchas gracias por la rapidez y por la ayuda tan efectiva, me quedó el
> gráfico de maravilla, al parecer el problema era que habían individuos en
> una sola categoría, así quedó la gráfica:
>
> [image: Imágenes integradas 1]
>
>
>
> El 4 de diciembre de 2013 11:04, Carlos J. Gil Bellosta <
> cgb en datanalytics.com> escribió:
>
> Hola, ¿qué tal?
>>
>> Tu problema es que tienes agrupaciones de países (p.e., ASIA_MERIDIOL)
>> con un solo respresentante. Las matrices con las que construyes las
>> elipses son degeneradas (más bien, no existen: sospecho que tienen que
>> ver con la covarianza... ¿de un único caso?).
>>
>> Elimina las agrupaciones triviales.
>>
>> Un saludo,
>>
>> Carlos J. Gil Bellosta
>> http://www.datanalytics.com
>>
>> El día 4 de diciembre de 2013 16:58, Carlos Ortega
>> <cof en qualityexcellence.es> escribió:
>> > Hola,
>> >
>> > He mirado por encima el conjunto y el posible error.
>> > Como alternativa, mira la solución que se propone aquí:
>> >
>> >
>> http://stats.stackexchange.com/questions/24450/how-to-highlight-predefined-groups-in-pca-individual-map
>> >
>> > Y como sugerencias:
>> >
>> >    - Limpiaría un poco el conjunto para quitar las filas que tienen
>> todas
>> >    sus columnas con NA. ¿Qué sentido tiene dejarlas?
>> >    - Y haría el ejercicio de PCA con un conjunto más reducido de países
>> >    para ver cómo funciona.
>> >
>> > Otro consejo, para esta lista, es que digas qué paquetes estás usando
>> en tu
>> > ejemplo. PCA es una función de "FactoMineR" que no es un paquete por
>> > defecto de R.
>> >
>> > Saludos,
>> > Carlos Ortega
>> > www.qualityexcellence.es
>> >
>> >
>> > El 4 de diciembre de 2013 14:46, Camilo Calle <camicalle en gmail.com
>> >escribió:
>> >
>> >> Cordial saludo, soy principiante en R y estoy haciendo un análisis
>> >> multivariado de datos, y necesito hacer las elipses de confianza donde
>> la
>> >> variable 19 es cualitativa suplementaria y el resto son variables
>> activas;
>> >> tengo el siguiente código:
>> >>
>> >> eu60 <- read.csv(file.choose(), header=T, sep=";", dec=",",
>> row.names=1)
>> >> ### capturar datos
>> >> eu60.pca <- PCA(eu60, quali.sup=19) ## genero el objeto PCA
>> >> eu60.data <- cbind.data.frame(eu60[1:182,19], eu60.pca$ind$coord) ###
>> creo
>> >> el data.frame con las coordenadas de los individuos y la variable
>> pasiva
>> >> eu60.ellipse <- coord.ellipse(eu60.data, bary=TRUE) ## genero el objeto
>> >> ellise (acá me aparece el siguiente error):
>> >> Error en if (scale[1] > 0) r <- r/scale[1] :
>> >>   valor ausente donde TRUE/FALSE es necesario
>> >>
>> >> no puedo continuar con la siguiente rutina porque tengo un error en
>> este
>> >> paso, no sé por que se presenta si alguien sabe como puedo arreglar
>> esto le
>> >> agradecería enormemente.
>> >>
>> >> adjunto base de datos.
>> >>
>> >> --
>> >> Camilo Alberto Calle Velásquez
>> >> Zootecnista UdeA.
>> >>
>> >> _______________________________________________
>> >> R-help-es mailing list
>> >> R-help-es en r-project.org
>> >> https://stat.ethz.ch/mailman/listinfo/r-help-es
>> >>
>> >>
>> >
>> >
>> > --
>> > Saludos,
>> > Carlos Ortega
>> > www.qualityexcellence.es
>> >
>> >         [[alternative HTML version deleted]]
>> >
>> >
>> > _______________________________________________
>> > R-help-es mailing list
>> > R-help-es en r-project.org
>> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>> >
>>
>
>
>
> --
> Camilo Alberto Calle Velásquez
> Zootecnista UdeA.
>



-- 
Camilo Alberto Calle Velásquez
Zootecnista UdeA.
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20131204/bdfe0ccf/attachment-0001.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: no disponible
Type: image/png
Size: 88748 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20131204/bdfe0ccf/attachment-0001.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: eu601.csv
Type: text/csv
Size: 15410 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20131204/bdfe0ccf/attachment-0001.bin>