Question

在SAS中，我创建了一个程序，该程序将从数据集中随机抽取50个观测值，并计算观测值的平均值。

data subset (drop=i samplesize);
samplesize=50;
obsleft=totobs;
do i=1 to samplesize;
  obsnum=ceil(ranuni(0)*totobs);
  set sashelp.baseball point=obsnum nobs=totobs;
  output;
end;
stop;
run;


proc sql;
select mean(nHome) from subset;
quit;

我想编辑代码，因此它将生成10个独立的随机样本而不是一个（我知道Proc SurverySelect中的reps =，但我不应该在这里使用它）。谢谢

Answer 1

k / n算法选择固定数量的样本，每个样本的概率为1 / n。

%let SEED = 1234;

data mySurveySelection;

  retain k 10; drop k;
  length select_n_ 8;

  set sashelp.baseball nobs=n;

  if (ranuni(&SEED) <= k/n) then do;
    k = k - 1;
    select_n_ = _n_;
    output;
  end;
  n = n - 1;

  if n = 0 then stop;
run;

你没有要求证明这些选择确实是1 / n，所以我不会证明这一点。

SurveySelect通常用于任何生产级别的研究或代码库。

SAS随机抽样

1 个答案: