powered by
create_syndata_mv(Z, W, max_var_arr = rep(1, length(W)))
Each dataset is a data.frame to gain compatibility with lm() and glm()