求助 SAS模拟双峰分布密度曲线及其抽样分布？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问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);

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：distribution histogram proc sql footnote Sampling SAS 双峰分布 抽样分布

你为什么要求2000个样本均值？

建议看看混合模型的经典著作：Finite Mixture Models_[McLachlanPeel2001]，第一章就说明了混合模型的涵义