/*
Number of in-hospital deaths at Edouard Herriot Hospital ,
and Daily Maximal Temperatures during summers of 2002 and 2003,
Lyon, France.

Source:
*/

options ls= 85 ps=50; run;

data a;
input year month day temp_C deaths;
T_gt_37 = (temp_C > 37);
year2_3 = (year - 2000);
year0_1 = (year - 2002);
aug_day = day; if month = 7 then aug_day = day - 31;
lines;
2002 7 26 27 5
2002 7 27 29 3
2002 7 28 31 7
2002 7 29 32 7
2002 7 30 29 3
2002 7 31 24 1
2002 8  1 21 2
2002 8  2 26 2
2002 8  3 21 5
2002 8  4 23 2
2002 8  5 26 0
2002 8  6 22 0
2002 8  7 24 4
2002 8  8 26 2
2002 8  9 21 1
2002 8 10 20 3
2002 8 11 19 2
2002 8 12 22 2
2002 8 13 25 4
2002 8 14 26 3
2002 8 15 29 3
2002 8 16 30 1
2002 8 17 30 1
2002 8 18 31 0
2002 8 19 32 1
2003 7 26 34 4
2003 7 27 23 2
2003 7 28 26 7
2003 7 29 28 2
2003 7 30 29 4
2003 7 31 24 3
2003 8  1 31 2
2003 8  2 33 5
2003 8  3 37 3
2003 8  4 39 3
2003 8  5 39 4
2003 8  6 39 6
2003 8  7 38 5
2003 8  8 38 9
2003 8  9 38 9
2003 8 10 39 11
2003 8 11 39 9
2003 8 12 40 13
2003 8 13 40 5
2003 8 14 33 4
2003 8 15 28 1
2003 8 16 31 2
2003 8 17 31 2
2003 8 18 30 2
2003 8 19 32 3
;
run;

proc plot data=a ;
 plot deaths*aug_day=year2_3 temp_c*aug_day=year2_3 / overlay;

proc means data=a n mean std; var deaths;  run;

proc means data=a n mean std;
class  year; var deaths; run;

title 2002; run;
proc genmod data=a ;
model deaths = / dist = poisson
                 link = identity;
where (year = 2002);
run;

proc genmod data=a;
model deaths = / dist = normal
                 link = identity;
where (year = 2002);
run;

proc genmod data=a;
model deaths = year0_1 / dist = normal
                        link = identity;
run;


title 2003; run;
proc genmod data=a;
model deaths = / dist = poisson
                 link = log;
where (year = 2003);

proc means data=a n mean std;
class  T_gt_37; var deaths;
where (year = 2003); run;

title 2003; run;
proc genmod data=a;
model deaths = T_gt_37 / dist = poisson
                 link = identity;
where (year = 2003);
run;

title 2003; run;
proc genmod data=a;
model deaths = T_gt_37 / dist = normal
                 link = identity;
where (year = 2003);
run;