Writing a new script to test the variational approximation on
idealized simulated data with a binary trait (second commit attempt).