sex<-factor(rep(c("M","F"),times=c(5,5))) mean.age<-c(23.2,36.5,54.3,69.2,79.5,23.2,36.5,54.3,69.2,79.5) taking<-c(21,32,70,43,19,46,89,169,98,51) total<-c(704,628,775,338,118,784,789,1016,434,247) drugs<-data.frame(sex=sex,mean.age=mean.age,taking=taking,total=total) response<-taking/total logodds<-log(response/(1-response))