[Solved] MAST30025 Lab 6

Questions 17 use the sleep dataset, which you can download from the course website. This dataset contains (among other things) data on the body weight (kg) and brain weight (g) of 62 mammals. Use the following commands to read the data:

> mammals <- read.csv(../data/sleep.csv)

> mammals$BodyWt <- log(mammals$BodyWt)

> mammals$BrainWt <- log(mammals$BrainWt)

This creates a data frame, mammals, with components (among others) named BodyWt and BrainWt, then applies a logarithmic transformation to both BodyWt and BrainWt.

1. Fit a linear model explaining brain weight from body weight, using the lm

Display the summary of the fitted model, and then create a scatter plot of the data and superimpose the fitted regression line on it. Does it look like a reasonable fit?

Use diagnostic plots to assess if the model assumptions are satisfied.

2. Using the fitted model or otherwise, obtain:
   • The least squares estimator of the parameters, b;
   • The vector of residuals, e;
   • The residual sum of squares, SS_Res;
   • The regression sum of squares, SS_Reg;
   • The estimator for the variance of the errors, s²;
   • The standardised residuals;
   • The leverages of the points;
   • The Cooks distances of the points;
   • 95% confidence intervals for each of the parameters.
3. Find a 95% confidence interval for a mammal weighing 50 kg.
4. Find a 95% prediction interval for a mammal weighing 50 kg.
5. Find and draw a 95% joint confidence region for the parameters.
6. Test the following hypotheses, using the anova
   • H₀ : = 0
   • H₀ : ₁ = 0
   • H₀ : ₀ = 0
   • H₀ : = (2,1)
7. By visualising the raw data, justify the use of a double logarithmic transformation. Write downthe final model for the (untransformed) brain weight vs. body weight.
8. In this question we consider the hypothesis H₀ : = . The test statistic for this hypothesis is

.

1

• Show that

(b)^TX^TX(b) = (y X)^T(y X) (y Xb)^T(y Xb).

That is, it is the SS_Res for the null model minus the SS_Res for the full model. Also show that

(b)^TX^TX(b) 6= y^TX(X^TX)¹X^TyX^TX.

That is, in this case we can not write it as the SS_Reg for the full model minus the SS_Reg for the model under H₀.

• Show directly that (b)^TX^TX(b) and SS_Res are independent, that is without using our existing results that b and SS_Res are independent.

Hint: set q = y X then

1. Show that (b)^TX^TX(b) = q^TX(X^TX)¹X^Tq.
2. Show that SS_Res = q^T[I X(X^TX)¹X^T]q and hence that these two quadratic forms are independent.

   

   6