Dear Ksharp, thanks for the tips, but I failed to apply. Below is the entire codes. I'd like to RMS's to appear in the work library. The below case has 10 replications, this means that I will have 10 RMS values for each replication. Thanks a lot... Serpil data one; call streaminit(156789); do by = 1 to 10; do x2 = 1 to 30; x1 = 30 * ranuni( 512); y = x1 + 2*x2 + rand('chisquare',5); output; end; end; run; proc entropy data=one GME outest=parm1; model y = x1 x2; by; run;
... View more
