R help es - Apr 2012 - Resumen de R-help-es, Vol 38, Envío 13

2012/4/7 <r-help-es-request@r-project.org>
> Envíe los mensajes para la lista R-help-es a
>        r-help-es@r-project.org
>
> Para subscribirse o anular su subscripción a través de la WEB
>        https://stat.ethz.ch/mailman/listinfo/r-help-es
>
> O por correo electrónico, enviando un mensaje con el texto "help"
en
> el asunto (subject) o en el cuerpo a:
>        r-help-es-request@r-project.org
>
> Puede contactar con el responsable de la lista escribiendo a:
>        r-help-es-owner@r-project.org
>
> Si responde a algún contenido de este mensaje, por favor, edite la
> linea del asunto (subject) para que el texto sea mas especifico que:
> "Re: Contents of R-help-es digest...". Además, por favor, incluya
en
> la respuesta sólo aquellas partes del mensaje a las que está
> respondiendo.
>
> Asuntos del día:
>
>   1. grafica modelo logit y probit (sofia andrea ospina valencia)
>   2. Re: grafica modelo logit y probit (Carlos Ortega)
>
Hola Sofía:
                   no he seguido toda la línea de tu pregunta así que
disculpame si la información que te brinde no sea la adecuada o la
respuesta a tu pregunta la haya dado otro usuario.
                  Si te fijás en el libro: "Extending the linear model with
R" de Julian Faraway en la página 31, en el famoso ejemplo del challenger.
Utiliza la sentencia:
#paquete faraway
plot(damage/6 ~ temp, orings, xlim=c(25,85))
#damage/6 dado que da la proporción de fallos respecto de 6 pruebas
corresponde a una distribución binomial
x <- seq(25,85,1)
lines(x, ilogit(11.66-0.21*x))#la función ilogit del paquete faraway es la
inversa del logit y los valores de intercepto y pendiente corresponden al
ajuste del modelo binomial con función de enlace logit.
De esta manera  grafica el logit o digamos el logit estimado del modelo

si deseas graficar el probit en el mismo gráfico utiliza:
lines(x,pnorm(5.59-0.1058*X),lty=2)#el pnorm es la inversa del probit y los
valores de intercepto y pendiente son los valores estimados del modelo
probit.

los modelos ajustados de donde se obtuvieron los valores de ordenada y
pendiente anteriormente nombrados fueron:
logitmod <- glm (cbind(damage,6-damage)~temp,family=binomial,orings)
probitmod<-glm(cbind(damage,6-damage)~temp,family=binomial(link=probit),orings)

saludos!
>
>
> ---------- Mensaje reenviado ----------
> From: sofia andrea ospina valencia <zphya@yahoo.es>
> To: "r-help-es@r-project.org" <r-help-es@r-project.org>
> Cc:
> Date: Fri, 6 Apr 2012 20:00:06 +0100 (BST)
> Subject: [R-es] grafica modelo logit y probit
> Buenas noches
>
> Por favor me pueden indicar cual es la sentencia para realizar la grafica
> de un modelo logit y probit
>
> llevo horas buscando
>
> Gracias
>        [[alternative HTML version deleted]]
>
>
>
>
> ---------- Mensaje reenviado ----------
> From: Carlos Ortega <cof@qualityexcellence.es>
> To: sofia andrea ospina valencia <zphya@yahoo.es>
> Cc: "r-help-es@r-project.org" <r-help-es@r-project.org>
> Date: Sat, 7 Apr 2012 00:27:55 +0200
> Subject: Re: [R-es] grafica modelo logit y probit
> Hola,
>
> ¿ Has probado con los paquetes: "effects" y "glmulti" ?
>
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>
>
> El 6 de abril de 2012 21:00, sofia andrea ospina valencia
> <zphya@yahoo.es>escribió:
>
> > Buenas noches
> >
> > Por favor me pueden indicar cual es la sentencia para realizar la
grafica
> > de un modelo logit y probit
> >
> > llevo horas buscando
> >
> > 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]]
>
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>

-- 
Ing. Agr. Gustavo Giménez

	[[alternative HTML version deleted]]