glasses <- read.table(hh("datasets/glasses.dat"), header=FALSE)
glasses <- matrix(unlist(glasses),
nrow=2,
dimnames=list(
c("glasses","no.glasses"),
c("delinq","non.del")))
## draw the iso-odds ratio plot with 50\% CI and 95\% CI,
## invisibly return the 95\% CI
plot.odds.ratio(glasses)
Run the code above in your browser using DataLab