《sas程序总结》由会员分享，可在线阅读，更多相关《sas程序总结（8页珍藏版）》请在金锄头文库上搜索。
1、1计数资料：RXC 表检验：freq 卡方检 验（do r 行 do c 列；table r*c/chisq noprecent nocol exact trend chm ;weight freq ） 回归方程：plot 散点图（plot x1*x2= *） ，model x2=x1/p; output out=r R=RES; corr 相关系数（var with）, reg 回归方程（model y=x/noint ,plot y*x） 方差分析：单因素方差分析： 正态性检验 univariate data=dat5_1 normal;var x; by group;方差齐性检验 glm
2、 方差分析（class group;model x=group ;means group/hovtest）均数间多重比较（means group/hovtest;means group/snk bon dunnett2;means group /snk ALPHA=0.01;contrast 1 2 vs 3group -0.5 -0.5 1;contrast1vs 2 group 1 -1 0;）随即区组设计 do row ;do colum; 总体均数估计和假设检验：总体区间 估计（proc means data=数据集名 n mean std clm）t 检验 （d=x- ？；proc
3、means data=数据集名 n mean std stderr t prt）配对设计样本（同上）成组设计的两样本均 数比较（proc ttest data=数据集名； class group;）先看方差齐性检验的结 果，然后再根据齐不齐判断 P 值大小。计量资料的统计描述：means 过程 （proc means data=数据集名 n mean std stdrr cv min max;）Proc sort 将数 据排序 by group;Univariate data= plot freq normal;2p228 配对设计资料的符号检验配对设计资料的符号检验 （ch9_1） data
4、 dat1;input a b;d=a-b; cards; 1 0 1 0 1 0 0 1 1 0 1 0 1 0 1 0 1 0 ; proc univariate data=dat1;var d; run;P230 配对设计资料的符号和检验配对设计资料的符号和检验 (ch9_2) data dat2;input before after;d=before-after; cards; 336 258 371 291 386 300 364 285 377 298 292 303 288 312 304 260 333 339 302 290 ; proc univariate data=da
5、t2 normal;var d; run;p235 两独立样本秩和检验两独立样本秩和检验 （ch9_4）团体）团体t检验检验 data dat4;do group=1 to 2;input n;do i=1 to n;input x ;output;end;end; cards; 8 0.32 0.47 0.57 2.21 0.64 3.08 0.67 2.13 10 0.26 0.24 0.59 0.37 0.58 0.21 0.33 0.42 0.67 0.45 ; run; proc univariate data=dat4 normal; class group; var x; run
6、; proc npar1way data=dat4 wilcoxon;class group;var x; run;p236 未配对两样本比较未配对两样本比较(ch9_5) data dat5;do rank=1 to 20;input group;output;end; cards; 1 1 2 1 1 2 2 1 2 1 2 2 1 1 1 2 2 1 2 23; run; proc npar1way data=dat5 wilcoxon;class group;var rank; run;p242 等级资料（显效等级资料（显效）两样本比）两样本比 较的秩和检验较的秩和检验(ch9_7)
7、data dat7;do group=1 to 2;do rank=1 to 4;input freq ;do i=1 to freq;output; end;end;end; cards; 4 15 32 13 13 26 21 8 ; proc npar1way data=dat7 wilcoxon;class group;var rank; run;p 244 完全随机设计多个独立样本完全随机设计多个独立样本 秩和检验（秩和检验（再障患者 CD8 水平） (ch9_8) data dat8;do group=1 to 3;input n;do i=1 to n;input x ;outp
8、ut;end;end; cards; 12 293 409 392 244 213 409 57 97 244 254 352 168 10 441 538 390 589 244 409 72 168 254 374 9 807 833 409 914 380 883 254 993 667 ; proc glm data=dat8; class group; model x=group; means group/hovtest=levene; run; proc npar1way data=dat8 wilcoxon;class group;var x; run; proc rank da
9、ta=dat8 out=a; var x; ranks r; run; proc glm data=a; class group; model r=group; means group/snk; run;p250 随机单位组设计秩和检验（伤随机单位组设计秩和检验（伤 口）口）(ch9_10) data dat10;do block=1 to 12;do treat=1 to 3;input x ;output;end;end;cards;48.02 71.9066.2752.70 56.3560.5960.22 70.0866.1244.49 86.60 55.3649.31 68.2553.
10、39446.23 63.3652.3455.16 66.1255.1642.48 70.0258.6450.84 66.9744.0139.38 67.0552.4945.16 69.8959.9953.47 61.0861.08;proc freq data=dat10;tables block*treat*x/noprint scores=rank cmh2;run;p267 秩相关秩相关 Spearman 等级相关分析等级相关分析 碘碘 甲状腺肿甲状腺肿 (ch9_17) data dat17;input x y ;cards;1.1 41.22.3 37.62.5 38.83.7 24
11、.33.8 23.04.0 32.44.6 16.14.9 18.47.6 7.18.1 8.08.2 9.18.4 4.28.8 5.318.7 0.223.5 0.0;proc corr data=dat17 spearman;var x;with y;run;(ch3_2) p54 means 过程（过程（n mean std cv min max） data shg;input x ;cards;proc means data= shg n mean std cv min max;var x;run;(ch3_3) VA VB分组分组 bygroup P 55data a;input
12、group VA VB1;cards;1 1.8 1.4 2 1.7 1.1 1 2.2 1.5 3 1.9 1.2 2 2.5 1.0 1 2.7 1.6 2 2.3 1.3 2 2.8 0.9 3 3.0 1.1 1 2.6 1.4 1 2.4 1.2 2 1.9 1.3 3 2.9 0.8 1 3.2 1.7 3 3.1 1.5 2 2.6 1.9 3 3.5 1.6 3 3.3 1.5;proc sort data=a ;by group;proc means mean std max min maxdec=3;by group;var VA VB1;run;（ch3_4）p56 平均
13、抗体滴度平均抗体滴度 data a;input f x ;y=log10(x);cards;1 4 3 8 8 16 13 32 21 64 9 128 4 256 1 512;proc means noprint;var y;freq f;5output out=b mean=meany;run;data c;set b;meanx=10*(meany);run;proc print;run;（ch4_1 ）p85 总体均数区间估计总体均数区间估计 Data ;input x ;cards; proc means data=shg n mean std clm（总体均数的置信区间）;var
14、x;run;（ch4_2）p86 样本均数与总体均数样本均数与总体均数 的的 t 检验检验 data mb;input x ;d=x-72;cards;74 73 68 75 75 82 80 69 72 74 83 72 71 74 76 79 67 73 81 70 67 70 78 69 70 72 67 74 80 66;proc means data=mb mean std stderr t prt;var x d;run;（ch4_3）配对设计）配对设计 t 检验检验 p87 data ch4_7;input after before ;d=after-before;cards;7
15、0.5564.2988.6064.0768.4445.8861.6445.2364.7350.4074.6861.5969.1551.8560.5160.1365.5964.2969.0451.93;proc means data=ch4_7 mean std stderr t prt;var d;run; （ch 4_4）成组设计的两样本均数）成组设计的两样本均数 比较的比较的 t 检验检验 p88 data hb;input group x ;cards;1 26 1 32 1 25 1 22 1 20 1 28 1 24 1 19 1 29 1 17 1 34 1 21 1 20 1 2
16、3 1 272 21 2 23 2 18 2 24 2 23 2 19 2 16 2 22 2 20 2 25 2 23 2 17 2 15 2 26 2 22;data hb;do group=1to 2;input n;do i=1 to n;input x ;output;end;end;6cards; nN;Run;proc ttest data=hb;class group;var x;run; 方差齐行检验方差齐行检验 0.05 方差齐方差齐（ch5_1）p100 正态性检验正态性检验 方差齐方差齐 行检验行检验 proc sort data=dat5_1;by group;run;proc univari
