论坛币求大神编个sas程序

如上图，要对480个人进行实验对照分组，按顺序4人一组，实验组：对照组=3：1。这四个人的处理方式随机分配，比如1到4号分组实验组，对照组，对照组，对照组，对照组；5-8号分组 对照组，对照组，对照组，实验组 。
具体要求如下

假设stratum为分层因素，分4层。

proc plan seed=20151223;
      factors stratum=4 block=30 rand_code=4/noprint;
      output out=a;
run;
data b;
   set a;
   id=_n_;
   if rand_code in(1 2 4) then do;ARM=0;TRT="试验组";end;
   else do;ARM=1;TRT="对照组";end;
run;

3点提醒
1、试验组：对照组=3 : 1
2、没有涉及到分层
3、'试验组'而非'实验组'

proc plan seed=20151223;
      factors block=120 rand_code=4/noprint;
      output out=a;
run;
data b;
   set a;
   id=_n_;
   if rand_code in(1 2 4) then do;ARM=0;TRT="试验组";end;
   else do;ARM=1;TRT="对照组";end;
run;

孤单的我们 发表于 2015-12-23 16:48
3点提醒
1、试验组：对照组=3 : 1
2、没有涉及到分层

额 不显示结果。。。
那要是分层区组随机的怎么编写呢？谢谢

孤单的我们 发表于 2015-12-23 16:48
3点提醒
1、试验组：对照组=3 : 1
2、没有涉及到分层

额 不显示结果。。。
那要是分层区组随机的怎么编写呢？谢谢

孤单的我们 发表于 2015-12-23 16:07
假设stratum为分层因素，分4层。

proc plan seed=20151223;

分数给你了。采纳了。谢谢了
分四层？那就是每层三十人，那就不行了，实验有五个中心，每个中心的人数是不一样的。估计只能按照你给我的第一种来了。这样真的不算事分层了

terrorist002 发表于 2015-12-23 19:54
分数给你了。采纳了。谢谢了
分四层？那就是每层三十人，那就不行了，实验有五个中心，每个中心的人数是 ...

这方面虽然我并没有什么经验，不过可以提个思路。

5层不均的话，可以选取最大的一层人数N,把N适当的增加几个数，变成M（4的倍数）。

factors stratum=5 block=M/4 rand_code=4

每个中心按照入组顺序给随机号，多的随机号就用不到了。

这个办法是否在实际编盲中可行，有待大神指教。

孤单的我们 发表于 2015-12-25 09:20
这方面虽然我并没有什么经验，不过可以提个思路。

5层不均的话，可以选取最大的一层人数N,把N适当的增 ...

https://bbs.pinggu.org/thread-4150796-1-1.html
我有发了个帖子，写的详细些了。这个要求做出来好像真的不算分层了。但是要求这么干，只好。。。