SAS编程:频数汇总时如何处理分析分组种类不全的情况?

在项目Table的输出中,少不了对分类变量进行频数汇总,这一整个计算输出过程并不复杂。一般研究要求,分类变量的所有类别都需要输出,这里有一点需要注意,分析时很有可能出现个别分组分类计数为0的情况。

为了保证分析分组变量的种类完整,计数为0的类别的结果也需要输出。并且,编程中一些同事会机械地将分类计数结果赋值到宏变量中,如果不将计数为0的类别输出,后续调用宏变量会出现Error。下面介绍两个具体的处理方法。

Dummy数据集

处理分类计数为0的情况,在日常工作中想必大家非常熟悉。常见的做法是生成一个包含所有种类的Dummy数据集,与有缺失种类的分析结果数据集进行拼接,补全计数为0的那一个类别的结果。这里用SASHELP.Class数据集简单演示一下:

**Get data from sashelp.class;
data class;
  set sashelp.class;
  where sex = "M";
run;

**Count for each sex category;
proc freq data = class;
  tables sex / out = result1(keep = sex count);
run;

**Create a dummy dataset for missing category;
data dummy;
  length sex $1;
  sex = "F"; output;
  sex = "M"; output;
run;

**Merge result dataset with dummy dataset;
data result;
  merge dummy result1;
  by sex;

  if count = . then count = 0;
run;

Result

结果数据集中输出了Sex=F这一类别计数为0的情况了。

Preloadfmt选项

Preloadfmt是Means过程步中Class语句的选项之一,SAS官方文档中是这样描述Preloadfmt选项的——指定分类变量的所有格式都预先加载(specifies that all formats are preloaded for the class variables.,来源:SAS Help Center: Syntax: PROC MEANS CLASS Statement)。

使用Preloadfmt选项时,分组分析的所有种类的Format都会输出,包括计数为0的种类。只要提前设置好分组变量的Format,就可以输出计数为0的那一组类别。官方文档中也说明,如果想要输出计数为0的种类,PROC语句中要使用COMPLETETYPES选项。

需要注意的是,Means过程步的分析变量只能是数值变量,所以计数分类变量数目时需要一个数值的标记变量,通常我是新建一个变量flag = 1;

具体的演示代码如下:

**Get data from sashelp.class;
data class2;
  set sashelp.class;
  where sex = "M";

  *Create count flag;
  flag = 1;
run;

**Set format for class var;
proc format;
  value $sex
    "F"="Female"
    "M"="Male"
  ;
run;

**Get count;
proc means data = class2 completetypes;
  class sex/ preloadfmt;
  var flag;
  output n=count out = result2;

  format sex $sex.;
run;

输出结果如下,我们可以看到Sex=F的那一组类别计数为0的记录也输出出来。

Result2

但是,这个结果也多输出了一行两组汇总的数据(_type_=0),在调用Proc语句时,使用nway选项就可以移除这一行记录。当然,这一行也可以在输出数据集的数据集选项中筛选删除(out = result2 (where=(_type_ ne 0)))。但是,遇到多个分组变量时,type变量的取值就不只是0和1这两个了,进行手动删除还需要查看数据集中的具体输出值,这里其实已经涉及到Hardcoding了。所以,建议直接使用nway选项,移除多余的汇总计数记录。

**Get count;
proc means data = class2 nway completetypes;
  class sex/ preloadfmt;
  var flag;
  output n=count out = result2;

  format sex $sex.;
run;
Result3

结语

计数处理在临床试验SAS编程中非常常见,Preloadfmt选项可以简化计数编程过程,这个简化效果在分组变量数目越多、分组变量种类越多时,越明显。关于这个选项的延伸应用,会在后续文章中介绍。

感谢阅读!若有疑问,欢迎评论区交流!

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 205,132评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 87,802评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,566评论 0 338
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,858评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,867评论 5 368
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,695评论 1 282
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,064评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,705评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 42,915评论 1 300
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,677评论 2 323
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,796评论 1 333
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,432评论 4 322
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,041评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,992评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,223评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,185评论 2 352
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,535评论 2 343

推荐阅读更多精彩内容