Writing a new script to test the variational approximation on
idealized simulated data with a binary trait.