发帖
每发1个新帖
可以获得0.5个丁当奖励
回帖

【原创】临床试验计算样本量的SAS程序---以文献为例

楼主

入门站友

1楼

这个帖子发布于5年零314天前，其中的信息可能已发生改变或有所发展。

如题，在研究《样本量估计及其在nQuery和SAS软件上的实现4——率的比较(三)》这篇文章时，对例子2-9进行的SAS编程实践，最终得出和文献一样的结果，自己写个程序不容易，以前有人发过相关SAS程序，但是根本无法运行，还是自己动手比较好，现在也把正确的SAS程序文档上传附件，供大家使用~~~

SAS9.2软件程序如下（都是我手动输入的，只需替换不同的a,s,p1,p2,power值即可）：
%macro PTT0(a,s,p1,p2,power);
data PTT0;
p1=&p1;
p2=&p2;
s=&s;
a=&a;
power=&power;
q1=1-p1;
q2=1-p2;
p=1/2*(p1+p2);
q=1-p;
za=probit(1-a/s);
zb=probit(power);
n=(za*sqrt(2*p*q)+zb*sqrt(p1*q1+p2*q2))**2/(p1-p2)**2;
n=ceil(n);
run;
proc printdata=PTT0 label;
var a s p1 p2power n;
label a='Test significance level'
s='1 or 2 sided test'
p1='proportion1'
p2='proportion2'
power='power'
n='n'n;quit;
%mend PTT0;
%PTT0(a=0.05,s=2,p1=0.45,p2=0.65,power=0.95);
输出的结果是：

临床试验计算样本量的SAS程序---以文献为例.docx(191.04k)

邀请讨论

搜索

2014-03-19 03:23 浏览 : 6444 回复 : 7

投票 2
收藏 11
打赏
引用
分享
- 微信扫一扫
- 新浪微博
- 丁香客
- 复制网址
举报
- 广告宣传推广
- 政治敏感、违法虚假信息
- 恶意灌水、重复发帖
- 违规侵权、站友争执
- 附件异常、链接失效
- 其他

• 2020，庚子，肾脏泌尿版给大家拜年了！

amtf_2018

铁杆站友

2楼

谢谢分享。
我把上面的程序拷到SAS中运行，没有出现运行结果，日志中显示
181: 无法确定 LINE 和 COLUMN。
NOTE: NOSPOOL 已打开。用 OPTION SPOOL 重新运行可恢复出错的 LINE 和 COLUMN。
ERROR 181-322: 过程名拼错。
NOTE: SAS 系统由于错误而停止了该步的处理。
NOTE: PROCEDURE PRINTDATA 所用时间（总处理时间）:
实际时间 0.15 秒
CPU 时间 0.03 秒

2014-03-19 08:28

投票
收藏
打赏
引用
分享
- 微信扫一扫
- 新浪微博
- 丁香客
- 复制网址
举报
- 广告宣传推广
- 政治敏感、违法虚假信息
- 恶意灌水、重复发帖
- 违规侵权、站友争执
- 附件异常、链接失效
- 其他

• 痛心！湖北一医生感染新型肺炎，刚刚不幸离世

楼主

雪雪的静

病案科

入门站友

3楼

gwzh008
谢谢分享。
我把上面的程序拷到SAS中运行，没有出现运行结果，日志中显示
181: 无法确定 LINE 和 COLUMN。
NOTE: NOSPOOL 已打开。用 OPTION SPOOL 重新运行可恢复出错的 LINE 和 COLUMN。
ERROR 181-322: 过程名拼错。
NOTE: SAS 系统由于错误而停止了该步的处理。
NOTE: PROCEDURE PRINTDATA 所用时间（总处理时间）:
实际时间 0.15 秒
CPU 时间 0.03 秒

我把帖子里的命令贴进去也是这样的结果，但是用附件word里的命令就好了，可能是排版的问题，会有一些偏差？

2014-03-20 09:04