# NOT RUN {
x <- model.matrix(~ factor(wool) + factor(tension), warpbreaks)
y <- warpbreaks$breaks
warm_up <- small.em(x = x, y = y, K = 2, b.init = "random", sample.size = 20)
select_best(warm_up)
# }
Run the code above in your browser using DataLab