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. ------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20131204/2f749c43/attachment.html> ------------ próxima parte ------------ A non-text attachment was scrubbed... Name: eu60.csv Type: text/csv Size: 16642 bytes Desc: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20131204/2f749c43/attachment.bin>
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@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@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]
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 >