统计分析 SAS练习.doc

第七章 t检验单样本t检测，用means、 univariatedata temp;input x;y=x-20.7;cards;20.99 20.41 20.10 20.00 20.91 22.6020.99 20.42 20.90 22.99 23.12 20.89;proc means mean stderr t probt;var y;run;data temp;input x;y=x-20.7;cards;20.99 20.41 20.10 20.00 20.91 22.6020.99 20.42 20.90 22.99 23.12 20.89;proc univariate normal;var y;run;无原始数据的单样本t检验:data temp;input u0 u1 s n;t=(u1-u0)/(s/sqrt(n);v=n-1;p=probt(t,v);cards;72 74.2 6.5 25;proc print;run;配对的t检测：有两个对比的检测 有ttest means univariatedata temp;input id a b;cards;1 2.41 2.802 2.90 3.043 2.75 1.884 2.23 3.435 3.67 3.816 4.49 4.007 5.16 4.448 5.45 5.419 2.06 1.2410 1.64 1.8311 1.06 1.4512 0.77 0.92;proc ttest;paired a*b;run;data temp;input id a b;y=a-b;cards;1 2.41 2.802 2.90 3.043 2.75 1.884 2.23 3.435 3.67 3.816 4.49 4.007 5.16 4.448 5.45 5.419 2.06 1.2410 1.64 1.8311 1.06 1.4512 0.77 0.92;proc univariate normal;var y;run;正态性检测:w p ，t=0.162 pr=0.870.05,接受H0；data temp;input id a b;y=a-b;cards;1 2.41 2.802 2.90 3.043 2.75 1.884 2.23 3.435 3.67 3.816 4.49 4.007 5.16 4.448 5.45 5.419 2.06 1.2410 1.64 1.8311 1.06 1.4512 0.77 0.92;proc means mean stderr t probt;var y;run;data temp;input d s n;t=d/(s/sqrt(n);v=n-1;p=probt(t,v);cards;0.625 0.78 8;proc print;run;独立样本t检测data temp;input group x;cards;1 134 1 146 1 104 1 119 1 124 1 161 1 107 1 83 1 113 1 129 1 97 1 1232 70 2 118 2 101 2 85 2 107 2 132 2 94proc ttest;class group;var x;run;data temp;input x1 x2 n1 n2 s1 s2;sc=sqrt(1/n1+1/n2)*(n1-1)*s1*2+(n2-1)*s2*2)/(n1+n2-2);t=(x2-x1)/sc;v=n1+n2-2;p=probt(t,v);cards;20.95 21.79 20 20 5.89 3.43;proc print;run;第九章 直线回归和相关data exam9_1;input x1 x2 ;cards;1.21 3.90 1.30 4.50 1.39 4.20 1.42 4.83 1.47 4.161.56 4.93 1.68 4.32 1.72 4.99 1.98 4.70 2.10 5.20;proc gplot; /*拓扑图*/plot x2*x1=*; /*图*/run;proc corr;var x1 x2;run;相关性 相关系数为：0.68073。 p为0.0303 data exm9_2;input x y;cards;79 45 80 30 91 16 90 24 70 28 87 25 92 14;proc corr spearman;var x y;run;直线回归data exm9_3;input x y;cards;25.5 9.2 19.5 7.8 24.0 9.420.5 8.6 25.0 9.0 22.0 8.821.5 9.0 23.5 9.4 26.5 9.723.5 9.4 26.5 9.7 23.5 8.822.0 8.5 20.0 8.2 28.0 9.9;proc reg;model y=x/r clm; /*针对每一条的观测的输出*/run;多元线性回归data temp;input x1 x2 x3 y;cards;51.3 73.6 36.4 2.9948.9 83.9 34.0 3.1142.8 78.3 31.0 1.9155.0 77.1 31.0 2.6345.3 81.7 30.0 2.8645.3 74.8 32.0 1.9151.4 73.7 36.5 2.9853.8 79.4 37.0 3.2849.0 72.6 30.1 2.5253.9 79.5 37.1 3.27;proc reg;model y=x1 x2 x3/tol vif collin selection=stepwise r;run;偏相关与复相关data exm10_2;input x1 x2 x3 x4 y;cards;7 26 6 60 78.51 29 15 52 74.3;proc corr;var x1 x2 x3 x4 y;run;proc corr nosimple;var x1 y;partial x2 x3 x4;run;proc corr nosimple;var x2 y;partial x1 x3 x4;run;proc corr nosimple;var x3 y;partial x2 x1 x4;run;proc corr nosimple;var x4 y;partial x2 x3 x1;run;proc corr nosimple;var x1 x2 x3 y;partial x4;run;data exm10_2;input x1 x2 x3 x4 y;cards;7 26 6 60 78.51 29 15 52 74.311 56 8 20 104.311 31 8 47 87.67 52 6 33 95.9;proc corr;var x1 x2 x3 x4 y;run;proc corr nosimple;var x1 y;partial x2 x3 x4;run;proc corr nosimple;var x2 y;partial x1 x3 x4;run;proc corr nosimple;var x3 y;partial x2 x1 x4;run;proc corr nosimple;var x4 y;partial x2 x3 x1;run;proc corr nosimple;var x1 x2 x3 y;partial x4;run;第十四章 非参数检验data exm10_2;input x1 x2 x3 x4 y;cards;7 26 6 60 78.51 29 15 52 74.311 56 8 20 104.311 31 8 47 87.67 52 6 33 95.9;proc corr;var x1 x2 x3 x4 y;run;proc corr nosimple;var x1 y;partial x2 x3 x4;run;proc corr nosimple;var x2 y;partial x1 x3 x4;run;proc corr nosimple;var x3 y;partial x2 x1 x4;run;proc corr nosimple;var x4 y;partial x2 x3 x1;run;proc corr nosimple;var x1 x2 x3 y;partial x4;run;data exm4_2;input wt;devi=wt-2.15;cards;2.15 2.10 2.20 2.1.22.42 2.52 2.62 2.722.99 3.19 3.37 4.57;proc univariate; /单变量/var devi;run;P0.05，小概率事件发生，拒绝H0，接受H1。配对的非参数的检验：data exm10_2;input x1 x2 x3 x4 y;cards;7 26 6 60 78.51 29 15 52 74.311 56 8 20 104.311 31 8 47 87.67 52 6 33 95.9;proc corr;var x1 x2 x3 x4 y;run;proc corr nosimple;var x1 y;partial x2 x3 x4;run;proc corr nosimple;var x2 y;partial x1 x3 x4;run;proc corr nosimple;var x3 y;partial x2 x1 x4;run;proc corr nosimple;var x4 y;partial x2 x3 x1;run;proc corr nosimple;var x1 x2 x3 y;partial x4;run;data exm4_2;input wt;devi=wt-2.15;cards;2.15 2.10 2.20 2.1.22.42 2.52 2.62 2.722.99 3.19 3.37 4.57;proc univariate; /单变量/var devi;run;data ex14_4;input x1 x2;d=x1-x2;cards;39 5542 5451 5543 4755 5345 6322 5248 4440 4845 5540 3249 57;proc univariate;var d;run;两种非参数检验data ex14_6;input group veci;cards;1 2.78 1 3.23 1 4.2 1 4.87 1 5.121 6.21 1 7.18 1 8.05 1 8.56 1 9.62 3.23 2 3.5 2 4.04 2 4.15 2 4.282 4.34 2 4.47 2 4.64 2 4.75 2 4.822 4.95 2 5.10;proc nparlway wilcoxon;class group;var veci;run;data ex14_9;do group=1 to 3;input rate;output; end;cards;32.5 16.0 6.535.5 20.5 9.040.5 22.5 12.546.0 29.0 18.049.0 36.0 24.0;proc nparlway wilcoxon;class group;var rate;run;data ex14_12;do x=1 to 4;do g=1 to 2;input f;output;end;end;cards;65 4218 630 2313 11;proc nparlway wilcoxon;calss g;var x;freq f ;run;data ex14_12;do x=1 to 4;do g=1 to 4;input f;output;end;end;cards;0 3 5 32 5 7 59 5 3 36 2 2 0 ;proc nparlway wilcoxon;calss g;var x;freq f ;run;data ex14_17;input x y ;cards;79 45 80 30 91 16 90 2470 28 87 25 92 14;proc corr spearman;var x y;run;data exm14_15;do block=1 to 8;do treat=1 to 4;input rate;output;end; end;cards;8.4 9.6 9.8 11.711.6 12.7 11.8 12.09.4 9.1 10.4 9.89.8 8.7 9.9 12.08.3 8.0 8.6 8.68.6 9.8 9.6 10.68.9 9.0 10.6 11.47.8 8.2 8.5 10.8;proc sort data=exm14_15 out=a;by block;run;proc rank data=a out=b;var rate;by block;ranks r;proc glm data=b;class treat block;model r=treat block