data("USAArrests")

str(USAArrests)
'data.frame':   50 obs. of  4 variables:
 $ Murder  : num  13.2 10 8.1 8.8 9 7.9 3.3 5.9 15.4 17.4 ...
 $ Assault : int  236 263 294 190 276 204 110 238 335 211 ...
 $ UrbanPop: int  58 48 80 50 91 78 77 72 80 60 ...
 $ Rape    : num  21.2 44.5 31 19.5 40.6 38.7 11.1 15.8 31.9 25.8 ...


hist(USAArrests$murder)

train = USArrests[1:30,]

train$Murder = log(train$Murder)

test = USArrests[31:50,]

fit = lm(Murder~., data = train)

pred = predict(fit, test)

exp(pred)

library(bestNormalize)

train$Murder = yeojohnson(train$Murder)$x.t