STAT7030 ANU Generalised Linear Models

Question 1. [48 Marks]. The file Doctor Data, available on the course Wattle page in theDatasets section, is concerned with the number of deaths by heart disease in a sample of
doctors from different age groups and with different smoking status (smoker/non-smoker).
Further details are given below.
Column 1 Observation number – [Variable V1=1 to 10]
Column 2 Age – [Variable V2= 1,2,3,4,5] {1 = 35-44 years; 2 = 45-54 years;
3 = 55-64 years; 4 = 65-74 years; 5 = 75-84 years}.
Column 3 Deaths – [Varaible V3 = Number of Deaths] {Count Variable}.
Column 4 Person-Years at Risk – [Variable V4 = Aggregate years] {Count Variable}.
Column 5 Smoker – [Variable V5 = 1,2] { Factor: 1 = Smoker; 2 = Non-Smoker}.
Further comment: do not forget to declare variable V 5 as a factor.
Please provide answers to the following.
(a) Fit a suitable set of Poission regression models to the data and select what you think is
the best model, using AIC to inform your selection. Please bear the following points
in mind:
• Treat log(V 4) as a covariate (as opposed to an offset) throughout part (a).
• Although variable V 2 could be treated as a factor, it is suggested that you treat
it as a numerical variable. If treated as a factor, Age will use up too many
parameters relative to the sample size. Also, consider the covariate V 2 ∗ V 2.
• Investigate interactions as well as main effects.
• Present the AIC values in a suitable table.
• Use suitable residual plots and other diagnostics to investigate the goodness of
fit of the model you select.
• In part (a) and in part (b) below, take the scale parameter φ to be 1.
[12 marks]
(b) Focusing on the model you selected in part (a), perform a change of deviance test to
assess whether it would be better to treat log(V 4) as an offset or as a covariate.
[6 marks]
(c) Repeat parts (a) and (b) using negative binomial regression models with a log link.
Comment on (i) similarities and (ii) differences in your results for part (c), compared
with parts (a) and (b).
[18 marks]
2
(d) Returning to the best Poisson regression model that you identified in part (a), estimate
the scale parameter φ and provide a Wald test of the null hypothesis that φ = 1,
presenting your answer in the form of a p-value. As part of your answer to part (d)
you should present the R code that you used.
[12 marks]
3
1 1 32 52407 1
2 2 104 43248 1
3 3 206 28612 1
4 4 186 12663 1
5 5 102 5317 1
6 1 2 18790 2
7 2 12 10673 2
8 3 28 5710 2
9 4 28 2585 2
10 5 31 1462 2