The following data set (in WinBUGS format) refers to a smoking cessation experiment. Several smokers (actual number not of interest) were entered into the study and 261 of them quit. These 261 quitters were monitored and checked when they relapsed. RelapseT is the variable giving the time to relapse for an individual and is the dependent variable (like survival time). Age (age of the individual), AgeStart (age when he/she started smoking), SexF (Female=1, Male=0), SIUC (whether the individual received an intervention or not) and F10Cigs (the average number of cigarettes smoked per day) are independent covariates. In addition, censored.time is the vector of censorship times as is specified by the WinBUGS syntax. Write a WinBUGS code to carry out survival analysis, where survival is "no relapse" and failure is "relapse into smoking" . You may want to look closely at the mice example given in Vol I of the Examples, which has also been discussed in class. Write a brief report saying how many iterations you ran and present your summary results for the parameter estimates. model { #Your code here } #Your Inits here list() #datainput supplied list(Age = c(51, 44, 49, 47, 50, 55, 59, 41, 55, 42, 51, 49, 49, 50, 42, 54, 44, 44, 36, 43, 51, 54, 47, 47, 51, 49, 58, 42, 41, 36, 53, 35, 44, 54, 37, 49, 45, 45, 44, 46, 41, 55, 46, 40, 58, 37, 47, 53, 55, 36, 54, 43, 51, 57, 58, 50, 50, 36, 49, 36, 55, 46, 41, 59, 53, 46, 57, 44, 48, 44, 59, 42, 54, 42, 58, 40, 35, 43, 45, 49, 47, 59, 39, 41, 55, 45, 44, 58, 52, 42, 41, 54, 50, 51, 50, 36, 53, 57, 36, 43, 54, 39, 40, 35, 50, 56, 40, 48, 57, 52, 53, 59, 47, 39, 49, 45, 43, 53, 59, 42, 52, 47, 58, 44, 57, 42, 35, 44, 56, 43, 45, 37, 41, 36, 50, 42, 58, 45, 47, 58, 53, 58, 58, 46, 57, 45, 50, 50, 42, 51, 53, 41, 45, 52, 52, 54, 41, 48, 44, 55, 57, 48, 58, 45, 58, 59, 49, 49, 50, 46, 57, 44, 44, 39, 58, 56, 58, 56, 58, 42, 50, 52, 48, 46, 55, 51, 46, 42, 46, 51, 40, 50, 45, 51, 50, 41, 35, 43, 59, 54, 38, 52, 54, 49, 55, 41, 53, 51, 47, 40, 50, 58, 47, 55, 59, 48, 49, 41, 49, 54, 54, 52, 38, 45, 48, 51, 46, 49, 54, 53, 53, 48, 50, 52, 51, 58, 42, 50, 50, 48, 55, 59, 58, 46, 47, 48, 56, 52, 48, 51, 59, 51, 47, 49, 48, 42, 56, 48, 49, 40, 37, 57, 58), SexF = c(0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 1, 0), AgeStart = c(13, 16, 18, 14, 19, 15, 18, 16, 15, 18, 18, 18, 18, 17, 17, 22, 18, 18, 18, 20, 16, 20, 13, 16, 18, 16, 18, 17, 17, 19, 18, 19, 10, 21, 25, 17, 24, 15, 18, 25, 23, 21, 15, 17, 16, 17, 20, 12, 14, 17, 21, 17, 15, 15, 25, 18, 17, 21, 20, 16, 17, 18, 17, 13, 20, 17, 20, 16, 17, 15, 17, 19, 14, 16, 18, 20, 17, 26, 17, 19, 20, 26, 20, 19, 17, 18, 19, 19, 20, 21, 17, 21, 32, 20, 13, 18, 19, 21, 18, 18, 20, 20, 15, 20, 15, 20, 17, 15, 20, 15, 19, 17, 18, 16, 16, 18, 20, 21, 18, 18, 19, 17, 18, 17, 20, 23, 18, 17, 20, 18, 20, 19, 15, 20, 18, 16, 30, 19, 21, 28, 17, 13, 18, 20, 20, 20, 15, 18, 18, 22, 21, 18, 12, 20, 17, 35, 15, 14, 18, 16, 18, 23, 18, 13, 20, 20, 20, 18, 18, 14, 24, 19, 15, 19, 17, 19, 14, 20, 16, 15, 16, 16, 18, 18, 18, 18, 18, 17, 18, 18, 18, 17, 21, 18, 38, 17, 17, 15, 25, 18, 20, 18, 22, 21, 18, 18, 17, 15, 19, 18, 19, 20, 19, 21, 21, 18, 24, 16, 17, 16, 18, 17, 18, 15, 19, 18, 17, 19, 17, 18, 17, 17, 20, 20, 12, 16, 16, 15, 7, 13, 18, 20, 20, 14, 21, 18, 18, 23, 19, 24, 22, 14, 20, 12, 17, 19, 10, 15, 19, 16, 18, 17, 21), SIUC = c(1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0), F10Cigs = c(40, 20, 30, 20, 12, 40, 20, 40, 18, 40, 20, 18, 30, 30, 20, 60, 30, 30, 30, 25, 30, 30, 50, 40, 15, 50, 30, 23, 30, 25, 17, 18, 35, 25, 20, 30, 20, 40, 20, 5, 5, 40, 40, 16, 20, 20, 30, 25, 20, 25, 10, 10, 25, 20, 20, 30, 27, 20, 40, 35, 40, 25, 30, 40, 25, 50, 20, 20, 30, 15, 15, 10, 35, 30, 30, 25, 20, 25, 25, 40, 15, 20, 15, 20, 30, 20, 40, 15, 15, 20, 40, 30, 25, 40, 30, 35, 40, 20, 20, 25, 30, 20, 20, 40, 20, 35, 25, 20, 40, 30, 20, 15, 30, 30, 50, 20, 40, 20, 7, 30, 40, 15, 40, 20, 15, 30, 20, 40, 35, 24, 20, 30, 50, 20, 15, 20, 20, 60, 20, 14, 40, 20, 30, 20, 20, 20, 20, 30, 17, 30, 10, 60, 25, 20, 20, 15, 30, 10, 30, 10, 30, 20, 30, 16, 20, 30, 25, 10, 50, 40, 20, 20, 20, 40, 40, 30, 30, 40, 60, 45, 60, 25, 40, 20, 20, 20, 30, 20, 20, 35, 20, 50, 30, 60, 8, 40, 20, 20, 40, 40, 35, 6, 35, 48, 40, 25, 20, 30, 30, 20, 20, 20, 40, 20, 30, 25, 30, 25, 50, 20, 15, 20, 30, 20, 25, 30, 20, 40, 25, 20, 35, 25, 40, 15, 40, 30, 30, 20, 35, 40, 25, 20, 25, 16, 35, 20, 25, 12, 30, 20, 30, 25, 20, 20, 20, 30, 30, 35, 25, 10, 20, 6, 20), RelapseT = c(1.941, NA, 0.9855, NA, NA, NA, NA, NA, NA, NA, NA, 1.05, NA, 3.9865, 1.043, 3.0555, 0.977, 2.976, NA, 0.987, NA, 0.9405, 1.035, 1.0335, NA, NA, NA, 0.9375, NA, 1.013, 2.0135, NA, NA, 0.958, NA, 1.955, NA, NA, NA, NA, NA, NA, NA, NA, 2.048, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.9405, 1.006, 1.9895, NA, NA, NA, NA, NA, NA, 1.073, NA, NA, 1.0745, NA, NA, NA, 2.003, NA, NA, 0.912, NA, NA, NA, NA, 1.9995, NA, NA, NA, NA, 0.98, 3.04, NA, NA, NA, 0.776, NA, 0.9375, NA, NA, NA, NA, NA, NA, 1.944, NA, 2.1, NA, NA, 1.0155, 1.043, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 1.032, 0.9445, NA, NA, 0.995, NA, NA, 0.9635, NA, NA, NA, NA, NA, NA, 0.9625, NA, 1.013, NA, NA, 0.9035, 0.9895, NA, NA, 4.0235, NA, 1.12, 1.017, 1.9645, 1.9715, NA, 0.932, NA, NA, 3.9915, NA, 2.884, NA, 1.046, NA, 2.1615, 3.005, NA, NA, NA, NA, NA, 2.9475, NA, NA, NA, 3.0745, NA, NA, 0.9885, NA, NA, NA, NA, NA, NA, NA, NA, NA, 1.0225, NA, NA, NA, NA, NA, NA, NA, 0.9775, NA, NA, NA, 0.946, 2.0395, NA, 4.0455, NA, NA, NA, NA, NA, NA, 0.9375, NA, NA, NA, 0.9965, 1.0635, 1.005, NA, NA, NA, NA, NA, NA, NA, 2.0575, NA, NA, 1.0035, NA, NA, NA, NA, NA, NA, NA, NA, 1.0145, 0.9785, 1.0265, NA, NA, NA, 0.9695, NA, NA, NA, NA, 2.962, NA, NA, NA, NA, NA, NA, 0.9965, 1.002, 0.9965, NA, NA, 1.0145, 1.017, 0.994, NA, NA, NA, NA, NA, NA, NA), censored.time = c(0, 4.4055, 0, 1.4565, 2.3825, 0.576, 2.549, 0.5255, 4.42, 0.5435, 4.516, 0, 4.4775, 0, 0, 0, 0, 0, 4.4845, 0, 4.382, 0, 0, 0, 1.401, 4.4265, 4.535, 0, 0.556, 0, 0, 1.4905, 3.451, 0, 4.424, 0, 0.4695, 1.4725, 0.4395, 4.5925, 2.5385, 1.602, 2.593, 0.4585, 0, 4.412, 4.427, 3.474, 1.533, 4.471, 4.5245, 0.5545, 4.416, 0.4995, 3.8155, 1.445, 0, 0, 0, 4.4765, 4.6835, 3.5335, 2.474, 3.4915, 3.112, 0, 4.4075, 1.4475, 0, 4.586, 0.49, 4.393, 0, 3.518, 4.582, 0, 4.5435, 0.4585, 4.482, 3.4235, 0, 2.53, 4.563, 4.424, 2.6645, 0, 0, 4.4365, 4.424, 4.669, 0, 0.587, 0, 3.483, 4.575, 3.6785, 3.623, 0.5845, 4.5705, 0, 4.5485, 0, 4.518, 3.5335, 0, 0, 2.6025, 2.587, 0.4585, 2.369, 4.5115, 4.5035, 0.4325, 0.5095, 4.476, 4.3735, 4.4325, 4.418, 1.4135, 0, 0, 4.326, 4.523, 0, 4.53, 4.4665, 0, 4.6675, 3.552, 2.439, 2.5965, 1.4245, 4.505, 0, 0.49, 0, 4.4345, 4.5815, 0, 0, 4.4885, 4.463, 0, 4.461, 0, 0, 0, 0, 4.385, 0, 4.4835, 4.4395, 0, 3.477, 0, 0.4805, 0, 4.4525, 0, 0, 2.618, 4.4165, 0.509, 4.469, 3.5855, 0, 4.6765, 4.5785, 4.531, 0, 4.433, 4.5055, 0, 3.407, 4.515, 4.475, 4.441, 4.404, 3.3635, 4.5175, 4.5705, 4.4235, 0, 4.493, 4.4875, 4.6315, 4.3995, 1.5495, 4.5585, 4.469, 0, 0.4915, 4.538, 4.5965, 0, 0, 3.478, 0, 4.5115, 4.4545, 2.4195, 4.4365, 0.509, 4.578, 0, 4.424, 4.496, 4.393, 0, 0, 0, 4.449, 4.5, 0.4955, 4.5525, 4.5445, 4.472, 4.4805, 0, 4.512, 4.4135, 0, 3.5715, 4.5365, 0.5295, 0.505, 4.509, 4.4225, 1.4225, 3.1475, 0, 0, 0, 0.4535, 4.3865, 4.426, 0, 4.3805, 4.363, 4.436, 4.4385, 0, 4.42, 4.44, 4.553, 0.5735, 4.409, 4.4105, 0, 0, 0, 4.3875, 4.3665, 0, 0, 0, 4.4095, 4.479, 4.415, 4.37, 4.5765, 4.2725, 1.587))