NiftyRscript1

creating some test data of invented body temperatures

bodytemp1 <- rnorm(n = 100, mean = 37.1) # create normally distributed random data
bodytemp2 <- rnorm(n = 100, mean = 36.5)

tempdata <- data.frame(c(bodytemp1, bodytemp2), c(rep("group1", 100), rep("group2", 100))) # create data frame

colnames(tempdata)[1] <- "bodytemp"
colnames(tempdata)[2] <- "group"

boxplots with individual data points using ggplot

require(ggplot2)

## Loading required package: ggplot2

print(qplot(group, bodytemp, data=tempdata, geom=c("boxplot")) + geom_point(alpha = 0.4) + theme_bw())