Sample code from *R CookBook* (awesome book btw) *11.12: Finding the Best
Power Transformation (Box-Cox) Procedure*

require(MASS)
x <- 10:100
eps <- rnorm(length(x), sd = 5)
y <- (x + eps)^(-1 / 1.5)

m <- lm(y ~ x)

# ***************** What does the *which* in this line do???
********************************************
plot(m, which = 1)

	[[alternative HTML version deleted]]

