powered by
genCorrelatedData(N = 100, means = c(50, 50), sds = c(10, 10), rho = 0, stde = 1, beta = c(0, 0.2, 0.2, 0))