y = .5 + 1.2*x2 + 5*x10 + rnormal()