topic How to add stratified variables in proc freq code when finding risk diff between treatments in SAS Programming

How to add stratified variables in proc freq code when finding risk diff between treatments

knveraraju91 — Wed, 21 Dec 2022 14:23:14 GMT

Dear

I am using below proc freq code to find risk diff between treatments. Now I need to calculate stratified riskdiff. The stratification variable are age ss ph. Please suggest me how to update code. Thank you

data one;
input trt res$ age ss ph $;
datalines;
1 cr 55 1 Y
1 cr 53 1 N
1 cri 52 1 Y
1 cri 58 1 N
1 cri 50 2 N
1 cr 53 2 N
1 cri 52 2 Y
1 cri 58 2 N
2 cri 56 2 Y
2 cri 51 2 N
2 cr 50 2 Y
2 cr 56 2 N
;
 proc freq data = one;
                tables trt * res / riskdiff cl alpha = 0.025 ;
                ods output riskdiffcol1 = rd (keep = control risk lowercl uppercl where = (control)) ;
        run ;

Re: How to add stratified variables in proc freq code when finding risk diff between treatments

Ksharp — Thu, 22 Dec 2022 11:45:00 GMT

data one;
input trt res$ age ss ph $;
datalines;
1 cr 55 1 Y
1 cr 53 1 N
1 cri 52 1 Y
1 cri 58 1 N
1 cri 50 2 N
1 cr 53 2 N
1 cri 52 2 Y
1 cri 58 2 N
2 cri 56 2 Y
2 cri 51 2 N
2 cr 50 2 Y
2 cr 56 2 N
;
 proc freq data = one;
                tables ss*ph* trt * res / riskdiff(common cl=wald)  alpha = 0.025 ;
        run ;

Re: How to add stratified variables in proc freq code when finding risk diff between treatments

knveraraju91 — Thu, 22 Dec 2022 17:03:37 GMT

Thank you for the help. Can i also include age in

 tables ss*ph* trt * res / riskdiff(common cl=wald)  alpha = 0.025 ;

tables statement. How to output dataset containing risk diff values. Thank you

Re: How to add stratified variables in proc freq code when finding risk diff between treatments

Ksharp — Fri, 23 Dec 2022 11:56:20 GMT

/*
Sure you can add age in it ( if your data was enough to analysis ).
The following could get what you need.
*/
ods select none;
ods output CommonPdiff=CommonPdiff RiskDiffPlot=RiskDiffPlot;
 proc freq data = one;
                tables age*ss*ph* trt * res / riskdiff(common cl=wald)  alpha = 0.025 ;
        run ;
ods select all;