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 >
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. ------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20131204/5f70afa6/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/5f70afa6/attachment-0001.png>
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>