data(CoalMiners)
## Fourfold display, both margins equated
fourfoldplot(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
fourfoldplot(CoalMiners, std = "ind.max", mfcol = c(2,4))
Run the code above in your browser using DataLab