Hi, I have the same mistake. I want to make a table that gives computed type I error rate from the comparison repeat at 50 times. From the Van der Waerden test p values should be drawn and the type 1 error calculated.How can I get output as a single table. %mend; %macro simu (r=, n1=, n2=, n3=, mu1=, mu2=, mu3=, sigma1=, sigma2=, sigma3=, times=); data Simulation(drop=i); do SampleID = 1 to &Times; do i = 1 to &n1; group=1; y=rand("normal",&mu1,&sigma1);/* Group 1: x ~ N(mu1,sigma1^) */ output; end; do i = 1 to &n2; group=2; y=rand("normal",&mu2,&sigma2); /* Group 2: x ~ N(mu2,sigma2^) */ output; end; do i = 1 to &n3; group=3; y=rand("normal",&mu3,&sigma3); /* Group 3: x ~ N(mu3,sigma2^) */ output; end; end; run; proc npar1way VW data=Simulation; by SampleID; class group; ods graphics off; run; %mend; %simu (n1=5, n2=5, n3=5, mu1=10, mu2=10, mu3=10, sigma1=1, sigma2=1, sigma3=1, times=50); run;
... View more