求救，SOS，关于重复随机抽取5个变量进行处理的问题？

3149

收藏 2008-02-20

各位武林高手，有两个集合S1、S2，各有48个序列，在两个集合中随机抽样5个序列（如在S1抽取第2、8、17、34、41个序列，那么在S2中也是抽取第2、8、17、34、41序列，两两序列相减，如S1中第2个序列－S2中第2个序列，S1中第2个序列－S1中第8个序列―――，共产生25个序列T1－T25，求每个序列的各个值的平方和，如T1有120个观察值r1-r120，那序列T1的平方和为r1*r2+r2*r2+---+r120*r120，如共有25个值，取最小值。将该过程重复500次，有500个值生成。1．&nbsp;&nbsp; 先生成500个序列 data shuij2．&nbsp;&nbsp; 生成data new1，应有500个观察值，可始终只有1个观察值，请各位帮忙，非常感谢。&nbsp;data shuij;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /* 在1－48 中随机抽取5个数字，不放回抽样，共抽取500次*/&nbsp; array reed(5);&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; */&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; m=5;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; */&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;data shuij(keep=reed1-reed5);&nbsp;&nbsp;&nbsp; /* 在1－48 中随机抽取5个数字，不放回抽样，共抽取500次*/&nbsp;&nbsp; do i=1 to 500;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /* 循环500次 */ &nbsp;&nbsp;&nbsp; &nbsp;array reed(5)；&nbsp;&nbsp;&nbsp; &nbsp;reed(1)=ranuni(13)*48+1;&nbsp;&nbsp; /* 从1－49中随机抽取1个数，可以是小数 */&nbsp;&nbsp; &nbsp;&nbsp;reed(1)=int(reed(1));&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /* 取整数*/&nbsp;&nbsp;&nbsp; &nbsp;j=2;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;do until(j&gt;=6);&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;reed(j)=int((ranuni(13)*48)+1);&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;do k=1 to j-1;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;if reed(j)=reed(k) then go to skip;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;end;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;j=j+1;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;skip: j=j;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /* 再进行4次抽样，如果有重复，不要重复数字，加多一次抽样 */&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;end;&nbsp;&nbsp;&nbsp; output;&nbsp;&nbsp;&nbsp;&nbsp; /*每次循环都输出*/&nbsp;&nbsp;end;run;&nbsp;以上生成5个变量，500个观察值是的数据集shuij，没有问题，可是下面却出了大问题&nbsp;data new1;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; format&nbsp; _all_ 22.16;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; set shuij(drop=time);&nbsp;&nbsp; array reed(5) reed1-reed5;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; if _n_=1 then do until(last);&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;set sss.hmonth;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /* hmonth是一个48个变量，120个观察值的数据集 */&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;set sss.ssmonth end=last;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /* ssmonth是一个48个变量，120个观察值的数据集 */&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;array reh{48} return01-return48;&nbsp;&nbsp;&nbsp;&nbsp; /*将hmonth 48个变量return01-return48 定义为 reh{48}*/&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;array res{48} ret01-ret48;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /*将ssmonth 48个变量ret01-ret48 定义为 res{48}*/&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;array xz(26);&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;array tot(26);&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;f=1;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;do p=1 to 5;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;do q=1 to 5;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;xz(f)=reh{reed(p)}-res{reed(m)};&nbsp; /* 将hmonth和ssmonth按第一次的随机5个数各抽取5个变量，两两相减，得到25变量，即xz1－xz25 */ &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;xz(f)=xz(f)**2;&nbsp;&nbsp;&nbsp; /* 各自平方 */&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;tot(f)+xz(f);&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /* 得到25个变量，每个变量的所有观察值的和tot1－tot25*/&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;f=f+1;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;end;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;end;&nbsp;&nbsp;&nbsp;&nbsp; /* 将按第一次的随机5个数 */&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;end;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;y=20;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;do ac=1 to 25;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;if y&gt;tot(ac) then y=tot(ac);&nbsp;&nbsp; 比较tot1－tot25，取最小值&nbsp;&nbsp;&nbsp; &nbsp;end; run;&nbsp;&nbsp;关于生成new1，我还进行以下测试，都不行，&nbsp;1． data new1;&nbsp;&nbsp; do sbaa=1 to 500;&nbsp;&nbsp; set shuij point=sbaa;&nbsp;&nbsp; &nbsp;array reed(5) reed1-reed5;&nbsp;&nbsp; do sabb=1 to 120;&nbsp;&nbsp; &nbsp;&nbsp;set sss.hmonth point=sabb;&nbsp;&nbsp; &nbsp;&nbsp;set sss.ssmonth point=sabb;&nbsp; &nbsp;&nbsp;&nbsp;array reh{48} return01-return48;&nbsp; &nbsp;&nbsp;&nbsp;array res{48} ret01-ret48;&nbsp;&nbsp; array xz(26);&nbsp;&nbsp; array tot(26);&nbsp;&nbsp;&nbsp; f=1;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;do p=1 to 5;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;do q=1 to 5;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;xz(f)=reh{reed(p)}-res{reed(q)};&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;xz(f)=xz(f)**2;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;tot(f)+xz(f);&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;f=f+1;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;end;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;end;&nbsp;&nbsp;&nbsp; end;&nbsp;&nbsp;&nbsp; stop;&nbsp;&nbsp;&nbsp; y=20;&nbsp;&nbsp;&nbsp;&nbsp; do ac=1 to 25;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if y&gt;tot(ac) then y=tot(ac);&nbsp;&nbsp;&nbsp;&nbsp; end;output;end;stop;run; &nbsp;&nbsp;2．data new1;&nbsp;&nbsp; do sbaa=1 to 500;&nbsp;&nbsp; set shuij point=sbaa;&nbsp;&nbsp; &nbsp;array reed(5) reed1-reed5;&nbsp;&nbsp; if _n_=1 then do until(last);&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;set sss.hmonth (drop=time) end=last;;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;set sss.ssmonth (drop=time) ;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;array reh{48} return01-return48;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;array res{48} ret01-ret48;&nbsp;&nbsp; array xz(26);&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;array tot(26);&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;f=1;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;do p=1 to 5;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;do q=1 to 5;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;xz(f)=reh{reed(p)}-res{reed(q)};&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;xz(f)=xz(f)**2;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;tot(f)+xz(f);&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;f=f+1;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;end;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;end;&nbsp;&nbsp;&nbsp; end;&nbsp;&nbsp; &nbsp;y=20;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;do ac=1 to 25;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;if y&gt;tot(ac) then y=tot(ac);&nbsp;&nbsp;&nbsp;&nbsp; end;output;end;stop;run;

扫码加我拉你入群

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

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

全部回复

patriciabi

2008-2-20 22:56:00

193556.rar
大小:(85.94 KB)

马上下载

附件中两个序列hmonth和ssmonth，各位如果测试成功，将我的错误原因和测试好程序发给我，不胜感谢。

扫码加我拉你入群

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

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

patriciabi

2008-2-20 22:58:00

data shuij;         /* 在1－48 中随机抽取5个数字，不放回抽样，共抽取500次*/  array reed(5);          /*      */        m=5;                /*      */      这个内容重复了，不要

扫码加我拉你入群

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

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

lili83xu

2008-2-24 02:09:00

基于你的思路可修改如下，但不是很有效，运行时间长。 data shuij;                                                                                                                                                                 array reed(5);                                                                                                                                                              do i=1 to 500;                                                                                                                                                      reed(1)=int(ranuni(13)*48+1);                                                                                                                                       j=2;                                                                                                                                                                do until(j>=6);                                                                                                                                                     reed(j)=int((ranuni(13)*48)+1);                                                                                                                                             do k=1 to j-1;                                                                                                                                                      if reed(j)=reed(k) then go to skip;                                                                                                                                 end;                                                                                                                                                                j=j+1;                                                                                                                                                              skip: j=j;                                                                                                                                                  end;                                                                                                                                                            output;                                                                                                                                                         end;                                                                                                                                                        run;                                                                                                                                                                libname sss 'd:/sas/';                                                                                                                                              data month(drop=time);                                                                                                                                                      merge sss.hmonth sss.ssmonth ; by time;                                                                                                                             run;  %macro test(m);                                                                                                                                                    %do k=1 %to &m;                                                                                                                                                             data a(keep=reed1-reed5);                                                                                                                                              set shuij ; if i=&k;                                                                                                                                          run;                                                                                                                                                                data b&k(keep=y);                                                                                                                                                           set month; if _n_=1 then do until(last); set a end=last; end;                                                                                                       array reh{48} return01-return48;                                                                                                                                    array res{48} ret01-ret48;                                                                                                                                          array reed{5} reed1-reed5;                                                                                                                                          array xz(25);                                                                                                                                                       array tot(25);                                                                                                                                                      f=1;                                                                                                                                                                do p=1 to 5;                                                                                                                                                                do q=1 to 5;                                                                                                                                                        xz(f)=reh{reed(p)}-res{reed(q)};                                                                                                                                    tot(f)+xz(f)*xz(f); f+1;                                                                                                                                            end;                                                                                                                                                        end;                                                                                                                                                                if _n_=120;                                                                                                                                                         y=20;                                                                                                                                                               do ac=1 to 25;                                                                                                                                                               if y > tot(ac) then y=tot(ac);                                                                                                                                end;                                                                                                                                                        run;                                                                                                                                                        %end;                                                                                                                                                                data c;                                                                                                                                                                    set %do k=1 %to &m; b&k %end;                                                                                                                                        ;                                                                                                                                                              run;                                                                                                                                                                 proc datasets lib=work;                                                                                                                                                delete %do k=1 %to &m; b&k %end; / memtype=data;                                                                                                                   run;                                                                                                                                                                 quit;                                                                                                                                                                %mend test(m);                                                                                                                                                       %test(500); /***开始可用10或50测试下，用500运行时间较长。***/

[此贴子已经被作者于2008-2-24 4:19:23编辑过]

扫码加我拉你入群

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

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

lili83xu

2008-2-24 04:29:00

之前的问题在于：set shuij;                                                                                                                                                                     array reed(5) reed1-reed5;                                                                                                                                                   if _n_=1 then do until(last);                                                                                                                                                   set sss.hmonth;                                                                                                                                                             set sss.ssmonth end=last; 运行下就可发现数据集只包含hmonth和ssmonth在shuiji第500次的对应值。

扫码加我拉你入群

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

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

lili83xu

2008-2-24 04:34:00

也可以这样做，应该会快些。%macro merge(no1,no2);                                                                                                                                                      data shuij2(keep=reed1-reed5 i);                                                                                                                                                    set %do k=1 %to &no1; shuij %end;                                                                                                                                           ;                                                                                                                                                                   run;                                                                                                                                                                        data month2(drop=time);                                                                                                                                                             set %do k=1 %to &no2; month %end;                                                                                                                                         ;                                                                                                                                                                     %mend merge(no1,no2);                                                                                                                                                       %merge(120,500);                                                                                                                                                                                                                                                                                                                                        proc sort data=shuij2;                                                                                                                                                              by i;                                                                                                                                                               run;                                                                                                                                                                                                                                                                                                                                                    data test(keep=y);                                                                                                                                                                  merge shuij2 month2;                                                                                                                                                        array reh{48} return01-return48;                                                                                                                                            array res{48} ret01-ret48;                                                                                                                                                  array reed{5} reed1-reed5;                                                                                                                                                  array xz(25);                                                                                                                                                               array tot(25);                                                                                                                                                              f=1;                                                                                                                                                                        do p=1 to 5;                                                                                                                                                                        do q=1 to 5;                                                                                                                                                                xz(f)=reh{reed(p)}-res{reed(q)};                                                                                                                                            if mod(_n_-1,120)=0  then                                                                                                                                                           do;                                                                                                                                                                         tot(f)=xz(f)*xz(f);                                                                                                                                                         end;                                                                                                                                                                else                                                                                                                                                                                do;                                                                                                                                                                         tot(f)+xz(f)*xz(f);                                                                                                                                                         end;                                                                                                                                                                f+1;                                                                                                                                                                        end;                                                                                                                                                                end;                                                                                                                                                                      keep=tot1-tot25;                                                                                                                                                            if mod(_n_,120)=0;                                                                                                                                                            y=20;                                                                                                                                                                       do a=1 to 25;                                                                                                                                                                       if y > tot(a) then y=tot(a);                                                                                                                                                end;                                                                                                                                                        run;

扫码加我拉你入群

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

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

点击查看更多内容…

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群