* ----------------------------------------------------------------
* ----------------------------------------------------------------
*
* PROJECT: World Championship in Oberstdorf (30 km), Male
* Broadcasted on Thursday, February 12, TV 1 10.55 a.m
*
* Male employment,   1 kv 1987=2,130,500
* Female employment, 1 kv 1987=1,957,500
*
* ----------------------------------------------------------------
* ----------------------------------------------------------------

************************** GLOBALS *******************************
*global femp "12692"     /* Female 20-44 Quart 87, from Statistics Sweden  */
*global memp "13784"     /* Male   20-44 Quart 87, from Statistics Sweden */

*global femp "6881"       /* Female 45-64 Quart 87, from Statistics Sweden */
*global memp "7520"       /* Male   45-64 Quart 87, from Statistics Sweden */
*********************** END GLOBALS ******************************


* PRODUCE GRAPH DATA FOR AGE<44
* ----------------------------------------------------------------
set more off

clear
use sj87_1             /* Population of Sick Cases 1987 */

gen man=(sex==0)       /* Note: sex=1 if female */
gen  week_d=dow(start)
keep if week_d==4      /* Keep only Thursdays */
keep if age<44         /* Keep if age less than 44 */

collapse (count) pnr, by(man start sday smonth)  /* Number of cases per thursday and sex */

by man: gen weeknr=_n            /* Creates a week count week identifier */

reshape wide pnr, i( weeknr) j(man) /* Reshape to wide form */

gen empy0=12692     /* Employment for female aged 20-44 */ 
gen empy1=13784     /* Employment for male aged 20-44 */ 

gen sickr0=pnr0/empy0   /* Female sick rate */
gen sickr1=pnr1/empy1   /* Male sick rate */

gen sickr_young=sickr1/sickr0 /* Male sick rate/Female sick rate */

********* drop weekend and summer periods *********************
drop if weeknr==1 /* nyarsdagen */
drop if weeknr==2 /* tisd trettd jul */
drop if weeknr==16 | weeknr==17 | weeknr==18 /* langfredag + pasklov + 1 maj (fred) */
drop if weeknr==22 /* kristihim. */
drop if weeknr>=23 & weeknr<=35 /* June-August */
drop if weeknr>=49 /* Christmas 17,24,31 dec */

sort weeknr 
drop weeknr
gen weeknr=_n /* Creates a new count week variable ranging from 1-29 */

* Oberstdorf = 5 (weeknr=1,2,...,29)

lab var weeknr "Thursdays 1987"
lab var sickr_y "Male sick rate/Female sick rate, Age 20-44"

graph sickr_y weeknr, c(l) xline(5) ylabel(0.35,0.40,0.45,0.50) t1("Figure 1. Thursday Male/Female Sick-rates, Oberstdorf, x=5") xlabel(1,3,5,7,9,11,13,15,17,19,21,23,25,27,29) saving(Oberst_y, replace) gap(3)
save ographdata_y, replace


* PRODUCE GRAPH DATA FOR AGE>=45
* ----------------------------------------------------------------
set more off

clear
use sj87_1 /* Population of Sick Cases 1987 */

gen man=(sex==0) /* sex=1 if female */
gen  week_d=dow(start)
keep if week_d==4  /* Keep only Thursdays */
keep if age>=45    /* Keep if age less than or equal to 45 */

collapse (count) pnr, by(man start sday smonth)  /* Number of cases per thursday and sex */

by man: gen weeknr=_n         /* Creates a week count week identifier */

reshape wide pnr, i( weeknr) j(man) /* Reshape to wide form */

gen empo0=6881 /* Employment for female aged 45-64 */ 
gen empo1=7520 /* Employment for male aged 45-64 */ 

gen sickr0=pnr0/empo0 /* Female sick rate */
gen sickr1=pnr1/empo1 /* Male sick rate */

gen sickr_o=sickr1/sickr0  /* Male sick rate/Female sick rate */


********* drop weekend and summer periods *********************
drop if weeknr==1 /* nyarsdagen */
drop if weeknr==2 /* tisd trettd jul */
drop if weeknr==16 | weeknr==17 | weeknr==18 /* langfredag + pasklov + 1 maj (fred) */
drop if weeknr==22 /* kristihim. */
drop if weeknr>=23 & weeknr<=35 /* June-August */
drop if weeknr>=49 /* Christmas 17,24,31 dec */

sort weeknr 
drop weeknr
gen weeknr=_n /* Creates a new count week variable ranging from 1-29 */

* Oberstdorf = 5 (weeknr=1,2,...,29)

lab var weeknr "Thursdays 1987"
lab var sickr_o "Male sick rate/Female sick rate, Age 45-64"

graph sickr_o weeknr, c(l) xline(5) ylabel(0.35,0.40,0.45,0.50) t1("Figure 1. Thursday Male/Female Sick-rates, Oberstdorf, x=5") xlabel(1,3,5,7,9,11,13,15,17,19,21,23,25,27,29) saving(Oberst_o, replace) gap(3)
save ographdata_y, replace



* PRODUCE GRAPH DATA FOR ALL AGES (for explanations, see above)
* ----------------------------------------------------------------
set more off

clear
use sj87_1 

gen man=(sex==0) /* sex=1 if female */
gen  week_d=dow(start)
keep if week_d==4  /* KEEP THURSDAYS */

collapse (count) pnr, by(man start sday smonth)  /* NR. SICK CASES PER THURSDAY AND SEX */

sort man start
by man: gen weeknr=_n            /* GEN COUNT WEEK VARIABLE */

reshape wide pnr, i( weeknr) j(man) /* RESHAPE (WIDE) */

gen emp0=19575
gen emp1=21305

gen sickr0=pnr0/emp0 /*Female employment, 1 kv 1987=1,957,500 */
gen sickr1=pnr1/emp1 /*Male employment,   1 kv 1987=2,130,500 */

gen sickr=sickr1/sickr0

drop if weeknr==1 /* nyarsdagen */
drop if weeknr==2 /* tisd trettd jul */
drop if weeknr==16 | weeknr==17 | weeknr==18 /* langfredag + pasklov + 1 maj (fred) */
drop if weeknr==22 /* kristihim. */
drop if weeknr>=23 & weeknr<=35 /* June-August */
drop if weeknr>=49 /* Christmas 17,24,31 dec */

sort weeknr 
drop weeknr
gen weeknr=_n /* Creates a new count week variable ranging from 1-29 */

* Oberstdorf = 5 (weeknr=1,2,...,29)

lab var weeknr "Thursdays 1987"
lab var sickr "Male sick rate/Female sick rate"

graph sickr weeknr, c(l) xline(5) ylabel(0.35,0.40,0.45,0.50) t1("Figure 1. Thursday Male/Female Sick-rates, Oberstdorf, x=5") xlabel(1,3,5,7,9,11,13,15,17,19,21,23,25,27,29) saving(Oberst, replace) gap(3)

translate Oberst.gph Oberst.ps, replace
save ographdata, replace

* ----------------------------------------------------------------
* Reshape from wide to long in order to perform grouped 
* probit estimations
* ----------------------------------------------------------------
reshape long pnr emp, i(weeknr) j(man) /* Reshape from wide to long */

