I have tried to do this two ways. I want to plot two regression lines. I am not sure I have the correct sg procedure for I am not getting the two regression lines I want and I don't know why the legend is not positioning itself in top left corner in proc gplot. No errors are indicated. Any ideas? Thanks. MM   data simpson; input week X Y percentLA; percentRecall=X/Y; datalines; 1 60 200 0 1 320 800 0 2 640 800 100 2 180 200 100 ; proc sgscatter data=simpson; title 'Simpson Paradox'; plot percentRecall*percentLA / group=week; run; data simpson2; input week percentRecall percentLA; cards; 1 .3 0 2 .4 0 1 .8 100 2 .9 100 ; PROC GPLOT DATA=simpson2; PLOT percentRecall*percentLA=week; symbol1 value=dot i=r ci=orange; symbol2 value=dot i=r ci=blue; legend mode=protect down=2 position=(top left); run; quit; ... View more

I am doing a cochran-mantel-haenszel test using proc freq. The output is very confusing. Is there an order to the options included with the tables statement ?  For example, the Breslow-Day statistic is not always output and there are two cmh tables with different values. What are those supposed to be? I expected one cmh table with cmh1 cmh2 and cmh3.   SAS CODE: data b; input confounder $ treatment $ response $ count; cards; 1 drug1 yes 19 1 drug1 no 3 1 drug2 yes 14 1 drug2 no 11 2 drug1 yes 27 2 drug1 no 18 2 drug2 yes 21 2 drug2 no 15 ; proc freq data=b order=data; table treatment*response / nopercent norow nocol chisq relrisk; weight count; run cancel; proc freq data=b order=data; tables treatment*response / nopercent norow nocol relrisk chisq cmh(mf); by confounder; weight count; run; Thanks.   MM ... View more

  Output from proc freq gives a Clopper-Pearson (exact) confidence interval but how do I compute the accompanying statistic of a test of p=po and its pvalue? Also when I am working with single binomial proportion, how is Fishers exact test computed. It is also listed on the output. Where can I find freq documentation about that. I must be using wrong search criteria for I am not seeing the answers for each of these questions. My notes have approximate and exact z statistics varying by denominator as seen in snaphot at beginning of this post.   title " Ho: p=.7 vs Ha: p ne .7"; data CellCounts; input R $ Count; np=40*.7; nq=40*.3; datalines; yes 35 no 5 ; run; proc freq data = CellCounts order=freq; tables R / binomial(p=0.7 exact score) alpha=.05; exact binomial; weight Count; run; quit; Thank you. MM   ... View more