临床试验非劣效试验的样本量模拟计算SAS宏

Rock2000

10046

收藏 2013-09-14

悬赏 20 个论坛币已解决

这是张高魁, 夏结来, 姚晨《相对率的非劣效性试验检验效能及样本量的模拟计算方法及SAS实现》的程序
例: 一种足癣标准治疗方案, 疗程结束时的真菌学清除率为80%。现提出一种相对简单的治疗方案, 为考察对比新方案与标准治疗方案的疗效, 拟设计一个非劣效性试验, 临床提出新方案的疗效达到标准方案清除率的90% 作为临床非劣效的相对标准, 即: 当新方案与标准方案在疗程末的真菌学相对清除率( 新方案P标准方案) 的95% 可信区间下限不低于0. 9时, 即可判断其临床非劣效于标准方案。问题: ( 1) 、病例为100 对, 检验的效能是多少? ( 2) 、要使效能达到0. 8, 需要多少病例?
问题（1）解决的SAS宏如下，运行没有问题。

复制代码

最佳答案

jingju11 查看完整内容

Appreciate it. 我看了原文,方法比较容易理解。遗憾的是没有找到其中几个公式的出处。原来的程序，可能因为中英文乱码的原因，很难一一纠查。索性自己编写。根据作者的思路，我差不多照搬了第一个程序%RR。不同的是我直接用BINOMIAL数据，尤其对于尺寸大的数据，效率应该有所改善。并且抑制输出。原来的第二个程序%PSIZE看起来比较复杂些，没有细读。条条大路通罗马。我使用二分算法，效率虽然不高，但是简单。 http://blog.sin ...

扫码加我拉你入群

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

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

全部回复

jingju11

2013-9-14 13:35:01

Appreciate it.
我看了原文,方法比较容易理解。遗憾的是没有找到其中几个公式的出处。原来的程序，可能因为中英文乱码的原因，很难一一纠查。索性自己编写。根据作者的思路，我差不多照搬了第一个程序%RR。不同的是我直接用BINOMIAL数据，尤其对于尺寸大的数据，效率应该有所改善。并且抑制输出。原来的第二个程序%PSIZE看起来比较复杂些，没有细读。条条大路通罗马。我使用二分算法，效率虽然不高，但是简单。
http://blog.sina.com.cn/s/blog_a3a926360101fp1q.html
京剧

扫码加我拉你入群

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

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

Rock2000

2013-9-14 13:41:43

现在问题是问题（2）的SAS宏，我安照文章编写，最后运行出错，程序如下：

复制代码

出错提示如下：

复制代码

哪哪位帮忙修正下？谢谢！！

扫码加我拉你入群

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

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

jingju11

2013-9-15 00:00:39

我觉得换个思路：
解决作者提出的问题，如果用sas提供的过程power，差不多就5-10行程序，而且提供的信息更全面。
针对以上的宏：在参数未知的情况下，一味使用%EVAL（）来估计和判断，不出错误只是运气。
京剧

扫码加我拉你入群

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

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

webgu

2013-9-15 09:10:15

jingju11 发表于 2013-9-15 00:00
我觉得换个思路：
解决作者提出的问题，如果用sas提供的过程power，差不多就5-10行程序，而且提供的信息更 ...

PROC POWER? 愿闻其详？

扫码加我拉你入群

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

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

yongyitian

2013-9-15 11:41:11

把一楼，31行改为

power=L_lgrrc1>=δ

后一楼和二楼的宏都可以运行. 只是二楼的宏要手动停止，还有一楼 25行处有一个 warning:
NOTE: Line generated by the invoked macro "RR".
5    distr /cmh2 scores=ridit; output out=TRT_AB N L_LGRRC1; run;
                                                      --------
                                                      1

WARNING 1-322: Assuming the symbol LGRRC1 was misspelled as L_LGRRC1.

NOTE: There were 410160 observations read from the data set WORK.DEPART.
   WHERE group in (1, 2);

扫码加我拉你入群

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

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

点击查看更多内容…

Rock2000

2013-9-15 17:42:49

yongyitian 发表于 2013-9-15 11:41
把一楼，31行改为

power=L_lgrrc1>=δ

谢谢回复。

其实是改成output out=TRT_AB N lgrrc1;
(lgrrc1=Logit common relative risk, column 1 .)

正确的宏程序如下：

复制代码

扫码加我拉你入群

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

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

Rock2000

2013-9-15 18:07:24

其实，在实际工作中，此例子是是比较少见的，此例假设如下：

复制代码

问题出现在“清除率的95%可信区间下限0.9”这里，在实际的临床试验中，不会知道这个下限的，仅知道非劣效界值，此值通常在5%~15%之间，不会超过15%，当然，界值约小需要样本量越多，比如某临床试验定了非劣效界值是15%，试验组及对照组的估计率是80%，这里的率95%可信区间下限不能等于80%-15%=65%。

如果本例中的“清除率的95%可信区间下限0.9”改为非劣效界值是15%，疗程结束时的真菌学清除率为80%，通过PASS计算得出每组样本量是112例。

当然，也可以通过公式计算，公式如下。

附件列表

GS.png

原图尺寸 7.1 KB

扫码加我拉你入群

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

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

webgu

2013-9-15 20:52:31

Rock2000 发表于 2013-9-15 18:07
其实，在实际工作中，此例子是是比较少见的，此例假设如下：

问题出现在“清除率的95%可信区间下限0.9” ...

楼上的大哥，PASS是买的？

扫码加我拉你入群

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

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

Rock2000

2013-9-15 22:41:48

webgu 发表于 2013-9-15 20:52
楼上的大哥，PASS是买的？

下载的 http://www.ncss.com/download/pass/updates/pass11/

扫码加我拉你入群

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

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

webgu

2013-9-16 08:51:39

Rock2000 发表于 2013-9-15 22:41
下载的 http://www.ncss.com/download/pass/updates/pass11/

谢谢，好像只能试用7天。

扫码加我拉你入群

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

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

Rock2000

2013-9-16 09:42:03

jingju11 发表于 2013-9-16 08:45
Appreciate it.
我看了原文,方法比较容易理解。遗憾的是没有找到其中几个公式的出处。原来的程序，可能因为 ...

拜读了你的程序，佩服。

正如我在上帖所说，本例子最大不足是假设知道率的95%下限来模拟某样本量情况下的Power值，但实际情况是，实际临床试验中提出的假设是知道非劣效界值，率的95%是不固定的，这样，模拟不同的样本量的情况下的Power用本例程序就不能计算。
实际上，在知道非劣效界值情况，模拟不同的样本量下的power值，用PASS模拟是很方便的，请看贴图。

扫码加我拉你入群

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

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

jingju11

2013-9-16 09:54:28

what is 率的95%? 你指的是TYPE I ERROR? JINGJU

扫码加我拉你入群

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

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

Rock2000

2013-9-16 10:04:11

webgu 发表于 2013-9-16 08:51
谢谢，好像只能试用7天。

这个你拿去。

k.rar
大小:(12.25 KB)

马上下载

本附件包括：

digerati.nfo
keygen.exe

扫码加我拉你入群

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

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

Rock2000

2013-9-16 10:07:00

jingju11 发表于 2013-9-16 09:54
what is 率的95%? 你指的是TYPE I ERROR? JINGJU

我没有说清楚，是非劣效率的95%可信区间（CI）下限

扫码加我拉你入群

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

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

jingju11

2013-9-16 21:13:06

我可能没有完全理解.
我想是否可以通过改变FREQ里的ALPHA=..京剧

扫码加我拉你入群

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

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

Rock2000

2013-9-16 21:43:55

jingju11 发表于 2013-9-16 21:13
我可能没有完全理解.
我想是否可以通过改变FREQ里的ALPHA=..京剧

不成的。

扫码加我拉你入群

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

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

sunny13gao

2013-9-25 16:07:36

thank you!

扫码加我拉你入群

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

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

jiweisi

2014-7-13 21:43:00

Rock2000 发表于 2013-9-14 13:41
现在问题是问题（2）的SAS宏，我安照文章编写，最后运行出错，程序如下：

出错提示如下：

请问一下这个最后怎么编程的，多谢指教

扫码加我拉你入群

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

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

最佳答案

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群