source("http://192.38.117.59/~linearpredictors/datafiles/readVitaminD.R")

## Factor with levels corresponding to normal resp. overweight
irlwomen$bmigroup2 <- (irlwomen$bmi > 0) + (irlwomen$bmi > 25)
irlwomen$bmigroup2 <- as.factor(irlwomen$bmigroup2)

## Scatterplot of Vitamin D against the two factor levels
offset <- c(5/4, 7/4)
stripchart(log10(irlwomen$vitd) ~ irlwomen$bmigroup2,
           vertical = TRUE,
           pch = 21,
           method = "jitter", jitter = 0.01,
           ylab = expression(paste(log[10],"(vitamin D)")),
           group.names = c("BMI < 25", "25 < BMI"),
           at = offset ## where to place groups at x axis
           )

## Marking groupwise means of log10 of Vitamin D
xbar <- tapply(log10(irlwomen$vitd), irlwomen$bmigroup2, mean)
r <- 0.03
segments(offset-r, xbar, offset+r, xbar)