Description
1. Cancer of Tongue.
Sickle-Santanello et al (1988)1 provide data on 80 males diagnosed
with cancer of the tongue. Data are provided in the file tongue.csv|dat|xlsx.
The variables
in the dataset are as follows:
• Tumor DNA profile (1 – aneuploid tumor, 2 – diploid tumor);
• Time to death or on-study time (in weeks); and
• Censoring indicator (0=censored, 1=observed)
Fit the regression with tumor profile as covariate. What is the 95% Credible Set for the
slope β1?
2. Airfreight Breakage with Missing Data.
A substance used in biological and medical
research is shipped by air freight to users in cartons of 2,000 ampules. The data below,
involving 15 shipments, were collected on the number of times a carton was transferred from
one aircraft to another over the shipment route (X) and the number of ampules found to be
broken upon arrival (Y ).
X 2 1 0 2 NA 3 1 0 1 2 3 0 1 NA NA
Y NA 16 9 17 12 22 13 8 NA 19 17 11 10 20 2
(a) Using OpenBUGS/WinBUGS, fit Y by Poisson regression, with X as a covariate.
Report the deviance of your fit.
(b) According to your model, how many packages on average are expected will be broken
if the number of shipment routes is X = 4? What is 95% CS for your estimate.
(c) For a particular shipment sent from Shenzhen you learned that it would involve X = 4
shipping routes. Predict the number of broken packages. What is here different from (b)?
(d) What are estimates for unobserved X5, X14, X15, Y1 and Y9?
Hint: Note that for missing X, you need to specify the distribution. It could be any nonnegative valued distribution, but since X is discrete, a good choice is Poisson(2), as in
for(i in 1:n){
x[i] ~ dpois(2)}
1Sickle-Santanello, B. J., Farrar, W. B., DeCenzo, J. F., Keyhani-Rofagha, S., Klein, J., Pearl, D.,
Laufman, H., and O’Toole R. V. (1988). Technical and statistical improvements for flow cytometric DNA
analysis of paraffin-embedded tissue. Cytometry, 9, 594–599.
1