**********************************************;
*f=sasmanual\mb_quest.sas                    *;
*                                            *;
*This program assumes that the temporary SAS *;
*data set called "test" has been created     *;
*from the simulated MB Health raw data.      *;
**********************************************;

options ls=min;

*--------------QUESTION 1a---------------;

proc print data=test (obs=20);
  var gender age los op01 diag01 diag02;
title1 'SIMULATED MB HEALTH DATA FOR SAS COURSE';
title2 'Question 1a';
run;

*--------------QUESTION 1b---------------;

proc sort data=test;
  by gender regionre;
run;

proc print data=test (obs=40);
  var ncase gender regionre icd17brk;
format gender $genderl. regionre $regionl.
       icd17brk $icd17l.;
title2 'Question 1b';
run;

*--------------QUESTION 2---------------;

proc format;
  value $wpgf  '8' = '1'
             other = '0';
  value $wpgl '1' = '1.Winnipeg'
              '0' = '0.non-Wpg';
run;

*--------------QUESTION 3a---------------;

proc means data=test n mean min max maxdec=2;
  var age los deathsep;
title2 'Question 3a';
run;

*--------------QUESTION 3b---------------;

proc means data=test n mean min max maxdec=2;
  var deathsep;
  where deathsep<9999;
title2 'Question 3b';
run;

*--------------QUESTION 3c---------------;

proc means data=test n mean min max maxdec=2;
  var age los;
  class regionre;
format regionre $regionl.;
title2 'Question 3c';
run;

*--------------QUESTION 4a---------------;

proc freq data=test;
  tables icd17brk * gender;
title2 'Question 4a';
run;

*--------------QUESTION 4b---------------;

proc freq data=test;
  tables icd17brk * gender;
  format icd17brk $icd17l. gender $genderl.;
title2 'Question 4b';
run;

*--------------QUESTION 5a---------------;

proc freq data=test;
  tables incdr * charyes;
  format incdr $incdrl. charyes $charl.;
title2 'Question 5a';
run;

*--------------QUESTION 5b---------------;

proc freq data=test;
  tables gender * incdr * charyes;
  format incdr $incdrl. charyes $charl.
         gender $genderl.;
title2 'Question 5b';
run;

*--------------QUESTION 6---------------;

proc format;
  value losf 0-30  = '1'
            31-365 = '2'
          366-high = '3';

  value $losl '1' = '0-30 days'
              '2' = '31-365 days'
              '3' = '366+ days';
run;

data test;
  set test;

loswks = (los/7);

losgroup = put (los,losf.);

wpgres = put (regionre,$wpgf.);

diag3x = substr(diag01,1,3);

op2x = substr(op01,1,2);

label loswks   = 'Length of stay in weeks'
      losgroup = 'LOS grouped into 3 categories'
      wpgres   = 'Winnipeg residence'
      diag3x   = '3-digit diagnosis'
      op2x     = '2-digit procedure';

run;

*--------------QUESTION 7---------------;

proc freq data=test;
  tables los * losgroup/list missing;
  tables regionre * wpgres /list missing;
format losgroup $losl. regionre $regionl. wpgres $wpgl.;
title2 'Question 7a: PROC FREQ check';
run;

proc means data=test n mean min max maxdec=2;
  var los loswks;
title2 'Question 7b: PROC MEANS check';
run;

proc print data=test (obs=30);
  var diag01 diag3x op01 op2x;
title2 'Question 7c: PROC PRINT check';
run;

proc contents data=test;
title2 'Question 7d: PROC CONTENTS check';
run;