classes <- c("1: None","2: 0-5","3: 6-14","4: 15-24","5: 25-49","6: 50+")
classes

lung.cancer.patients <- c(7,55,489,475,293,38)
control.patients <- c(61,129,570,431,154,12)

data <- cbind(lung.cancer.patients,control.patients)
row.names(data) <- classes

row.totals <- rowSums(data)
col.totals <- colSums(data)

emp.probs <- lung.cancer.patients/row.totals

control.patients/row.totals
odds.ratio.1 <- (55/184)/(7/68)
odds.ratio.5 <- (38/50)/(7/68)

odds.ratio.vector <- c(odds.ratio.1,odds.ratio.5)
comparison <- c("1: 0-5 v. none","5: 50+ v. none")
library(lattice)
dotplot(odds.ratio.vector ~ comparison,ylab="Odds ratio")