*****************************************
* Clinical Data Set Questions           *
* V. Adding Variables and Observations  *
*****************************************
options linesize=min;
libname mydir 'c:\My Documents\My SAS Files\sasmanual\data'; 

data nopreg;
set mydir.clinical;
if pregnant=0 and gender='F' then output;
run;

proc print data=nopreg;
title 'Females that are not pregnant';
run;

data male;
set mydir.clinical;
if gender='M' then output;
run;

proc print data=male;
title 'Males';
run;

data concat;
  set nopreg (in=m1)
        male (in=m2);
inone=m1;
intwo=m2;
run;

proc print data=concat;
format prim_dx sec_dx $dxcodes.
       dob visit mmddyy8.
       gender $gender.
       vitamins pregnant $yesno.;
title 'All that are not pregnant concatenated'; 
run;

proc sort data=male;
by dob;
run;

proc sort data=nopreg;
by dob;
run;

data interleave;
set nopreg
	male;
by dob;
run;

proc print data=interleave;
format prim_dx sec_dx $dxcodes.
       dob visit mmddyy8.
       gender $gender.
       vitamins pregnant $yesno.; 
title 'All that are not pregnant interleaved by dob';
run;
title;

data one;
set mydir.clinical (keep=id gender dob prim_dx sec_dx);
run;

data two;
set mydir.clinical (keep=gender hr);
run;

proc print data=one;
title 'Data=One';
run;

proc print data=two;
title 'Data=Two';
run;

proc sort data=one;
by gender;
run;

proc sort data=two;
by gender;
run;

data mer;
merge one (in=m1)
      two (in=m2);
inone=m1;
intwo=m2;
run;

proc print data=mer;
title 'Adding hr variable using Merge';
run;

data mer;
merge one (in=m1)
      two (in=m2);
by gender;
if hr>70 then output;
run;

proc print data=mer;
title 'All patients with heart rate over 70';
run;

proc means data=mydir.clinical noprint nway;
class gender;
var hr;
output out=mhr mean=mean_hr;
run;

proc print data=mhr;
title 'Mean heart rate by gender';
run;

proc sort data=mydir.clinical out=clin;
by gender;
run;

proc sort data=mhr;
by gender;
run;

data mer; 
merge clin (in=m1)
      mhr (in=m2);
by gender;
if hr>mean_hr then hi_hr=1;
else hi_hr=0;
label hi_hr='Heart Rate Greater than the Mean value for each sex';
run;

proc print data=mer;
title 'Data mer - adding mean heart rate by gender';
run;

proc means data=mer n sum mean;
title 'Proportion of observation with Heart rate greater than the mean heart rate';
class gender;
var hi_hr;
run;