关于本站
人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!
获取电子版《CDA一级教材》
完整电子版已上线CDA网校,累计已有10万+在读~ 教材严格按考试大纲编写,适合CDA考生备考,也适合业务及数据分析岗位的从业者提升自我。
TOP热门关键词
免费学术公开课,扫码加入![]() |
老师们帮忙看看该sas 宏错在哪里,运行总是出错!
谢谢!
%macro RR(delta,n1,n2,distr1,distr2,essai,report=1);
%global pow;
data -null-;
call symput("seed ",date());
run;
Data depart;
do numessai=1 to &essai;
do pat=1 to %eval(&n1+&n2);
if pat<=&n1 then do;
Distr=rantbl(&seed,&distr1);
group=1;
end ;
if pat> %eval(&n1) and pat <= %eval (&n1+&n2)
then do;
Distr=rantbl(&seed,&distr2);
group=2;
end;
output;
end;
end;
Run;
proc freq data=depart(where=(group in (1,2))) noprint;
by numessai;
tables group * distr /cmh2 scores=ridit;
output out=TRT-AB N LGRRC1;
run;
Data power1;
set TRT-AB;
if l-lgrrc1^=.;
power=l-lgrrc1>=δ
scenario="TRT-A/TRT-B";
Run;
proc univariate data=power1 noprint;
by scenario;
var power;
output out =power2 n=n sum=rep;
run;
Data power3;
length distr1 distr2 $ 30 scenario $ 20;
merge power2;
by scenario;
lcl=100-100*Betainv(0.975,n-rep+1,rep)+0;
ucl=100*Betainv(0.975,rep+1,n-rep)+0;
lC="["|| compress (put(lcl,best4.))||"%;"|| compress (put(ucl,best4.))||"%"||"]";
distr1=compress(symget("distr1")*100||"%");
distr2=compress(symget("distr2")*100||"%");
n1=symget("n1");
n2=symget("n2");
delta=symget("delta");
simul=symget("essai");
power=100*rep/n;
call symput('pow',put(power,best4.1));
Run;
%if %eval(&report)=1 %then %do;
Title2
Proc report data= power3 headline headskip nowindows split ="@" ls=120;
column("--" scenario simul delta power lc
("Expected sample size" "--" n1 n2)
("Expected rates" "--" distr1 distr2));
define scenario /order format=$18. left width=18 "Scenario";
define simul/display format=$5.center width=12 "N? of @ simulations";
define delta /display format=$5.center width=6"Delta";
define power/display format=4.1 center width=6"Power";
define lC/display format=$12.center width=12"95% cl of @ Power";
define n1/display format=$10.center width=12 "TRT-A";
define n2/display format=$10.center width=12 "TRT-B";
define distr1/display format=$12.center width=12 "TRT-A";
define distr2/display format=$12.center width=12 "TRT-B";
compute after;
line @31 108*"-";
endcomp;
fotnote "seed=&seed";
%end;
Run;
%mend RR;
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
您可能感兴趣的文章
- SAS软件培训 ... | 如何记录SAS逻辑库数据文件的详细 ...
- SAS软件培训 ... | sas9.2sid 2013年到期32位
- SAS软件培训 ... | [讨论]有关SAS90功能模块的问题
- SAS软件培训 ... | [求助]SAS9.13绿色版启动错误
- SAS软件培训 ... | 6周复习成功通过 SAS base and a ...
- SAS软件培训 ... | sas关于reg的问题
- SAS软件培训 ... | [分享]Using SAS in Financial R ...
- SAS软件培训 ... | sas marco programming made eas ...
人气文章
本文标题:高手帮忙看看该sas宏哪儿又错误
本文链接网址:https://bbs.pinggu.org/jg/ruanjianpeixun_sasruanjianpeixun_295967_1.html
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。



