data survey;
   infile "C:/WINDOWS/Desktop/hs survey.txt";
   do race = 'white','other';
      do gender = 'female', 'male';
        do marijuana = 'yes','no';
           do Alcohol = 'yes','no';
              do cigarette = 'yes','no';
                 input count @;
                 output;
              end;
           end;
        end;
      end;
   end;
run;
proc print data=survey;
   var race -- cigarette;
run;
title 'HS Survey Data';
proc genmod data=survey descending;
   title1 'Mutual independence + GR';
   class cigarette alcohol marijuana race gender;
   model count = race -- cigarette gender*race  / dist=poisson link=log lrci;
run;
proc genmod data=survey descending;
   title1 'Homogenous association';
   class cigarette alcohol marijuana race gender;
   model count = cigarette | alcohol | marijuana | race | gender @2  / dist=poisson link=log lrci;
run;
proc genmod data=survey descending;
   title1 'All three-factor terms';
   class cigarette alcohol marijuana race gender;
   model count = cigarette | alcohol | marijuana | race | gender @3  / dist=poisson link=log lrci;
run;
proc genmod data=survey descending;
   title1 'AC, AM, CM, AG, AR, CG, GM, MR, GR - CR';
   class cigarette alcohol marijuana race gender;
   model count = cigarette alcohol marijuana race gender
                 alcohol*cigarette alcohol*marijuana cigarette*marijuana
                 alcohol*gender alcohol*race cigarette*gender gender*marijuana
                 gender*race marijuana*race/ dist=poisson link=log lrci;
run;
proc genmod data=survey descending;
   title1 'AC, AM, CM, AG, AR, GM, MR, GR - CG';
   class cigarette alcohol marijuana race gender;
   model count = cigarette alcohol marijuana race gender
                 alcohol*cigarette alcohol*marijuana cigarette*marijuana
                 alcohol*gender alcohol*race gender*marijuana gender*race
                 marijuana*race/ dist=poisson link=log lrci;
run;
proc genmod data=survey descending;
   title1 'AC, AM, CM, AG, AR, GM, GR - MR';
   class cigarette alcohol marijuana race gender;
   model count = cigarette alcohol marijuana race gender
                 alcohol*cigarette alcohol*marijuana cigarette*marijuana
                 alcohol*gender alcohol*race gender*marijuana gender*race
                 / dist=poisson link=log lrci;
run;
proc genmod data=survey descending;
   title1 'AC, AM, CM, AG, AR, GR - GM';
   class cigarette alcohol marijuana race gender;
   model count = cigarette alcohol marijuana race gender
                 alcohol*cigarette alcohol*marijuana cigarette*marijuana
                 alcohol*gender alcohol*race gender*race
                 / dist=poisson link=log lrci;
run;
proc genmod data=survey descending;
   title1 'AC, AM, CM';
   class cigarette alcohol marijuana race gender;
   model count = cigarette alcohol marijuana race gender
                 alcohol*cigarette alcohol*marijuana cigarette*marijuana
                 / dist=poisson link=log lrci;
run;