(a) Use the function optim()to compute β^ using initial value (.25,.75). (b) Again, starting with (.25,.75) find the next value when using the Newton-Raphson algorithm. (c) Assume that β0=0, and plot the likelihood function L(β1) as a function of β1. (d) Again, assume β0=0 and compute β^1 using uniroot(), a grid search, and by the Newton-Raphson algorithm. You can use the plot in part (c) to find a good initial value.