[R] Inverse of Probit

Calum calummclean at hotmail.com
Tue Nov 19 19:59:01 CET 2013


Hi there,
I hope someone can help me.

I have a dataset of Concentration against Mortality, and I am trying to
compare the use of Logit and Probit models using this data.

The issue I am having is trying to back transform the data from the probit
model, to plot it in normal space instead of log space.
I know this can be done with a logit model using the code below, where
ilogit is a function for the inverse logit:

NEWCONC <- seq(0,0.6, length=25)
NEWMORT <- predict(LOGIT, Conc=NEWCONC, se=TRUE)

plot(data=DATA, Prop~Conc)
lines(NEWCONC, ilogit(NEWMORT$fit))

However, I can't seem to find a function equivalent to ilogit for a probit
model, that I could use in this code:

NEWCONC <- seq(0,0.6, length=25)
NEWMORT <- predict(PROBIT, Conc=NEWCONC, se=TRUE)

plot(data=DATA, Prop~Conc)
lines(NEWCONC,###INVERSE PROBIT### (NEWMORT$fit))


Any advice on this issue would be appreciated,
Thanks,
Calum



--
View this message in context: http://r.789695.n4.nabble.com/Inverse-of-Probit-tp4680752.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list