请问SAS是否有直接生成双峰分布密度函数的function? 小白没找到，试着自己理解写了一下，不知道对不对？求大神指导？


%macro DOUBLET(n = );
proc sql noprint;
        create table data_7(MEAN_ num);
quit;


%do i = 1 %to 2000;
        data data_21;
                do i = 1 to &n ;
                        X1 = rand('NORMAL',-8,1);
                        X2 = rand('NORMAL',8,1);
                        K = RAND('BINOMIAL', 0.4, 1) ;
                        output;
                end;
        run;


        data data_4_&i;
                set data_21 ;
                X = K*X1+(1-K)*X2;
        run;


        proc means data = data_4_&i noprint;
                var X;
                output out = data_5_&i mean = MEAN;
        run;


        proc sql noprint;
                insert into data_7
                select MEAN as MEAN_
                from data_5_&i;
        quit;
%end;


data data_8;
        length GROUP $200;
        set data_7;
        GROUP = 'n='||strip(put(&n,best.));
run;


proc sgplot data = data_8;
        histogram MEAN_ / group = GROUP NBINS = 1800;
        xaxis min = -10 max = 10;
        yaxis min = 0 max = 2;
        title ;
        footnote 'Doublet_Sampling Distribution';
run;
%mend;


        %DOUBLET(n = 30);