From: Anne_Yeh@hc-sc.gc.ca
Date: Fri 20 Jul 2001 - 03:16:52 EST
Message-id: <85256A8E.005EF24F.00@smta00.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@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
This archive was generated by hypermail 2.1.3 : Wed 10 Mar 2004 - 08:18:20 EST