rep(0 statt colMeans for the mean of mvrnorm noise