sas统计分析代码 - 图文 

/*统计描述*/ /*制作频数表*/ data a1; input x@@; low= ; dis= ;

z=x-mod(x-low,dis); cards; 1 2 3 4 ;

proc freq; tables z; run;

/*用meaNS过程计算例数、均数、标准差、标准误和95%可信区间*/ proc means n mean std stderr clm; var x; run;

/*用means过程对频数表计算例数、均数、标准差、最大值最小值*/ proc means; freq f; var x; run;

proc univariate; freq f; var x; run;

/*正态性检验*/

proc univariate normal plot; var x; run;

/*t检验*/ data aa3; input x;

y=x-比较的数；

/*配对资料d=x1-x2*/ cards; ;

proc univariate normal; var y; run;

/*用MEANS作配对资料两个样本均数比较的t检验*/ data ex3_6;

input x1 x2 @@; d=x1-x2; cards;

0.840 0.580 0.591 0.509 0.674 0.500 0.632 0.316 0.687 0.337 0.978 0.517 0.750 0.454

0.730 0.512 1.200 0.997 0.870 0.506 ;

proc means t prt; var d; run;

/*用UNIVARIATE过程作配对资料两样本均数比较的t检验*/ proc univariate data=ex3_6; var d; run;

proc ttest data=ex3_6; var d; run;

data ex3_7; input x @@;

if _n_<21 then c=1; else c=2; cards;

-0.70 -5.60 2.00 2.80 0.70 3.50 4.00 5.80 7.10 -0.50 2.50 -1.60 1.70 3.00 0.40 4.50 4.60 2.50 6.00 -1.40 3.70 6.50 5.00 5.20 0.80 0.20 0.60 3.40 6.60 -1.10 6.00 3.80 2.00 1.60 2.00 2.20 1.20 3.10 1.70 -2.00 ;

proc ttest; var x; class c; run;

/*方差分析*/ DATA RANDOM;

INPUT X GROUP @@; CARDS; 3.53 1 4.59 1 4.34 1 2.66 1 3.59 1 3.13 1 2.64 1 2.56 1 3.50 1 3.25 1 3.30 1 4.04 1 3.53 1 3.56 1 3.85 1 4.07 1 3.52 1 3.93 1 4.19 1 2.96 1 1.37 1 3.93 1

2.33 2.98 4.00 3.55 2.96 4.30 4.16 2.59 2.42 3.36 4.32 2.34 2.68 2.95 1.56 3.11 1.81 1.77 1.98 2.63 2.86 2.93 2.17 2.72 2.65 2.22 2.90 2.97 2.36 2.56 2.52 2.27 2.98 3.72 2.80 3.57 4.02 2.31 2.86 2.28 2.39 2.28 2.48 2.28 3.21 2.23 2.32 2.68 2.66 2.32 2.61 3.64 2.58 3.65 2.66 3.68 2.65 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3

3.02 3 3.48 3 2.42 3 2.41 3 2.66 3 3.29 3 2.70 3 3.04 3 2.81 3 1.97 3 1.68 3 0.89 4 1.06 4 1.08 4 1.27 4 1.63 4 1.89 4 1.19 4 2.17 4 2.28 4 1.72 4 1.98 4 1.74 4 2.16 4 3.37 4 2.97 4 1.69 4 0.94 4 2.11 4 2.81 4 2.52 4 1.31 4 2.51 4 1.88 4 1.41 4 3.19 4 1.92 4 2.47 4 1.02 4 2.10 4 3.71 4 ; RUN;

/*正态性检验*/

PROC UNIVARIATE NORMAL; CLASS GROUP; VAR X; RUN;

PROC ANOVA; CLASS GROUP; MODEL X=GROUP;

MEANS GROUP/HOVTEST SNK LSD DUNNETT; RUN;

/*hovtest方差齐性检验*/ /*方差分析*/ PROC GLM;

CLASS GROUP; MODEL X=GROUP;

MEANS GROUP/HOVTEST SNK LSD DUNNETT; RUN;

DATA A;

INPUT X TREAT BLOCK @@; CARDS;

0.82 1 1 0.73 1 2 0.43 1 3 0.41 1 4 0.68 1 5 0.65 2 1 0.54 2 2 0.34 2 3 0.21 2 4 0.43 2 5 0.51 3 1 0.23 3 2 0.28 3 3 0.31 3 4 0.24 3 5 ; RUN;

PROC PRINT;RUN; PROC ANOVA; CLASS TREAT BLOCK; MODEL X=TREAT BLOCK; MEANS TREAT/SNK ALPHA=0.05; MEANS TREAT; RUN;

DATA BLOCK;

INPUT BLOCK TREAT X@@; CARDS;

1 1 0.82 1 2 0.65 1 3 0.51 2 1 0.73 2 2 0.54 2 3 0.23 3 1 0.43 3 2 0.34 3 3 0.28 4 1 0.41 4 2 0.21 4 3 0.31 5 1 0.68 5 2 0.43 5 3 0.24 ;

PROC ANOVA;

CLASS BLOCK TREAT; MODEL X=BLOCK TREAT;

MEANS TREAT/HOVTEST SNK LSD DUNNETT; RUN; ;

PROC GLM;

CLASS BLOCK TREAT; MODEL X=BLOCK TREAT;

MEANS TREAT/HOVTEST SNK LSD DUNNETT; RUN;

/*roc曲线分析*/ data ex21;

input group num;