Hello,
I would appreciate if someone could help me with the SAS code to count the number of cases(ca case) and controls (pop cont) not exposed to any of the the 4 agents (a1, a2,a3 and a4 ) below. Exposed is 1 and unexposed is zero (0). From the dataset, ids os3, os6, os7 and os9 are unexposed to any of the agents so the freq should be 4. I would like the frequency to be split the frequence into ca case and pop cont. os3: pop cont, os6: ca case, o7: ca case and os9: ca case so in total there are 1 pop cont and 3 ca case (total 4 subjects) but I would like SAS to split the freq of 4 into their respective lung cancer type(ca case and pop cont). The output from attached table 2 is not correct. I need SAS to show and count only the unexposed subjects.
Please, can anyone help me to with the SAS code to solve this problem? My code and log are found below; results are attached, but actually Table 2 is incorrect.
Thank you.
ak.
data agents_exp;
input id$ a1 a2 a3 a4 lung$ 13-20;
datalines;
os1 1 0 0 1 ca case
os2 1 1 0 0 ca case
os3 0 0 0 0 pop cont
os5 1 0 0 1 pop cont
os6 0 0 0 0 ca case
OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
72
73 data agents_exp;
74 input id$ a1 a2 a3 a4 lung$ 13-20;
75 datalines;
NOTE: The data set WORK.AGENTS_EXP has 8 observations and 6 variables.
NOTE: DATA statement used (Total process time):
real time 0.02 seconds
cpu time 0.01 seconds
84 ;
85 run;
86 proc print data=agents_exp;
87 Title 'Table 1: Exposure of ids to 4 agents';
88
NOTE: There were 8 observations read from the data set WORK.AGENTS_EXP.
NOTE: PROCEDURE PRINT used (Total process time):
real time 0.23 seconds
cpu time 0.22 seconds
89 proc sql;
90 create table t as
91 select
92 id, a1, a2, a3,a4,lung,
93 sum(a1,a2,a3,a4)=0 as afinal
94 from agents_exp
95 ;
NOTE: Table WORK.T created, with 8 rows and 7 columns.
96 quit;
NOTE: PROCEDURE SQL used (Total process time):
real time 0.01 seconds
cpu time 0.01 seconds
97
98 proc freq data=t;
99 tables lung*afinal;
100 title 'Table 2: freq of ca case and pop cont for for unexposed agents';
101 run;
NOTE: There were 8 observations read from the data set WORK.T.
NOTE: PROCEDURE FREQ used (Total process time):
real time 0.21 seconds
cpu time 0.20 seconds
102
103 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
115
os7 0 0 0 0 ca case
os8 1 0 1 1 pop cont
os9 0 0 0 0 ca case
;
run;
proc print data=agents_exp;
Title 'Table 1: Exposure of ids to 4 agents';
proc sql;
create table t as
select
id, a1, a2, a3,a4,lung,
sum(a1,a2,a3,a4)=0 as afinal
from agents_exp
;
quit;
proc freq data=t;
tables lung*afinal;
title 'Table 2: freq of ca case and pop cont for for unexposed agents';
run;
