<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta name=Generator content="Microsoft Word 15 (filtered medium)"><!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EstiloCorreo17
{mso-style-type:personal;
font-family:"Calibri","sans-serif";
color:#1F497D;}
span.EstiloCorreo18
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 3.0cm 70.85pt 3.0cm;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=ES link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>Estimados<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>Al leer este interesante script lo corrí, pero al hacer la elipse me devuelve esa información, no encuentro el porque<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>saludos<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>library(FactoMineR)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>eu60<- read.csv(file.choose(), header=T, sep=";", dec=".", row.names=1)<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>eu60.pca <- PCA(eu60, quali.sup=19)<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>eu60.data <- cbind.data.frame(eu60[,19], eu60.pca$ind$coord)<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>eu60.ellipse <- coord.ellipse(eu60.data, bary=TRUE)<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>plot.PCA(ellipse=eu60.ellipse, cex=0.8)<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>#devuelve esta información<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>Error en plot.PCA(ellipse = eu60.ellipse, cex = 0.8) : <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'> el argumento "x" está ausente, sin valor por omisión<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><div><div style='border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'><o:p> </o:p></span></b></p><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'><o:p> </o:p></span></b></p><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'><o:p> </o:p></span></b></p><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'>De:</span></b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'> Dr. José A Betancourt Bethencourt [mailto:jbetancourt@iscmc.cmw.sld.cu] <br><b>Enviado el:</b> jueves, 5 de diciembre de 2013 04:45<br><b>Para:</b> 'Camilo Calle'<br><b>Asunto:</b> RE: [R-es] AYUDA CON ERROR CON LA LIBRERIA PCA<o:p></o:p></span></p></div></div><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'>Gracias!!<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'>De:</span></b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'> Camilo Calle [<a href="mailto:camicalle@gmail.com">mailto:camicalle@gmail.com</a>] <br><b>Enviado el:</b> miércoles, 4 de diciembre de 2013 06:28<br><b>Para:</b> Carlos J. Gil Bellosta; <a href="mailto:jbetancourt@iscmc.cmw.sld.cu">jbetancourt@iscmc.cmw.sld.cu</a><br><b>CC:</b> Lista R<br><b>Asunto:</b> Re: [R-es] AYUDA CON ERROR CON LA LIBRERIA PCA<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>Cordial saludo adjunto la base de datos y el script:<o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>eu60<- read.csv(file.choose(), header=T, sep=";", dec=".", row.names=1)<o:p></o:p></p></div><div><p class=MsoNormal>eu60.pca <- PCA(eu60, quali.sup=19)<o:p></o:p></p></div><div><p class=MsoNormal>eu60.data <- cbind.data.frame(eu60[,19], eu60.pca$ind$coord)<o:p></o:p></p></div><div><p class=MsoNormal>eu60.ellipse <- coord.ellipse(eu60.data, bary=TRUE)<o:p></o:p></p></div><div><p class=MsoNormal>plot.PCA(ellipse=eu60.ellipse, cex=0.8)<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div></div><div><p class=MsoNormal style='margin-bottom:12.0pt'><o:p> </o:p></p><div><p class=MsoNormal>El 4 de diciembre de 2013 17:35, Camilo Calle <<a href="mailto:camicalle@gmail.com" target="_blank">camicalle@gmail.com</a>> escribió:<o:p></o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt'><div><p class=MsoNormal>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:<o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><img border=0 width=691 height=368 id="_x0000_i1025" src="cid:image001.png@01CEF17A.7EEEB3C0" alt="Imágenes integradas 1"><o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div></div></div><div><p class=MsoNormal style='margin-bottom:12.0pt'><o:p> </o:p></p><div><p class=MsoNormal>El 4 de diciembre de 2013 11:04, Carlos J. Gil Bellosta <<a href="mailto:cgb@datanalytics.com" target="_blank">cgb@datanalytics.com</a>> escribió:<o:p></o:p></p><div><div><p class=MsoNormal><o:p> </o:p></p><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt'><p class=MsoNormal>Hola, ¿qué tal?<br><br>Tu problema es que tienes agrupaciones de países (p.e., ASIA_MERIDIOL)<br>con un solo respresentante. Las matrices con las que construyes las<br>elipses son degeneradas (más bien, no existen: sospecho que tienen que<br>ver con la covarianza... ¿de un único caso?).<br><br>Elimina las agrupaciones triviales.<br><br>Un saludo,<br><br>Carlos J. Gil Bellosta<br><a href="http://www.datanalytics.com" target="_blank">http://www.datanalytics.com</a><br><br>El día 4 de diciembre de 2013 16:58, Carlos Ortega<br><<a href="mailto:cof@qualityexcellence.es" target="_blank">cof@qualityexcellence.es</a>> escribió:<o:p></o:p></p><div><p class=MsoNormal>> Hola,<br>><br>> He mirado por encima el conjunto y el posible error.<br>> Como alternativa, mira la solución que se propone aquí:<br>><br>> <a href="http://stats.stackexchange.com/questions/24450/how-to-highlight-predefined-groups-in-pca-individual-map" target="_blank">http://stats.stackexchange.com/questions/24450/how-to-highlight-predefined-groups-in-pca-individual-map</a><br>><br>> Y como sugerencias:<br>><o:p></o:p></p></div><p class=MsoNormal>> - Limpiaría un poco el conjunto para quitar las filas que tienen todas<o:p></o:p></p><div><p class=MsoNormal>> sus columnas con NA. ¿Qué sentido tiene dejarlas?<o:p></o:p></p></div><p class=MsoNormal>> - Y haría el ejercicio de PCA con un conjunto más reducido de países<o:p></o:p></p><div><div><p class=MsoNormal>> para ver cómo funciona.<br>><br>> Otro consejo, para esta lista, es que digas qué paquetes estás usando en tu<br>> ejemplo. PCA es una función de "FactoMineR" que no es un paquete por<br>> defecto de R.<br>><br>> Saludos,<br>> Carlos Ortega<br>> <a href="http://www.qualityexcellence.es" target="_blank">www.qualityexcellence.es</a><br>><br>><br>> El 4 de diciembre de 2013 14:46, Camilo Calle <<a href="mailto:camicalle@gmail.com" target="_blank">camicalle@gmail.com</a>>escribió:<br>><br>>> Cordial saludo, soy principiante en R y estoy haciendo un análisis<br>>> multivariado de datos, y necesito hacer las elipses de confianza donde la<br>>> variable 19 es cualitativa suplementaria y el resto son variables activas;<br>>> tengo el siguiente código:<br>>><br>>> eu60 <- read.csv(file.choose(), header=T, sep=";", dec=",", row.names=1)<br>>> ### capturar datos<br>>> eu60.pca <- PCA(eu60, quali.sup=19) ## genero el objeto PCA<br>>> eu60.data <- cbind.data.frame(eu60[1:182,19], eu60.pca$ind$coord) ### creo<br>>> el data.frame con las coordenadas de los individuos y la variable pasiva<br>>> eu60.ellipse <- coord.ellipse(eu60.data, bary=TRUE) ## genero el objeto<br>>> ellise (acá me aparece el siguiente error):<br>>> Error en if (scale[1] > 0) r <- r/scale[1] :<br>>> valor ausente donde TRUE/FALSE es necesario<br>>><br>>> no puedo continuar con la siguiente rutina porque tengo un error en este<br>>> paso, no sé por que se presenta si alguien sabe como puedo arreglar esto le<br>>> agradecería enormemente.<br>>><br>>> adjunto base de datos.<br>>><br>>> --<br>>> Camilo Alberto Calle Velásquez<br>>> Zootecnista UdeA.<br>>><br>>> _______________________________________________<br>>> R-help-es mailing list<br>>> <a href="mailto:R-help-es@r-project.org" target="_blank">R-help-es@r-project.org</a><br>>> <a href="https://stat.ethz.ch/mailman/listinfo/r-help-es" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-help-es</a><br>>><br>>><br>><br>><br>> --<br>> Saludos,<br>> Carlos Ortega<br>> <a href="http://www.qualityexcellence.es" target="_blank">www.qualityexcellence.es</a><br>><o:p></o:p></p></div></div><p class=MsoNormal>> [[alternative HTML version deleted]]<o:p></o:p></p><div><div><p class=MsoNormal>><br>><br>> _______________________________________________<br>> R-help-es mailing list<br>> <a href="mailto:R-help-es@r-project.org" target="_blank">R-help-es@r-project.org</a><br>> <a href="https://stat.ethz.ch/mailman/listinfo/r-help-es" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-help-es</a><br>><o:p></o:p></p></div></div></blockquote></div></div></div><div><div><p class=MsoNormal><br><br clear=all><o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><p class=MsoNormal>-- <br>Camilo Alberto Calle Velásquez<br>Zootecnista UdeA. <o:p></o:p></p></div></div></div></blockquote></div><p class=MsoNormal><br><br clear=all><o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><p class=MsoNormal>-- <br>Camilo Alberto Calle Velásquez<br>Zootecnista UdeA. <o:p></o:p></p></div></div></body></html>