《農(nóng)業(yè)試驗(yàn)統(tǒng)計(jì)》SAS軟件使用方法簡(jiǎn)介

1、SAS軟件使用方法簡(jiǎn)介 第二十講1西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作SAS的特點(diǎn)實(shí)用性強(qiáng)功能完善使用方便簡(jiǎn)單易學(xué)2西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作系統(tǒng)的啟動(dòng)將系統(tǒng)時(shí)間調(diào)至1998年6月15 日前 從桌面上找到“The SAS system for windows v6.12”，雙擊啟動(dòng)3西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作管理系統(tǒng)(PROGRAM EDITOR)窗口窗口窗口程序的輸入、修改、運(yùn)行、讀入、打印和存儲(chǔ)等操作存放程序運(yùn)行過(guò)程，包括程序語(yǔ)句、注釋、警告、錯(cuò)誤顯示等信息存放程序運(yùn)行產(chǎn)生的結(jié)果4西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作程 序一個(gè)SAS程序一般包括兩個(gè)部分：

2、 數(shù)據(jù)步 (DATA step) 建立SAS數(shù)據(jù)集 過(guò)程步 (PROC step) 調(diào)用一個(gè)或多個(gè)程序模塊完成數(shù)據(jù)分析5西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作數(shù)據(jù)步 (DATA step)Data exam2;INPUT sex$ x y;CARDS;男 18 178 女 18 167 女 21 165 女 17 158男 20 180;Data exam1;INPUT sex$ x y;CARDS;男 18 178 女 18 167女 21 165女 17 158男 20 180;6西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作過(guò)程步 (PROC step)PROC 程序模塊；RUN; PROC

3、 程序模塊；PROC 其它程序模塊；RUN; 7西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作常用的過(guò)程Sort: 對(duì)給定變量進(jìn)行排序Print: 輸出SAS數(shù)據(jù)集中的觀察值Tabulate: 制作統(tǒng)計(jì)報(bào)表Means: 對(duì)數(shù)值變量計(jì)算簡(jiǎn)單的描述統(tǒng)計(jì)量Freq: 對(duì)分類變量進(jìn)行頻數(shù)分布計(jì)算Ttest: 對(duì)變量進(jìn)行T檢驗(yàn)8西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作Anova: 對(duì)指定變量進(jìn)行方差分析Glm: 對(duì)指定變量進(jìn)行方差分析、協(xié)方差分析和 回歸分析Reg: 對(duì)指定變量進(jìn)行回歸分析Corr: 計(jì)算變量間的相關(guān)系數(shù)Plan: 進(jìn)行試驗(yàn)設(shè)計(jì)Nlin: 對(duì)指定變量進(jìn)行非線性回歸分析9西南科技大學(xué)生命科學(xué)與工

4、程學(xué)院周海廷制作描述統(tǒng)計(jì)常用的過(guò)程 MEANS CHART UNIVARIATE10西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作PROC MEANS 選擇項(xiàng)計(jì)算描述統(tǒng)計(jì)量的MEANS 過(guò)程VAR 變量表BY 變量表CLASS 變量WEIGHT 變量OUTPUT 選擇項(xiàng)FREQ 變量1.N 13.STDERR2.NMISS 14.T3.MEAN 15.PRT4.STD 16.SUMWGT5.MIN 17.SKEWNESS6.MAX 18.KURTOSIS7.RANGE8.SUM9.VAR10.USS11.CSS12.CV11西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作PROC CHART 選擇項(xiàng)打印各種

5、統(tǒng)計(jì)圖的CHART 過(guò)程VAR 繪圖變量表畫縱向直方圖HBAR 繪圖變量表畫水平直方圖BLOCK 繪圖變量表畫三維直方圖PIE 繪圖變量表畫餅圖STAR 繪圖變量表畫星座圖12西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作PROC UNIVARIATE 選擇項(xiàng)描述統(tǒng)計(jì)量的UNIVARIATE 過(guò)程VAR 變量表BY 變量表FREQ 變量WEIGHT 變量OUTPUT OUT=1.DATA=SAS數(shù)據(jù)集2.PLOT 做莖葉圖、盒式圖、正態(tài)概率圖3.FREQ 做頻數(shù)分布表4. 進(jìn)行正態(tài)性檢驗(yàn)13西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作 /*LT3*/OPTION PS=30 LS=100;DATA A;

6、INPUT YIELD;YIELD=INT(YIELD-67.5)/13)*13+67.5;CARDS;177 215 197 97 123 159 245 119 119 131 149 152167 104 161 214 125 175 219 118 192 176 175 95136 199 116 165 214 95 158 83 137 80 138 151187 126 196 134 206 137 98 97 129 143 179 174159 194 136 108 101 141 148 168 163 176 102 158145 173 75 130 149 1

7、50 161 155 111 205 131 189 定義輸出格式:每頁(yè)30行,每行100個(gè)字符14西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作91 142 140 154 152 163 123 187 149 155 131 209183 97 119 181 149 254 131 215 111 186 118 150155 197 116 184 239 160 172 179 151 198 124 179135 151 168 169 173 181 188 211 197 175 122 159171 166 175 143 190 213 192 231 163 159 158

8、159177 147 194 227 141 169 124 165;PROC MEANS MEAN STD CV STDERR MAXDEC=2;PROC CHART;VAR YIELD;PROC FREQ;PROC UNIVARIATE NORMAL PLOT;RUN;數(shù)值保留兩位小數(shù)15西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作統(tǒng)計(jì)推斷常用的過(guò)程 MEANS單個(gè)樣本平均數(shù)及兩個(gè)樣本平均數(shù)成對(duì)比較 TTEST兩個(gè)樣本平均數(shù)成組比較PROC TTEST ；CLASS 變量表；VAR 變量表；PROC MEANS T PRT N MEAN STD ；RUN；16西南科技大學(xué)生命科學(xué)與工程學(xué)院周海

9、廷制作 /*LT5-1*/ DATA A;INPUT WEIGHT;DIEF=WEIGHT-34;CARDS;35.6 37.6 33.4 35.1 32.7 36.8 35.9 34.6;PROC MEANS T PRT MAXDEC=2 STD;VAR DIEF;RUN;17西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作 /*LT5-6*/DATA T1;INPUT Y1 Y2;DIFF=Y1-Y2;CARDS;10 25 13 12 8 14 3 15 5 12 20 27 6 18;PROC MEANS T PRT MAXDEG=2;VAR DIFF;RUN;18西南科技大學(xué)生命科學(xué)與工程學(xué)

10、院周海廷制作 /*LT5-4*/DATA X;INPUT D GR$;CARDS;160 A 170 B 160 A 270 B 200 A 180 B160 A 250 B 200 A 270 B 170 A 290 B150 A 270 B 210 A 230 B 170 B;PROC TTEST;CLASS GR;VAR D;RUN;19西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作進(jìn)行方差分析的ANOVA 過(guò)程PROC ANOVA 選擇項(xiàng)FREQ 變量表MODLE 因變量=自變量CLASS 變量BY 變量表MEANS 變量/選項(xiàng)CLASS ACLASS A BMODEL Y=AMODEL Y

11、=A BMODEL Y=A B A*BY/LSD Y/DUNCAN Y/ SNKMODEL y=i j(i)20西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作/*LT6-1*/DATA Aa;DO treat=1 To 4; DO r=1 to 4; INPUT y; OUTPUT; END;END;CARDS;18 21 20 13 20 24 26 22 10 15 17 14 28 27 29 32;PROC ANOVA;CLASS treat;MODEL y=treat ;MEANS treat/LSD (DUNCAN、SNK) ;run;21西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作/*LT

12、6-11*/DATA Aa;DO i=1 To 4; DO j=1 to 8; INPUT y; OUTPUT; END;END;CARDS;12 13 14 15 15 16 17 . 14 10 11 13 14 11 . . 9 2 10 11 12 13 12 11 12 11 10 9 8 10 12 .;PROC ANOVA;CLASS i;MODEL y=i;MEANS i/LSD;run;22西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作/*LT6-12*/data aa;do i=1 to 4; do j=1 to 3; do l=1 to 4; input y; output;

13、end; end; end;cards;50 55 40 35 35 35 30 40 45 40 40 5023西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作50 45 50 45 55 60 50 50 55 45 65 5585 60 90 85 65 70 80 65 70 70 70 7060 55 35 70 60 85 45 75 65 65 85 75;proc anova;class i j;model y=i j(i);means i /duncan e=j(i) alpha=0.01;run;24西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作用REG過(guò)程進(jìn)行回歸分析PROC REG DATA=輸入數(shù)據(jù)集 選項(xiàng); VAR 可參與建模的變量列表; MODEL 因變量自變量表 / 選項(xiàng); PRINT 輸出結(jié)果; PLOT 診斷圖形; RUN; 25西南科技大學(xué)生命科學(xué)與工程學(xué)院周海廷制作 /*LT9-1*