如何利用SAS进行统计分析并自动将所需结果输出到数据集中(以T检验为例)

/**************************************************************** * 功能: 两样本T检验
* 具体用途：本文档所提供的程序可用于进行T检验（T test）时将需要的结果自动一个* 数据集，包括平均值(mean)，标准差(sd)，样本量(sample size)，
* T统计量，自由度以及P值(该P值是程序根据方差齐性检验结果自动选出来* 的。

)
* 版权所有:
* 日期: 20120314
* 版本号: 1.0
****************************************************************/
/************************************************
* Function: Mock two sample T-test
* Copy Right:
* Date: 20120314
* Edition: 1.0
*************************************************/
/* Generate mock data */
data Math;
do i=1to30;
Score=60.0+sqrt(81)* RANnor(1);
Sex="Female";
output;
end;
do i=31to60;
Score=80.0+sqrt(10)* RANnor(1);
Sex="Male";
output;
end;
run;
/*Find the names of output objects using the ODS Trace statement,
this step can be omitted if you are familar with the name of output objects of T test*/
ods trace on;
proc ttest data=Math;
var SCORE;
class sex;
run;
ods trace off;
/*Get Mean and SD */
ods output Statistics=Mean1;
proc ttest data=Math ;
var score;
class sex;
run;
ods output close;
/*Get results of Variance Equality test(Folded)*/
ods output Equality=EP;
proc ttest data=Math ;
var score;
class sex;
run;
ods output close;
/*Get P values of T test */
ods output TTests=TP;
proc ttest data=Math ;
var score;
class sex;
run;
ods output close;
/*Keep Variable Class N Mean StdDev */
DATA Mean(KEEP=Subject Variable Class N Mean StdDev );
set Mean;
where N GT 0;
Subject="Math";
RUN;
/*Change the format of P vales*/
DATA TP;
format ProbT 10.8;
set TP;
Subject="Math";
RUN;
DATA EP(KEEP=Subject Variable ProbF );
format ProbF 10.8;
set EP;
Subject="Math";
RUN;
/*Choose P value according to the Variances equality test */ PROC SORT DATA=EP;
BY Subject;
run;
PROC SORT DATA=TP;
BY Subject;
run;
DATA P;
merge ep tp;
by Subject;
run;
data p;
set p;
if ProbF <0.05 and Variances="Equal"then delete;
else if ProbF >0.05 and Variances="Unequal"then delete; where ProbT>0;
run;
/*Combine Mean(SD) and P Value*/
PROC SORT DATA=P;
BY Subject;
run;
PROC SORT DATA=Mean;
BY Subject;
run;
Data ALL(keep=Subject Variable Class N Mean StdDev tValue DF ProbT );
retain Subject Variable Class N Mean StdDev tValue DF ProbT; MERGE MEAN P;
BY Subject;
RUN;。