Thank you Reeza for proposing the alternative solution, but would you please care to elaborate why is that: "BTW the way you're doing this would only leave you with the results from the last run, not all runs."? Because when I am looking at the output in the SAS window, I see the results corresponding to the specific variable during each iteration. And the estimates are also different for each loop. ... View more

amazing...this is far more elaborate and easy to understand -- thanks a ton for providing a great example ... View more

Okay, great advice, thank you, Tom...I will remove those two lines to make it more efficient ... View more

So the way I understood it:   data new_Results; set &dat.; proc genmod data=&dat.;   the new_Results is created based on the dataset parameter (dat) as defined in the original macro;   and then this, in turn, is being used as input dataset for PROC GENMOD   to generate the parameter estimates which is stored in the data new_Results  ... View more

So here is the updated macro I wrote.  It fixes the problem and generates the output I was looking for. Thank you all for providing valuable feedback and input to get me to write this one.   %macro rate1(dat,vars); %do i = 1 %to %sysfunc(countw(&vars.)); %put %scan(&vars., &i.); *%LET newv=%scan(&vars., &i.); data new_Results; set &dat.; proc genmod data=&dat.; model %scan(&vars., &i.) = / offset=ln dist=poisson lrci; estimate 'Mean' intercept 1; title " cases for in &newv"; run; data new_Results2; set &dat.; IR=%scan(&vars., &i.)/n; IR_1K=IR*1000; LC=quantile('chisq',0.025, %scan(&vars., &i.)*2)/(n*2); UC=quantile('chisq',0.975, (%scan(&vars., &i.)+1)*2)/(n*2); res_l=LC*1000; res_u=UC*1000; run; proc print data=new_Results2; id n %scan(&vars., &i.); title "rate per 1K another way"; run; %END; %mend rate1; ... View more

thank you for clarifying --you are absolutely right. In fact, it is because you pointed it out, I was able to find the bug. I am posting the updated code below. It fixes the problem. Thanks again! ... View more

When I run as a DATA step without using macro, it gives out the specific estimates. For example, instead of using "inputval", if I directly use the "c_tot") with the data=new_Results, it has those estimates. I wonder how can I access the value of "c_tot" by referencing it as a macro variable and perform the arithmetic operations? I guess this is a better way to rephrase my original question. ... View more

As a count data, I have the parameter estimates and without running the MACRO, in a DATA with PROC GENMOD, it is running without giving any error.  ... View more

My mistake: please see the code below:   %macro rate1(dat,vars); %do i = 1 %to %sysfunc(countw(&vars.)); %put %scan(&vars., &i.); %LET newv=%scan(&vars., &i.); data new_Results; set &dat.; proc genmod data=&dat.; model %scan(&vars., &i.) = / offset=ln dist=poisson lrci; estimate 'Mean' intercept 1; title " cases for in &newv"; run; proc print data=new_Results; id n %scan(&vars., &i.); inputval=INPUT(SYMGET(%scan(&vars., &i.)); IR=inputval/n; IR_1K=IR*1000; LC=quantile('chisq',0.025, inputval*2)/(n*2); UC=quantile('chisq',0.975, (inputval+1)*2)/(n*2); res_l=LC*1000; res_u=UC*1000; title "rate per 1K another way"; run; %END; %mend rate1; %rate1(x_trans,c_tot); ... View more

This is the code I wrote.   %macro rate1(dat,vars); %do i = 1 %to %sysfunc(countw(&vars.)); %put %scan(&vars., &i.); %LET newv=%scan(&vars., &i.); data new_Results; set &dat.; proc genmod data=&dat.; model %scan(&vars., &i.) = / offset=ln dist=poisson lrci; estimate 'Mean' intercept 1; title " cases for in &newv"; run; proc print data=new_Results; id n %scan(&vars., &i.); inputval=INPUT(SYMGET(%scan(&vars., &i.)); IR=inputval/n; IR_1K=IR*1000; LC=quantile('chisq',0.025, inputval*2)/(n*2); UC=quantile('chisq',0.975, (inputval+1)*2)/(n*2); res_l=LC*1000; res_u=UC*1000; title "rate per 1K another way"; run; %END; %mend rate1; %run_pr(x_trans,c_tot); ... View more

872 data x_trans; 873 infile cards; 874 input n c_tot c_etopfa c_sb c_lb; 875 ln=log(n); 876 dummy=1; 877 datalines; NOTE: The data set WORK.X_TRANS has 1 observations and 7 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds 879 ; 880 run; 881 %macro rate1(dat,vars); 882 %do i = 1 %to %sysfunc(countw(&vars.)); 883 %put %scan(&vars., &i.); 884 %LET newv=%scan(&vars., &i.); 885 886 data new_Results; 887 set &dat.; 888 proc genmod data=&dat.; 889 model %scan(&vars., &i.) = / offset=ln dist=poisson lrci; 890 estimate 'Mean' intercept 1; 891 title " cases for in &newv"; 892 run; 893 894 895 896 proc print data=new_Results; 897 id n %scan(&vars., &i.); 898 inputval=INPUT(SYMGET(%scan(&vars., &i.)); 899 IR=inputval/n; 900 IR_1K=IR*1000; 901 LC=quantile('chisq',0.025, inputval*2)/(n*2); 902 UC=quantile('chisq',0.975, (inputval+1)*2)/(n*2); 903 res_l=LC*1000; 904 res_u=UC*1000; 905 title "rate per 1K another way"; 906 run; 907 908 %END; 909 910 %mend rate1; 911 %run_pr(x_trans,c_tot); c_tot NOTE: There were 1 observations read from the data set WORK.X_TRANS. NOTE: The data set WORK.NEW_RESULTS has 1 observations and 7 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds NOTE: Writing HTML Body file: sashtml4.htm NOTE: Fitting saturated model. Scale will not be estimated. NOTE: Algorithm converged. NOTE: The scale parameter was held fixed. NOTE: PROCEDURE GENMOD used (Total process time): real time 0.84 seconds cpu time 0.26 seconds ERROR: A character operand was found in the %EVAL function or %IF condition where a numeric operand is required. The condition was: NUVAR ERROR: The macro RUN_PR will stop executing.   This is the code I am getting ... View more

Okay, thank you Tom for pointing it out... I will edit my response. ... View more

okay...thanks a lot-- I will work through the %PUT and see what it does. Much appreciated Reeza for the advice ... View more

okay....so it was the data step during the iteration which I was missing earlier, thank you for explaining. But in the code, it only matches to the index in the array: Afaminc(year) and corresponding faminc value gets repleted. However, how does that get matched across three other columns--faminc96 -faminc98? If SAS matches it internally, then shouldn't it either match by name or by position? But the faminc value does not meet either of them--neither it matches "faminc96" column header nor its position. The same can be said for the other two columns as well. Or is it the SAS convention to replace a value across a row (matched by index) with a specific value? I am probably missing something here. ... View more