开源宏库计划_11：%stat_countbystrata

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

名称: %stat_countbystrata
简介: 为数据集中的变量计算缺失数和缺失率
依赖子宏: %var_exist
说明:
    1. 该宏采用的表查询方式速度较慢，将在未来的版本改进；
    2. 结果数据集中，例数(Number of Subjects)是当前分层组合下受试者(由SUBJECT=指定)去重后的计数，例次(Number of Events)是当前分层组合下受试者的直接计数；
    3. 结果数据集按照分层1的全部组别的例数降序排列，然后按分层1水平值升序，然后按分层2的全部组别的例数降序，然后按分层2水平值升序，依此类推；
    4. 在一个分层中，合计水平总是排在前面，未编码的水平总是排在后面，而其它水平排在中间；
    5. 该宏支持对多层数据的分析，原型是临床数据分析中“AE按照系统器官分类和首选术语的分层分析”。理论上，最多可以支持 1023 个层；
参考: SAS帮助文档、《A-P-P Advanced Data Management Functions》
源代码: stat_countbystrata.zip (2.7 KB) 本附件包括：

• stat_countbystrata.sas

2023-4-13 22:55:13 上传

%stat_countbystrata

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Count Management Managemen Functions function

用法示例：
*首先创建用于测试的数据集;

1. data adsl(drop=i);
   
2.   do i=1 to 300;
   
3.     usubjid='A167'||put(i,z3.);
   
4.     trt01an=rantbl(0,0.3333,0.3333,0.3334);
   
5.     output;
   
6.   end;
   
7. run;
   
8. 
   
9. data adae;
   
10.   length usubjid $20.  aeterm aebodsys aedecod aellt $200.;
    
11.   call streaminit(1);
    
12.   do i=1 to 300;
    
13.     p=rand('uniform',0,1);
    
14.     if p<=0.8 then do;
    
15.       usubjid='A167'||put(i,z3.);
    
16.       do j=1 to int(rand('uniform',0,7));
    
17.         aeterm='AE'||put(j,z1.);
    
18.         k=round(rand('uniform',1,5));
    
19.         a=round(rand('uniform',1,5));
    
20.         b=round(rand('uniform',1,5));
    
21.         aebodsys='SOC'||put(k*100,z3.);
    
22.         aedecod='PT'||put(k*100+a*10,z3.);
    
23.         aellt='LLT'||put(k*100+a*10+b,z3.);
    
24.         output;
    
25.       end;
    
26.     end;
    
27.   end;
    
28.   drop i p j k a b;
    
29. run;
    
30. 
    
31. data adae;
    
32.   set adae;
    
33.   do i=1 to rec until(trtan);
    
34.     set adsl(keep=usubjid trt01an rename=usubjid=tmpusubjid)nobs=rec point=i;
    
35.     if usubjid=tmpusubjid then trtan=trt01an;
    
36.   end;
    
37.   drop tmpusubjid i;
    
38. run;

复制代码

*调用宏%stat_countbystrata;

1. %stat_countbystrata(dsin=adae,dsout=myout,denom=adsl,strata=aebodsys aedecod aellt,class=trt01an,subject=usubjid);

复制代码

*结果截图;

谢谢分享！

谢谢分享！