<div dir="ltr">había un error en la ultima línea de código:<div><br><div>plot.PCA(eu60.pca, habillage=19, ellipse = eu60.ellipse, cex=0.8)<br></div></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
El 5 de diciembre de 2013 05:31, Dr. José A Betancourt Bethencourt <span dir="ltr"><<a href="mailto:jbetancourt@iscmc.cmw.sld.cu" target="_blank">jbetancourt@iscmc.cmw.sld.cu</a>></span> escribió:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div lang="ES" link="blue" vlink="purple"><div><p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Estimados<u></u><u></u></span></p><p class="MsoNormal">
<span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Al leer este interesante script lo corrí, pero al hacer la elipse me devuelve esa información, no encuentro el porque<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">saludos<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">library(FactoMineR)<u></u><u></u></span></p><div class="im"><p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">eu60<- read.csv(file.choose(), header=T, sep=";", dec=".", row.names=1)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">eu60.pca <- PCA(eu60, quali.sup=19)<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">eu60.data <- cbind.data.frame(eu60[,19], eu60.pca$ind$coord)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">eu60.ellipse <- coord.ellipse(eu60.data, bary=TRUE)<u></u><u></u></span></p><p class="MsoNormal">
<span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">plot.PCA(ellipse=eu60.ellipse, cex=0.8)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p></div><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">#devuelve esta información<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Error en plot.PCA(ellipse = eu60.ellipse, cex = 0.8) : <u></u><u></u></span></p><p class="MsoNormal">
<span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">  el argumento "x" está ausente, sin valor por omisión<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></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""><u></u> <u></u></span></b></p>
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri","sans-serif""><u></u> <u></u></span></b></p><p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri","sans-serif""><u></u> <u></u></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:<a href="mailto:jbetancourt@iscmc.cmw.sld.cu" target="_blank">jbetancourt@iscmc.cmw.sld.cu</a>] <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<u></u><u></u></span></p></div></div><div class="im"><p class="MsoNormal">
<u></u> <u></u></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Gracias!!<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></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" target="_blank">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" target="_blank">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<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Cordial saludo adjunto la base de datos y el script:<u></u><u></u></p><div><div class="h5"><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">
eu60<- read.csv(file.choose(), header=T, sep=";", dec=".", row.names=1)<u></u><u></u></p></div><div><p class="MsoNormal">eu60.pca <- PCA(eu60, quali.sup=19)<u></u><u></u></p></div><div><p class="MsoNormal">
eu60.data <- cbind.data.frame(eu60[,19], eu60.pca$ind$coord)<u></u><u></u></p></div><div><p class="MsoNormal">eu60.ellipse <- coord.ellipse(eu60.data, bary=TRUE)<u></u><u></u></p></div><div><p class="MsoNormal">plot.PCA(ellipse=eu60.ellipse, cex=0.8)<u></u><u></u></p>
</div><div><p class="MsoNormal"><u></u> <u></u></p></div></div></div></div><div><div class="h5"><div><p class="MsoNormal" style="margin-bottom:12.0pt"><u></u> <u></u></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ó:<u></u><u></u></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:<u></u><u></u></p>
<div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><img border="0" width="691" height="368" src="cid:image001.png@01CEF17A.7EEEB3C0" alt="Imágenes integradas 1"><u></u><u></u></p><div><p class="MsoNormal">
<u></u> <u></u></p></div></div></div><div><p class="MsoNormal" style="margin-bottom:12.0pt"><u></u> <u></u></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ó:<u></u><u></u></p>
<div><div><p class="MsoNormal"><u></u> <u></u></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ó:<u></u><u></u></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>><u></u><u></u></p></div><p class="MsoNormal">>    - Limpiaría un poco el conjunto para quitar las filas que tienen todas<u></u><u></u></p><div><p class="MsoNormal">>    sus columnas con NA. ¿Qué sentido tiene dejarlas?<u></u><u></u></p>
</div><p class="MsoNormal">>    - Y haría el ejercicio de PCA con un conjunto más reducido de países<u></u><u></u></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>><u></u><u></u></p></div></div><p class="MsoNormal">
>         [[alternative HTML version deleted]]<u></u><u></u></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>><u></u><u></u></p></div></div></blockquote></div></div></div><div><div><p class="MsoNormal">
<br><br clear="all"><u></u><u></u></p><div><p class="MsoNormal"><u></u> <u></u></p></div><p class="MsoNormal">-- <br>Camilo Alberto Calle Velásquez<br>Zootecnista UdeA. <u></u><u></u></p></div></div></div></blockquote></div>
<p class="MsoNormal"><br><br clear="all"><u></u><u></u></p><div><p class="MsoNormal"><u></u> <u></u></p></div><p class="MsoNormal">-- <br>Camilo Alberto Calle Velásquez<br>Zootecnista UdeA. <u></u><u></u></p></div></div></div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Camilo Alberto Calle Velásquez<br>Zootecnista UdeA.
</div>