Homework 13, Part (a) Due Friday March 7, 2008. Part (b) Due Monday March 10,2008. ============================================== Generate a 10000-observation dataset HW13Data.txt using the command: > HW13Data = rlogis(10000, runif(1,2,2.5), runif(1,.7,.85)) This gives a pseudo-random sample of size 10000 from a logistic distribution with location parameter A and scale-parameter B, which are unknown parameters, and which you are to estimate by numerical maximization of the log-likelihood. (For the sake of good starting values, you may use the fact that the true parameter values fall in the interval (2,2.5) for A , and (.7,.85) for B. Your task in this exercise is to find the Maximum Likelihood Estimators and Standard Errors for A, B. Do this by setting up a direct numerical maximization using nlm. First do it without providing "gradient" and "hessian" attributes, and then do it again providing these attributes coded as efficiently as you can. Do "unix.time" timing runs both ways, to compare the computing time with versus without analytical gradient and hessian. (b) Also do the minimization by a direct Gradient search using the same log-likelihood function (using your same analytical gradient function), and do a "unix.time" timing run on it. NEEDLESS TO SAY, YOUR DIFFERENT MAXIMIZATION ROUTINES SHOULD ALL REACH THE SAME MAXIMUM !!