proc freq计算客观缓解率中的几个选项

在肿瘤项目的table中,往往会有这样几张表计算最佳总体应答、客观缓解率、缓解持续时间等,比如今天我了解的这个客观缓解率,shell上有说通过Clopper Pearson方法计算ORR的95%置信区间(CI)。

好吧,什么方法能计算置信区间并且采用Clopper Pearson方法,去搜,proc freq可以,,proc freq我们用的最多的应该是去计算频数频率吧。



那什么选项能计算95%CI呢?看SAS官网,这个binomial选项可以,好的,接下来我们建立程序的大致框架。

data temp;

do i = 1 to 20; group='A'; response=0; output; end;

do i = 1 to 80; group='A'; response=1; output; end;

do i = 1 to 100; group='B'; response=0; output; end;

do i = 1 to 100; group='C'; response=1; output; end;

run;

proc freq data=temp;

by group;

    tables response / binomial(all);

ods output BinomialCLs=CLS;

run;

选项解释:

这个ALL代表用各种方法计算95%CI,包括下面几种方法,也是输出(type变量)

如果只是用Clopper Pearson方法,可以这样写:

tables response / binomial(exact);

exact就代表Clopper Pearson方法,注意exact的位置,它是放在tables后面的,只是作为一个tables的附属选项。

但是exact在tables还是一个语句,具体用法自己去查。

但是这样算出的好像不是我们想要的结果?为什么?因为ORR计算的是PR和CR=“Y”的受试者(在这个例子中可以看作是response=1,但是上面程序那样写好像把response=0也计算在内了,也就是把PR和CR=“N”的受试者计算在内了<但是实际上不是的>)

刚开始我想只筛选PR和CR=“Y”的受试者不就好了?但是结果算出来的结果却跟别人不一样,一直想不明白是为什么。虽然后面知道只筛选PR和CR=“Y的受试者会对统计结果产生影响,可能在tables这个过程发生了某些变化吧,所以不能只筛选特定受试者。

但是我看别人的程序添加了这样一个选项:

tables response / binomial(exact level=“Y”);这不也是只计算PR和CR=“Y”的受试者吗?

问题就出在这了!这个level不是代表只计算PR和CR=“Y”的受试者,而是把值为“Y”的先计算?或者说是当做重中之重去计算。看红框中官网的描述。

或者是网友的描述:

The LEVEL= option refers to the variable level (or order), not to the actual value itself. 

如果不加level,SAS默认把第一个出现的level(我不知道是不是第一条观测的值,还不确定)当作重中之重去计算。你可能还是会觉得这不就跟只筛选PR和CR=“Y”的受试者是一样的,怎么说呢?

你只筛选PR和CR=“Y”的受试者,那么avalc的值只有“Y”,但是用level,可以说PR和CR=“Y”和PR和CR=“N”都参与了95%CI的计算(不知道tables原理是什么样的),只是权重不一样。(上面是我自己的理解),

我们把上面数据集的结果加上level,并且只计算group=A看看结果

数据集先出现的response值是0,那么我们level 0看看是不是跟上面的数据集结果一样

发现结果确实一样!分别为0.1267和0.2918

然后我们level 1看看,发现结果变了。验证了上面我所讲的(至于不加level选项,到底是不是取观测第一条的值作为默认level我还是不确定)

所以总结一下,就是大家在求95%CI的时候不要提前筛选PR和CR=“Y”的受试者,虽然你知道你只需要PR和CR=“Y”的比例,但是会影响到结果。

还有这个ods output BinomialCLs=CLS;前面一个是统计量,后面是数据集名字(用这个选项你得提前在tables后面指定用什么方法,要不然像下面这样会报错)

proc freq过程创建了很多统计量,可以通过ods output输出,下面还有很多,几十个,这么多统计量,有时候真不知道输出哪些。

还有SAS官网说这个binomial只适合one-way table,但是上面3个group为什么都输出了?这个one-way table是怎么理解的?2X2列联表好理解,但是如果像上面的group=A,它的response有两个值=0,=1,它是one-way table,这个有时间也要去了解。

真的,到后面如果你主动去了解新知识,会越来越发现你怎么还有这么多东西不知道?,怎么还有这么多东西要学!!!然后什么都想学,一下想学这个,一下想学那个,结果都学不精,然后陷入深深的焦虑中。所以有时候我告诉自己,不要急,不要慌,慢慢来,学一点你就多知道一点,至少你肯主动去学去了解,这句话也送给你们。

封面的这两只狗真的太可爱了,我尤其喜欢左边这一只!嗯嘛

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

推荐阅读更多精彩内容

  • 《英国医学会杂志》(BMJ)自2008年9月开始至2015年由两位流行病与统计学专家不间断地出了300多期stat...
    简浅zx阅读 547评论 0 0
  • 16宿命:用概率思维提高你的胜算 以前的我是风险厌恶者,不喜欢去冒险,但是人生放弃了冒险,也就放弃了无数的可能。 ...
    yichen大刀阅读 6,042评论 0 4
  • 公元:2019年11月28日19时42分农历:二零一九年 十一月 初三日 戌时干支:己亥乙亥己巳甲戌当月节气:立冬...
    石放阅读 6,877评论 0 2
  • 今天上午陪老妈看病,下午健身房跑步,晚上想想今天还没有断舍离,马上做,衣架和旁边的的布衣架,一看乱乱,又想想自己是...
    影子3623253阅读 2,909评论 1 8