move the subset statement to the munge code
add some simulation for logistic regression with quadratic interaction term