library(foreign)
daisy2<-read.dta("r:/daisy2red.dta")
library(MASS)
bc <- boxcox(milk120 ~ parity, data=daisy2, lambda=seq(0, 1, 0.01))
bc