data("CoalMiners")
## Fourfold display, both margins equated
fourfold(CoalMiners, mfcol = c(2,4))
## Log Odds Ratio Plot
summary(l <- oddsratio(CoalMiners))
g <- seq(25, 60, by = 5)
plot(l,
xlab = "Age Group",
main = "Breathlessness and Wheeze in Coal Miners")
m <- lm(l ~ g + I(g^2))
lines(fitted(m), col = "red")
## Fourfold display, strata equated
fourfold(CoalMiners, std = "ind.max", mfcol = c(2,4))
Run the code above in your browser using DataLab