求助 利用SAS实现 cluster bootstrap

shiyiming

to jingju

I tried your code
but it didn't work, is that available for 9.2 ?? mine is 9.1

here is the log:

   proc surveyselect data =class out =SampleOut reps =10 n =10 method =urs seed =555;
9       samplingunit cluster;
        ------------
        180
ERROR 180-322: 语句无效或未按正确顺序使用。
10   run;


Thanks!

shiyiming

to caicaisas
你的意思似乎是要先对CLUSTER做一个URS，然后给定CLUSTER号码以后，对CLUSTER内部的ID再做URS，对不对？
刚吃早饭，目前只能想到一个很笨的方法，容我想想有没有更有效率的。

shiyiming

to OL:

[quote:15z1w7rs](samplingunit) ....is that available for 9.2 ?? mine is 9.1[/quote:15z1w7rs]

You are right!

[quote:15z1w7rs]Beginning with SAS/STAT 9.22 in SAS 9.2 TS2M3, use the SAMPLINGUNIT or CLUSTER statement to name variables that identify the sampling units as groups of observations (clusters).[/quote:15z1w7rs]

shiyiming

to OLO, 我不需要对cluster内的ID再做URS，但是在merge的时候出现了问题

BY the way , URS 是无放回抽样吗？
其实我要做的是有放回抽样
所以cluster做完之后有重复的值，导致merge的时候就出错了

shiyiming

to caicaisas
URS=unrestricted sampling=有放回抽样
try the following code


data temp;
     do cluster=1 to 10;
            do id=1 to 40;
                   output;
                end;
         end;
run;

proc freq data=temp noprint;
     table cluster/out=cluster(drop=percent);
run;

proc surveyselect data=cluster  method=urs  sampsize=10 out=cluster_samp;
run;

data temp2;
     merge temp  cluster_samp;
         by cluster;
         if NumberHits>=1 then do;
            do seq=1 to NumberHits;   output; end;
         end;
run;

proc sort data=temp2; by cluster  seq; run;

shiyiming

to OLO~

已经试出来了～～非常感谢！！

shiyiming

to oloolo
这段代码确实在功能上实现了cluster bootstrap～
但是关于cluster bootstrap 不知道还有没有改进的可能？程序最后只实现了一个样本，如果我想得到B个bootstrap的样本，就需要套一个宏，将它跑B次，然后把得到的B个数据集append在一起。这样特别不高效，因为我还需要把这个结果套在另一个宏里面，所以今天一直在想如何能不用宏解决多次的cluster bootstrap， 不过真的水平有限，一直没弄出来，而且我的9.2 确实没有 samplingunit 功能。。。。可能还没到9.22。。。所以还得再次请教各位高手了！！
不胜感激啊！

shiyiming

to caicaisas
very simple
add one more loop for BY variable

shiyiming

to oloolo
但是最后merge的时候出现问题了，假设sampnum是指示sample的变量， by sampnum cluster,会出问题，因为temp这个原数据里面没有sampnum这个变量，所以merge不上。。。。
唉，真的比较愚笨 <!-- s:? --><img src="{SMILIES_PATH}/icon_confused.gif" alt=":?" title="Confused" /><!-- s:? --> ，，，能不能稍微展示一下code?  多谢了！

shiyiming

原始程序里利用index key来寻到cluster的方法的确很好，为什么不采用呢？