Her kommer en indledning

data information_bias;
input exposure$ outcome$ sex$ freq;
cards;
Yes Case men 240
Yes Case women 140
No Case men 200
No Case women 208
Yes Control men 240
Yes Control women 187
No Control men 600
No Control women 560
;
run;

ods output parameterestimates=information_bias_or2;
proc genmod data=information_bias descending;
class outcome sex;
freq freq;
model exposure=outcome*sex/error=bin link=logit type3;
run;

data information_bias_or2;
set information_bias_or2;
or=exp(estimate);
or_lcl=exp(lowerwaldcl);
or_ucl=exp(upperwaldcl);
if parameter in ('Intercept','Scale') then delete;
format or or_lcl or_ucl 6.2;
keep level1 level2 or or_lcl or_ucl;
run;

ods output parameterestimates=information_bias_or3;
proc genmod data=information_bias descending;
class outcome sex;
freq freq;
model exposure=outcome sex outcome*sex/error=bin link=logit type3;
run;

data information_bias_or3;
set information_bias_or3;
or=exp(estimate);
or_lcl=exp(lowerwaldcl);
or_ucl=exp(upperwaldcl);
if parameter in ('Intercept','Scale') then delete;
format or or_lcl or_ucl 6.2;
keep level1 level2 or or_lcl or_ucl;
run;

ods output parameterestimates=information_bias_or4;
proc genmod data=information_bias descending;
class sex outcome;
freq freq;
model exposure=sex*outcome sex outcome/error=bin link=logit type3;
run;

data information_bias_or4;
set information_bias_or4;
or=exp(estimate);
or_lcl=exp(lowerwaldcl);
or_ucl=exp(upperwaldcl);
if parameter in ('Intercept','Scale') then delete;
format or or_lcl or_ucl 6.2;
keep level1 level2 or or_lcl or_ucl;
run;

ods output parameterestimates=information_bias_or5;
proc genmod data=information_bias descending;
class sex outcome;
freq freq;
model exposure=sex sex*outcome outcome/error=bin link=logit type3;
run;

data information_bias_or5;
set information_bias_or5;
or=exp(estimate);
or_lcl=exp(lowerwaldcl);
or_ucl=exp(upperwaldcl);
if parameter in ('Intercept','Scale') then delete;
format or or_lcl or_ucl 6.2;
keep level1 level2 or or_lcl or_ucl;
run;