data(aacost)
levels(aacost$gc) <- c("low G+C", "mid G+C", "high G+C")
stripchart(aacost$tot~aacost$gc, pch = 19, ylim = c(0.5,3.5),
xlim = c(0, max(aacost$tot)),
xlab = "Metabolic cost (high-energy phosphate bonds equivalent)",
main = "Metabolic cost of the 20 amino-acids
as function of their G+C class" )
boxplot(aacost$tot~aacost$gc, horizontal = TRUE, add = TRUE)
Run the code above in your browser using DataLab