Using R and please show the codes.

Question 2: Use the commands data ("beavers"); y-beaver1$temp to load a dataset on the body temperature of a beaver from Wisconsin. Fit a Generalized Gamma model to the data. Plot a density estimator of the data and overlay a predictive density plot from the fitted model. This model clearly does not fit the data correctly and so we probably should not use it to compute a confidence interval for the mean. Generate 2000 bootstrap samples and compute their sample means Compute the naive 95% bootstrap confidence interval L, U, which directly uses the quantiles of the bootstrap sample means. Compute also the pivot interval (Lp, Up where Lp-2y -Un and Up - 2y - Ln. Comment on these two intervals. Is there one that you would prefer to the other? Question 2: Use the commands data ("beavers"); y-beaver1$temp to load a dataset on the body temperature of a beaver from Wisconsin. Fit a Generalized Gamma model to the data. Plot a density estimator of the data and overlay a predictive density plot from the fitted model. This model clearly does not fit the data correctly and so we probably should not use it to compute a confidence interval for the mean. Generate 2000 bootstrap samples and compute their sample means Compute the naive 95% bootstrap confidence interval L, U, which directly uses the quantiles of the bootstrap sample means. Compute also the pivot interval (Lp, Up where Lp-2y -Un and Up - 2y - Ln. Comment on these two intervals. Is there one that you would prefer to the other