統(tǒng)計(jì)分析SAS練習(xí)

1、第七章 t檢驗(yàn)單樣本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 nor

2、mal;var y;run;無原始數(shù)據(jù)的單樣本t檢驗(yàn):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檢測：有兩個(gè)對比的檢測 有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

3、.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;正態(tài)性檢測:w p ，t=0.162 pr=0.870.05,接受H0；d

4、ata 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

5、;獨(dú)立樣本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

6、);cards;20.95 21.79 20 20 5.89 3.43;proc print;run;第九章 直線回歸和相關(guān)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; /*拓?fù)鋱D*/plot x2*x1=*; /*圖*/run;proc corr;var x1 x2;run;相關(guān)性 相關(guān)系數(shù)為：0.68073。 p為0.0303 data exm9_2;i

7、nput 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;多元線性

8、回歸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;偏相關(guān)與復(fù)相關(guān)d

9、ata 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

10、 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 n

11、osimple;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;第十四章 非參數(shù)檢驗(yàn)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 cor

12、r 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.42 2.5

13、2 2.62 2.722.99 3.19 3.37 4.57;proc univariate; /單變量/var devi;run;P0.05，小概率事件發(fā)生，拒絕H0，接受H1。配對的非參數(shù)的檢驗(yàn)：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

14、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.42 2.52 2.62 2.722.99 3.19 3.37 4.57;proc univariate

15、; /單變量/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;兩種非參數(shù)檢驗(yàn)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

16、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

17、;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

18、 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