/*----------------------------------------------------------------
ex3.sas  Code to run exemplar 3
Written by Gillian Raab
 November 2004

proc contents data=exemp3.ex3; run;
first tell sas where the data for this exemplar is stored
  Change this to the directory  where YOU have the data
-----------------------------------------------------------------*/
libname exemp3 'j:\depts\fhls\peas\ex3datafiles\data' ;

*libname exemp3 'c:\documents and settings\gillian raab\my documents\aprojects\peas\web\exemp3\data\' ;
/*-------------------------------------------------
NOW open and run the program 

ex3_formats.sas

YOU WILL NEED THIS FOR THE NEXT BITS TO RUN CORRECTLY

Now use surveymeans to get standard errrors of proportions
of different smoking categories

NOTICE THAT THE CLASS STATEMENT TELLS SAS TO DO PROPORTIONS
----------------------------------------------------*/
title 'PROPORTION SMOKING';
proc surveymeans data=exemp3.ex3 mean stderr nobs clm; * no strata;
weight WEIGHTA;
cluster psu;
CLASS CIGST1 ;
FORMAT CIGST1 CIGSTA.;
var CIGST1;
strata stratum;
run;
/*-------------------------------------------------------
nOW RECODE THIS TO GET A VARIABLE FOR SMOKERS
-----------------------------------------------------------------------*/
DATA EX3; * NOTE THIS IS TEMP FILE ;
SET EXEMP3.EX3;
if cigst1=4 then smoker=1;
else if cigst1=. then smoker=.;
else smoker=0;
if nofad>5 then nofad=5; * recode all large adults >5 to 5;
run;
title 'Proportion smoking by number of adults';
proc surveymeans data=ex3 mean stderr nobs clm; * smoking by number of adults;
weight WEIGHTA;
cluster psu;
CLASS smoker ;
domain nofad;
var smoker;
strata stratum;
run;
title 'Proportion smoking by sex';
proc surveymeans data=ex3 mean stderr nobs clm; * now by sex;
weight WEIGHTA;
cluster psu;
CLASS smoker ;
domain sex;
format sex sex.;
var smoker;
strata stratum;
run;
/*---------------------------------------------------------------
to compare smoking rates by groups need to use regression
--------------------------------------------------------------*/
proc surveyreg data=ex3;
weight WEIGHTA; cluster psu;strata stratum;
model smoker=sex;
format sex sex.;
run;
title 'Proportion smoking by health board';
proc surveymeans data=ex3 mean stderr nobs clm; * now by health board;
weight WEIGHTA;
cluster psu;
CLASS smoker ;
domain hboard;
format hboard hboard.;
var smoker;
strata stratum;
run;
proc surveyreg data=ex3 ; * now a complicated model;
weight WEIGHTA;cluster psu;strata stratum;
CLASS  hboard  ;
format  hboard hboardn.;
model smoker= hboard  /solution;
estimate 'fife - lothian' hboard  0 0 0 0 1 0 0 0 0 0 0 -1 0 0 0;
run;
title 'Proportion smoking more complicated model';
proc surveyreg data=ex3 ; * now a complicated model;
weight WEIGHTA;cluster psu;strata stratum;
CLASS age sex hboard  ;
format sex sex. age age. hboard hboardn.;
model smoker=nofad age sex hboard age*sex nofad*sex nofad*age ;
run;
title 'Proportion smoking final model';
proc surveyreg data=ex3 ; * now a complicated model;
weight WEIGHTA;cluster psu;strata stratum;
CLASS age sex hboard  ;
format sex sex. age age. hboard hboardn.;
model smoker=nofad age sex hboard /solution ; * this gives estimates;

run;