搜索
人大经济论坛 附件下载

附件下载

所在主题:
文件名:  模拟数据.xls
资料下载链接地址: https://bbs.pinggu.org/a-1539753.html
附件大小:
最近需要用到sas来对我的一个数据进行1:1抽样,我在网上找了些资料看,没怎么看懂,于是自己想了一个思路来编写这个抽样程序,现分享给大家。
首先,数据库js如下,含有唯一变量ID,需要匹配的变量 shcool、grade、class和age,case是定义变量,1是病例,0是对照。抽样要求school、grade、class要相同,age相差不超过2岁。

先说我的思路,首先是生成病例库和对照库, 然后依次从病例库中抽选1条,然后按照条件从对照库选择符合的样本组成待抽选库,然后从这个待抽选库中随机抽选1条,最后把抽出的这条从对照库中剔除,然后重复上述过程直到每条病例都匹配完成,最后合并抽出来的样本。
代码如下:
  1. proc sort data=js out=a1;
  2. by school grade class case age;
  3. data a1;set a1;
  4. pp=compress(school||grade||class); *班级匹配变量pp;
  5. run;
  6. data a_case a_control;set a1;
  7. if case=1 then output a_case; *病例数据集a_case;
  8. if case=0 then output a_control; *对照数据集a_control;
  9. run;

  10. %macro ss;
  11. proc datasets lib=work;delete sample;run;*清除前次抽选样本,方便重复运行宏抽样程序;
  12. proc sql noprint;
  13. select count(*) into: num from a_case; *统计病例个数;
  14. %do i=1 %to #
  15. %let k=%eval(&i-1);
  16. proc sql noprint;
  17. select id into: idx separated by ' ' from a_case;*病例ID号 到宏idx;
  18. %let xx=%scan(&idx,&i,' ') ;
  19. data case&i;set a_case;
  20. if id=&xx then call symput("pp_v",pp); *相同匹配条件pp;
  21. if id=&xx then call symput("age_v",age); run; *范围匹配条件age;
  22. data control&i;set a_control;
  23. if pp="&pp_v" and abs(age-&age_v)<3; run; *生成符合匹配条件control样本;
  24. proc surveyselect noprint data=control&i method=srs n=1 out=sample&i seed=1000; *n=1 ,1:1匹配;
  25. data a_control;set a_control sample&i; run;
  26. proc sql noprint;
  27. create table a_control asselect * from a_controlgroup by id having n(id)=1; *剔除已选择control样本;
  28. %end;
  29. data sample;set sample:;run; *生成匹配样本库;
  30. data hb;set a_case sample;run; *生成总库;
  31. proc datasets lib=work;save a1 js a_case a_control sample hb; run; *清除过程文件;
  32. %mend;
  33. %ss;
复制代码
最后,说说程序缺陷,首先是运行效率低,因为是循环语句反复生成数据库,病例较多时,运行肯定慢; 其次未考虑到匹配未成功的情况;第三,程序还是有点啰嗦。 仅供参考。




    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。
(如有侵权,欢迎举报)
二维码

扫码加我 拉你入群

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

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

GMT+8, 2025-12-27 08:53