(a) To investigate the distributions, simulate from this model 2500 times. To do so, calculate P[Y = 1 | X =x] for an observation, and then make a random draw from a Bernoulli distribution with that success probability. (Note that a Bernoulli distribution is a Binomial distribution with parameter n = 1. There is no direction function in R for a Bernoulli distribution.) Each time, fit the model: log D(x) 1 - p(x) = Bo + BIT1 + B212 + B313 Store the test statistics for two tests: . The Wald test for Ho : B2 =0, which we say follows a standard normal distribution for "large" samples