sofia andrea ospina valencia
2012-Mar-19 14:42 UTC
[R-es] elipse de confianza mod sin intercepto
Buenas tardes Estoy intentando representar graficamente la elipse de confianza de dos pendiente en un modelo sin intercepto La formula que introduzo esplot(ellipse(recta,c(1,1)),type="l"), pero elipse me sale una recta Gracias [[alternative HTML version deleted]]
Hola. Adjunta código reproducible de ejemplo y tal vez alguien pueda ayudarte. Un saludo Gregorio R. Serrano El 19 de marzo de 2012 15:42, sofia andrea ospina valencia <zphya@yahoo.es>escribió:> Buenas tardes > > Estoy intentando representar graficamente la elipse de confianza de dos > pendiente en un modelo sin intercepto > La formula que introduzo esplot(ellipse(recta,c(1,1)),type="l"), pero > elipse me sale una recta > > Gracias > [[alternative HTML version deleted]] > > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >-- Dr. Gregorio R. Serrano Dpto. Economía Cuantitativa (UCM) Voz:+34 91394 2361 Twitter: @grserrano_ http://www.grserrano.es [[alternative HTML version deleted]]
Hola Sofia, Revisa el primer ejemplo de require(car) ?ellipse Creo que lo unico que necesitas es tu modelo y listo. Siguiendo con el ejemplo, podrias en tu caso considerar algo similar a: confidenceEllipse(lm(prestige~income+education, data=Prestige), Scheffe=TRUE) Un saludo, Jorge.- PD: Si esto no es en definitiva lo que buscas, podrias por favor incluir en un proximo mensaje lo que has hecho (e.g., codigo en R), asi como los paquetes que has utilizado? Toma algo de tiempo encontrar las funciones :( 2012/3/19 sofia andrea ospina valencia <>> Buenas tardes > > Estoy intentando representar graficamente la elipse de confianza de dos > pendiente en un modelo sin intercepto > La formula que introduzo esplot(ellipse(recta,c(1,1)),type="l"), pero > elipse me sale una recta > > Gracias > [[alternative HTML version deleted]] > > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >[[alternative HTML version deleted]]
Hola, Además de lo que ya te han comentado, añadir que como tal "esplot()" no es una función soportada en ningún paquete de R. La única referencia que he podido encontrar [1] es el código que en su momento creó Deepayan (autor de lattice) y por ejemplo que incluye, veo que la capacidad de dibujar las distribuciones, así como un "rug" en el eje X existe ya disponible en otros paquetes/librerías como las que te han indicado. [1]: http://pages.stat.wisc.edu/~deepayan/771/esplot.R Saludos, Carlos Ortega www.qualityexcellence.es El 19 de marzo de 2012 15:42, sofia andrea ospina valencia <zphya@yahoo.es>escribió:> Buenas tardes > > Estoy intentando representar graficamente la elipse de confianza de dos > pendiente en un modelo sin intercepto > La formula que introduzo esplot(ellipse(recta,c(1,1)),type="l"), pero > elipse me sale una recta > > Gracias > [[alternative HTML version deleted]] > > > _______________________________________________ > 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]]
Sofia, Si el modelo en el que estas interesada es _sin_ intercepto, la sintaxis correcta (para el modelo en la parte inferior de este mensaje) es confidenceEllipse(lm(prestige~ income + education - 1, data = Prestige), Scheffe=TRUE) Espero no haber generado (aun mas) confusion. Jorge.- 2012/3/19 Jorge I Velez <>> Hola Sofia, > > Revisa el primer ejemplo de > > require(car) > ?ellipse > > Creo que lo unico que necesitas es tu modelo y listo. Siguiendo con el > ejemplo, podrias en tu caso considerar algo similar a: > > confidenceEllipse(lm(prestige~income+education, data=Prestige), > Scheffe=TRUE) > > Un saludo, > Jorge.- > > PD: Si esto no es en definitiva lo que buscas, podrias por favor incluir > en un proximo mensaje lo que has hecho (e.g., codigo en R), asi como los > paquetes que has utilizado? Toma algo de tiempo encontrar las funciones :( > > > 2012/3/19 sofia andrea ospina valencia <> > >> Buenas tardes >> >> Estoy intentando representar graficamente la elipse de confianza de dos >> pendiente en un modelo sin intercepto >> La formula que introduzo esplot(ellipse(recta,c(1,1)),type="l"), pero >> elipse me sale una recta >> >> Gracias >> [[alternative HTML version deleted]] >> >> >> _______________________________________________ >> R-help-es mailing list >> R-help-es@r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-help-es >> >> >[[alternative HTML version deleted]]