全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
2186 2
2015-03-14

〔例2-17〕在一项某中药治疗感染性腹泻的III 临床试验中,设计安排3 种干预,即高剂量药物组、低剂量药物组和安慰剂组,主要疗效指标为愈显率。根据

II 期临床试验数据,3个组的愈显率分别为60%50%40%,试以平衡设计、双侧检验水准0. 05、检验效能90%估计最低样本量。

SAS9. 2 软件实现:

%macro PGT0( pvalphagpower) ;

data PGT0;

g = g;

p1 = 0. 6;

p2 = 0. 5;

p3 = 0. 4;

alpha = alpha;

power = power;

if p = . then p = ( p1+ p2 + p3) /3;

else p = p;

if v = . then v = ( (p1-p) **2 + ( p2-p) **2 + ( p3-p) **2) /g;

else v = v;

h = v /( p* ( 1-p) ) ;

n = 0;

do until( abs( p1-power) = 0. 001) ;

p1 = 1-cdf ( 'CHISQUARE'cinv (1-alphag-1) g-1g * n* h) ;

end;

n = n + 1;

p1 = 0. 6;

run;

proc print data = PGT0 label;

var alpha p1 p2 p3 g power n;

label alpha = 'Test significance level'

p1 = 'proportion1'

p2 = 'proportion2'

p3 = 'proportion3'

power = 'Power'

n = 'n';

quit;

%mend PGT0;

%PGT0( p = . v = . alpha = 0. 05g = 3power =0. 90) ;

二维码

扫码加我 拉你入群

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

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

全部回复
2015-3-15 23:33:25
Some typos in the code. Plus, the loop has a flaw: power may not change continuously, which may cause the loop conditions have never been met and thus cause an infinite loop.
Jingju
replace the old code:
复制代码
by the new code:
复制代码
二维码

扫码加我 拉你入群

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

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

2015-3-20 08:51:19
jingju11 发表于 2015-3-15 23:33
Some typos in the code. Plus, the loop has a flaw: power may not change continuously, which may caus ...
算出所需要的sample size, 下一个问题就是,如何验证计算的样本大小的准确性呢?通常讲,样本的计算依托于检验方法.这就是给出利用bootstrapping方法的可能性.
我在BLOG里给出方法.可以阅读 Power of R by 2 table by jingju11
京剧
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群