Here, Y = {y1, y2, ..., y10} = {10, 9, 21, 4, 14, 11, 22, 9, 6, 4}. Using the (log-)likelihood function, how I can demonstrate the maximum likelihood estimation in STATA ?