library(foreign)
daisy2<-read.dta('h:/vhm/vhm802/data_stata/daisy2red.dta')
library(MASS)
bc <- boxcox(milk120 ~ parity, data=daisy2, lambda=seq(0, 1, 0.01))
bc