<div dir="ltr">Cordial saludo adjunto la base de datos y el script:<div><br></div><div>eu60<- read.csv(file.choose(), header=T, sep=";", dec=".", row.names=1)<br></div><div>eu60.pca <- PCA(eu60, quali.sup=19)<br>
</div><div>eu60.data <- cbind.data.frame(eu60[,19], eu60.pca$ind$coord)<br></div><div>eu60.ellipse <- coord.ellipse(eu60.data, bary=TRUE)<br></div><div>plot.PCA(ellipse=eu60.ellipse, cex=0.8)</div><div><br></div></div>
<div class="gmail_extra"><br><br><div class="gmail_quote">El 4 de diciembre de 2013 17:35, Camilo Calle <span dir="ltr"><<a href="mailto:camicalle@gmail.com" target="_blank">camicalle@gmail.com</a>></span> escribió:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">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:<div>
<br></div><div><img src="cid:ii_142bfbfdca2d075a" alt="Imágenes integradas 1" width="691" height="368"><br>
<div><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">El 4 de diciembre de 2013 11:04, Carlos J. Gil Bellosta <span dir="ltr"><<a href="mailto:cgb@datanalytics.com" target="_blank">cgb@datanalytics.com</a>></span> escribió:<div>
<div class="h5"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">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ó:<br>
<div>> 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>
><br>
</div>>    - Limpiaría un poco el conjunto para quitar las filas que tienen todas<br>
<div>>    sus columnas con NA. ¿Qué sentido tiene dejarlas?<br>
</div>>    - Y haría el ejercicio de PCA con un conjunto más reducido de países<br>
<div><div>>    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>
><br>
</div></div>>         [[alternative HTML version deleted]]<br>
<div><div>><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>
</div></div></blockquote></div></div></div><div><div class="h5"><br><br clear="all"><div><br></div>-- <br>Camilo Alberto Calle Velásquez<br>Zootecnista UdeA.
</div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Camilo Alberto Calle Velásquez<br>Zootecnista UdeA.
</div>