#your code to regress y on x and find b1 bivec [i] = b1 #your code to construct a 95%, confidence interval for betal #your code to compute mse msevec [i] = mse #the confidence interval is named as ci_b1 if (ci_b1 [1] = beta1) { coverage = coverage+1 overage . rate = coveragesim (a) The total number of simulations is denoted by nsim. Use nsim = 20, 500, 1000, and check the coverage rates. Are they close to 95%? (b) For each nsim value, make a histogram of bj's. What is the mean of bj's? Is the distribution symmetric around the mean? (c) Compute Sex, and compute o by the mean of values kept in msevec, i.e., E(MSE). Then compute the estimated standard deviation of bj by its math form for nsim = 20, 500, 1000. For the three simulation scenarios: nsim = 20, 500, 1000, rank the value of the estimated standard deviation of b1 from highest to lowest. (d) Also, rank the absolute value of the mean of bi's - B1, i.e., [E(b1) - B1| from highest to lowest for the three simulation scenarios