[R] probit plots

[Anne_Yeh@hc-sc.gc.ca]

Hello,

I'm trying to do a probit plot but am a little confused about how to do this.

Sample data:

 >  dat <-
data.frame(dose=c(0.18,0.48,0.78,0.18,0.48,0.78),n=c(103,120,123,60,110,100),r=c(19,53,83,14,54,81),

+  sample=c("s","s","s","t","t","t"))

> temp <- glm(cbind(r,n-r)~sample+log(dose),binomial(link='probit'),data=dat)

Does anyone know how to plot the probit of the response vs. log(dose) with the 2
 regression lines for the two treatment groups
superimposed onto one graph?

Any comments will be much appreciated,

Anne


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._