data eq; input ID y x z w; cards; 1 1 65 40 0 2 1 34 . 0 3 0 . 23 0 4 1 48 35 0 ; run; %macro em (missper, n_imp, my_seed); proc mi data=eq seed=&&my_seed&i simple nimpute=&n_imp out=EMIMP&missper ; em itprint ; var y x z w; run; %mend em; %macro mi (missper,n_imp,n_seed); %em (&missper,&n_imp,&n_seed); %mend mi; %macro large(missper,n_imp,n_seed,n_boot); data seed_long(drop=i); call streaminit(&n_seed); do i = 1 to &n_boot; x_seed = floor(rand("Uniform")*500000); output; end; run; data _null_; set seed_long; call symput('my_seed'|| compress(put(_n_,4.)), x_seed); run; %let i=1; %do %until (&i>&n_boot) ; %mi(&missper) %let i = %eval(&i+1); %end; %mend large; %large(20, 5, 3456,3);  Hello,  What I am trying to do is using a random fixed seed to create EM output. but I got the following error.   NOTE 137-205: Line generated by the invoked macro "EM". 1 proc mi data=eq seed=&&my_seed&i simple nimpute=&n_imp out=EMIMP&missper ; em --- 22 1 ! itprint ; var y x z w; run; ERROR 22-322: Syntax error, expecting one of the following: an integer constant, PCTMISSING.     I know %large has no problem. The problem is the writing of EM and MI macro variable. But I don't know what is the right way to do it.   Thanks ... View more

data eq1; input ID y x z w; cards; 1 1 30 . 25 3 1 . 45 . 4 1 48 25 56 ; run; %macro genmod(method,missper); proc genmod data = &&method&missper; class y(ref='0') id; model y= x z w/ dist=bin link=logit pred covb; repeated subject=id / type=ar(1) PRINTMLE; ods output ParameterEstimates=&&method.lgs&missper; run; %MEND genmod; %genmod(eq,1); Hello, I am writing some MACRO. I know for above code example ( eq=method and missper=1 ), if I want to write 'eq1', I could use &&method&missper. But I want to know more complex writing way. For 'eqlgs1', I use &&method.lgs&missper. what about 'leq.gs1' or 'leq_1gs'? something like those. Any other way to write eqlgs1?    Thanks ... View more

Hi all, I made some codes including macro. They went well when I used last week, but I always found errors when using them this week %macro try(missper,n_boot); %MACRO mi(missper); data rmse&missper; set rmse&missper; index = &i; run; %if &i=1 %then %do; data poolrmseem&missper rmse_old rmse&missper; set rmse&missper; index = &i; run; %end; %else %if &i>1 %then %do; data poolrmseem&missper;set rmse_old rmse&missper;run; data rmse_old;set poolrmseem&missper;run; %end; %mend mi; %mi(&missper); %let i=1; %do %until (&i>&n_boot); %mi(&missper); %let i = %eval(&i+1); %end; %mend try; %try (20, 5); . Pretty weird! I have listed the code and log message below. Could anybody help me to find out the reason? Thanks.       22: LINE and COLUMN cannot be determined. NOTE 242-205: NOSPOOL is on. Rerunning with OPTION SPOOL might allow recovery of the LINE and COLUMN where the error has occurred. ERROR 22-322: Syntax error, expecting one of the following: a name, a quoted string, a numeric constant, a datetime constant, a missing value, INPUT, PUT. WARNING: Apparent symbolic reference I not resolved. NOTE: The SAS System stopped processing this step because of errors. WARNING: The data set WORK.RMSE20 may be incomplete. When this step was stopped there were 0 observations and 5 variables. WARNING: Data set WORK.RMSE20 was not replaced because this step was stopped. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds WARNING: Apparent symbolic reference I not resolved. ERROR: A character operand was found in the %EVAL function or %IF condition where a numeric operand is required. The condition was: &i=1 ERROR: The macro MI will stop executing. ... View more

Hello,   I want to merge the mean and std results only for y and x, created by proc means from dataset eq1, eq2, eq3. And add a new column with dataset name.  Could anyone show me how to do it?  Thanks data eq1; input ID y x z w; cards; 1 1 1 27 40 1 2 3 . 29 1 3 5 30 . 1 4 7 38 38 2 1 1 23 45 2 2 3 32 20 2 3 5 67 . 2 4 7 . 27 3 1 1 33 23 3 2 3 21 12 3 3 5 78 . 3 4 7 13 45 4 1 1 56 45 4 2 3 67 13 4 3 5 . 35 4 4 7 48 35 ; run; data eq2; input ID y x z w; cards; 1 1 27 40 8 1 2 . 29 37 1 3 30 . 25 1 4 38 38 23 2 1 23 45 19 2 2 32 20 . 2 3 67 . . 2 4 . 27 . 3 1 33 23 46 3 2 21 12 56 3 3 78 . 34 3 4 13 45 . 4 1 56 45 23 4 2 67 13 67 4 3 . 35 13 4 4 48 35 56 ; run; data eq3; input y x ; cards; 1 1 2 3 3 5 4 7 ; run; ... View more

Hello,   I want to calculate RMSE and MAE. With data data4rmse I could get one RMSE and MAE using following code.   Now I have data EQ with repeating the experiment for 3 times. How could I get RMSE and MAE for each repeat and output them into one table?    Thanks data data4rmse; input ID bmi dm bmi_me; cards; 1 13 40 8 2 13 29 37 3 30 14 25 4 38 38 23 ; run; data rmse; retain square_error_sum abs_error_sum; set data4rmse end=last /* Flag for the last observation */ ; error = BMI - BMI_me; /* Calculate simple error */ square_error = error * error; /* error^2 */ if _n_ eq 1 then do; /* Initialize the sums */ square_error_sum = square_error; abs_error_sum = abs(error); end; else do; /* Add to the sum */ square_error_sum = square_error_sum + square_error; abs_error_sum = abs_error_sum + abs(error); end; if last then do; /* Calculate RMSE and MAE and store in SAS data set. */ mae = abs_error_sum/_n_; rmse = sqrt(square_error_sum/_n_); end; run; data eq; input repeat ID bmi dm bmi_me; cards; 1 1 13 40 8 1 2 13 29 37 1 3 30 14 25 1 4 38 38 23 2 1 30 45 19 2 2 30 20 14 2 3 67 13 14 2 4 13 27 13 3 1 33 23 46 3 2 13 12 56 3 3 13 13 34 3 4 13 45 13 ; run; ... View more