关于本站
人大经济论坛-经管之家：分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!
经管之家新媒体交易平台
提供"微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯"等虚拟账号交易，真正实现买卖双方的共赢。【请点击这里访问】
TOP热门关键词
非常有用的sas程序，拿来即可用，从基本的数据整理到高级数据分析，应有尽有。样例：/*-------------------------------------------------------------------*//*MultipleComparisonsandMultipleTestsUsingtheSAS(r ...
免费学术公开课,扫码加入 |
样例：
/*-------------------------------------------------------------------*/
/* Multiple Comparisons and Multiple Tests Using the SAS(r) System */
/* */
/* Program 2.4: Bonferroni and Sidak Adjusted p-values Using the DATA Step */
data one;
input test pval @@;
bon_adjp = min(1,10*pval);
sid_adjp = 1 - (1-pval)**10;
datalines;
1 0.0911 2 0.8912
3 0.0001 4 0.5718
5 0.0132 6 0.9011
7 0.2012 8 0.0289
9 0.0498 10 0.0058
;
proc sort data=one out=one;
by pval;
proc print data=one;
run;
/* Program 2.5: Bonferroni and Sidak Adjusted p-values Using PROC MULTTEST */
data one; set one;
rename pval=raw_p;
drop bon_adjp sid_adjp;
proc multtest pdata=one bon sid out=outp;
proc sort data=outp out=outp;
by raw_p;
proc print data=outp;
run;
/* Program 2.6: Conservative Simultaneous Confidence Intervals with Multivariate Data */
data _null_;
call symput('bonalpha',0.05/9 );
call symput('sidalpha',1-(1-0.05)**(1/9));
data HusbWive;
input HusbQ1-HusbQ4 WifeQ1-WifeQ4 @@;
DiffQ1 = HusbQ1-WifeQ1;
DiffQ2 = HusbQ2-WifeQ2;
DiffQ3 = HusbQ3-WifeQ3;
DiffQ4 = HusbQ4-WifeQ4;
DiffQAvg = sum(of HusbQ1-HusbQ4)/4 - sum(of WifeQ1-WifeQ4)/4;
DiffComp = sum(of HusbQ1-HusbQ2)/2 - sum(of WifeQ1-WifeQ2)/2;
DiffPass = sum(of HusbQ3-HusbQ4)/2 - sum(of WifeQ3-WifeQ4)/2;
DiffFFP= sum(of HusbQ1 HusbQ3)/2 - sum(of WifeQ1 WifeQ3)/2;
DiffFFY= sum(of HusbQ2 HusbQ4)/2 - sum(of WifeQ2 WifeQ4)/2;
datalines;
2 3 5 5 4 4 5 5 5 5 4 4 4 5 5 5 4 5 5 5 4 4 5 5
4 3 4 4 4 5 5 5 3 3 5 5 4 4 5 5 3 3 4 5 3 3 4 4
3 4 4 4 4 3 5 4 4 4 5 5 3 4 5 5 4 5 5 5 4 4 5 4
4 4 3 3 3 4 4 4 4 4 5 5 4 5 5 5 5 5 4 4 5 5 5 5
4 4 4 4 4 4 5 5 4 3 5 5 4 4 4 4 4 4 5 5 4 4 5 5
3 3 4 5 3 4 4 4 4 5 4 4 5 5 5 5 5 5 5 5 4 5 4 4
5 5 4 4 3 4 4 4 4 4 4 4 5 3 4 4 4 4 4 4 5 3 4 4
4 4 4 4 4 5 4 4 3 4 5 5 2 5 5 5 5 3 5 5 3 4 5 5
5 5 3 3 4 3 5 5 3 3 4 4 4 4 4 4 4 4 4 4 4 4 5 5
3 3 5 5 3 4 4 4 4 4 3 3 4 4 5 4 4 4 5 5 4 4 5 5
;
proc glm;
model HusbQ1-HusbQ4 WifeQ1-WifeQ4 = / nouni;
repeated Spouse 2, Question 4 identity;
run;
proc means alpha=0.05 n mean lclm uclm;
title "Unadjusted Confidence Intervals";
var DiffQ1-DiffQ4 DiffQAvg DiffComp DiffPass DiffFFP DiffFFY;
proc means alpha=&sidalpha n mean lclm uclm;
title "Simultaneous Sidak Intervals";
var DiffQ1-DiffQ4 DiffQAvg DiffComp DiffPass DiffFFP DiffFFY;
proc means alpha=&bonalpha n mean lclm uclm;
title "Simultaneous Bonferroni Intervals";
var DiffQ1-DiffQ4 DiffQAvg DiffComp DiffPass DiffFFP DiffFFY;
run;
/* Program 2.7: Multiple Tests with Multivariate Data */
proc means data=HusbWive n mean std prt;
title "Tests of Hypotheses With Husband/Wife Data";
var DiffQ1-DiffQ4 DiffQAvg DiffComp DiffPass DiffFFP DiffFFY;
run;
/* Program 2.8: HOLM Adjusted p-values Using PROC MULTTEST */
data one;
set one;
rename pval=raw_p;
drop bon_adjp sid_adjp;
proc multtest pdata=one bon stepbon out=outp;
proc sort data=outp out=outp;
by raw_p;
proc print data=outp;
run;
/* Program 2.9. Sidak-Holm Adjusted p-values Using PROC MULTTEST */
data one; set one;
rename pval=raw_p;
drop bon_adjp sid_adjp;
proc multtest pdata=one sid stepsid out=outp;
proc sort data=outp out=outp;
by raw_p;
proc print data=outp;
run;
「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
您可能感兴趣的文章
人气文章
1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责；
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性，不作出任何保证或承若；
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性，不作出任何保证或承若；
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
