%macro seedgen(n);
data pseudo(drop=i seed);
seed=1;
do i=1 to &n;
call rannor(seed,x);
output;
end;
proc univariate noprint data=pseudo;
var x;
output out=test normal=w probt=p;
run;
%mend;
上面是产生随机数的一个宏,没产生一组随机数会计算P 和 W 两个统计量。 我现在想得到100组W和P的值,然后分别画histograms,code该怎么写,本人菜鸟。下面是我写的code,但跑不出结果。求大神帮助。
data alltest;
input w p;
do i=1 to 100;
%seedgen(25)
w=w;
p=p;
output;
end;
run;
谢谢!