搜索
人大经济论坛 › 附件下载

附件下载

所在主题：
文件名:  k.rar
资料下载链接地址: https://bbs.pinggu.org/a-1401850.html
本附件包括：
  • digerati.nfo
  • keygen.exe
附件大小:
12.25 KB   举报本内容
这是张高魁, 夏结来, 姚晨《相对率的非劣效性试验检验效能及样本量的模拟计算方法及SAS实现》的程序
例: 一种足癣标准治疗方案, 疗程结束时的真菌学清除率为80%。现提出一种相对简单的治疗方案, 为考察对比新方案与标准治疗方案的疗效, 拟设计一个非劣效性试验, 临床提出新方案的疗效达到标准方案清除率的90% 作为临床非劣效的相对标准, 即: 当新方案与标准方案在疗程末的真菌学相对清除率( 新方案P标准方案) 的95% 可信区间下限不低于0. 9时, 即可判断其临床非劣效于标准方案。问题: ( 1) 、病例为100 对, 检验的效能是多少? ( 2) 、要使效能达到0. 8, 需要多少病例?
问题（1）解决的SAS宏如下，运行没有问题。
  1. %macro RR(delta,n1,n2,distr1,distr2,essai,report=1);
  2. %global pow;
  3. data _null_;
  4. call symput("seed ",date());
  5. run;
  6. Data depart;
  7. do numessai=1 to &essai;
  8. do pat=1 to %eval(&n1+&n2);
  9. if pat<=&n1 then do;
  10. Distr=rantbl(&seed,&distr1);
  11. group=1;
  12. end ;
  13. if pat> %eval(&n1) and pat <= %eval(&n1+&n2) then do;
  14. Distr=rantbl(&seed,&distr2);
  15. group=2;
  16. end;
  17. output;
  18. end;
  19. end;
  20. Run;
  21. /**对每组随机数进行CMH检验，并计算相对率；*/
  22. proc freq data=depart(where=(group in (1,2))) noprint;
  23. by numessai;
  24. tables group * distr /cmh2 scores=ridit;
  25. output out=TRT_AB N l_LGRRC1;
  26. run;
  27. /** 计算达到非劣效标准的次数；*/
  28. Data power1;
  29. set TRT_AB;
  30. if l_lgrrc1^=.;
  31. power=l_lgrrc1>=δ
  32. scenario="TRT-A/TRT-B";
  33. Run;
  34. proc univariate data=power1 noprint;
  35. by scenario;
  36. var power;
  37. output out =power2 n=n sum=rep;
  38. run;
  39. /** 计算效能，并准备报告；*/
  40. Data power3;
  41. length distr1 distr2 $ 30 scenario $ 20;
  42. merge power2;
  43. by scenario;
  44. lcl=100-100*Betainv(0.975,n-rep+1,rep)+0;
  45. ucl=100*Betainv(0.975,rep+1,n-rep)+0;
  46. lC="["|| compress (put(lcl,best4.))||"%;"|| compress (put(ucl,best4.))||"%"||"]";
  47. distr1=compress(symget("distr1")*100||"%");
  48. distr2=compress(symget("distr2")*100||"%");
  49. n1=symget("n1");
  50. n2=symget("n2");
  51. delta=symget("delta");
  52. simul=symget("essai");
  53. power=100*rep/n;
  54. call symput('pow',put(power,best4.1));
  55. Run;
  56. /**当报告开关为'开'时，报告模拟结果；*/
  57. %if%eval(&report)=1 %then %do;
  58. Title2"以相对率为指标时检验效能的计算结果";
  59. Proc report data= power3 headline headskip nowindows split ="@" ls=120;
  60. column("--" scenario simul delta power lc("Expected sample size" "--" n1 n2)("Expected rates" "--" distr1 distr2));
  61. define scenario /order format=$18. left width=18 "Scenario";
  62. define simul/display format=$5.center width=12 "N of @ simulations";
  63. define delta /display format=$5.center width=6"Delta";
  64. define power/display format=4.1 center width=6"Power";
  65. define lC/display format=$12.center width=12"95% cl of @ Power";
  66. define n1/display format=$10.center width=12 "TRT-A";
  67. define n2/display format=$10.center width=12 "TRT-B";
  68. define distr1/display format=$12.center width=12 "TRT-A";
  69. define distr2/display format=$12.center width=12 "TRT-B";
  70. compute after;
  71. line @31 108*"-";
  72. endcomp;
  73. footnote "seed=&seed";
  74. %end;
  75. Run;

  77. %mend RR;

  79. %RR(0.9, 100, 100, 0.80, 0.80, 1000);
复制代码


    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载，请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品，拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知，将积极的采取必要措施；同时，本站也将在技术手段和能力范围内，履行版权保护的注意义务。
(如有侵权，欢迎举报)
二维码

扫码加我 拉你入群

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

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

京ICP备16021002号-2 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问：王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2025-12-26 08:50