【文章內(nèi)容簡介】 9。女 39。)。 /*先生成一個取值為 0或 1的值，性別為女時為 1，否則為 0*,然后賦值給變量 isffem*/ newv = .。 /*將變量賦了缺失值 */ 注意想試驗上述語句要把它們放入數(shù)據(jù)步中，并且等號右邊的表達式中的各變量應該是存在的，否則會得到缺失值結(jié)果。 SAS數(shù)據(jù)步的輸出一般是數(shù)據(jù)集，用賦值語句計算的結(jié)果會自動寫入數(shù)據(jù)集。 SAS也提供了一個 PUT語句，可以象其它語言程序的 PRINT、 WRITE(*,*)、 printf等語句一樣立即打印輸出結(jié)果。 PUT的基本語法很多 ,只舉幾例： data。 x=。 y=sin(x)。 put 39。Sine function value of 39。 x 39。is 39。 y。 run。 結(jié)果將在運行記錄窗口顯示一行 Sine function value of is 在 PUT語句中使用“變量名＝”來指定輸出項可以顯示帶有變量名的輸出結(jié)果，比如把上程序中的 PUT語句改為 put x= y=。 則結(jié)果在 LOG窗口顯示為 X= Y= 分支結(jié)構(gòu) 如果需要在某條件滿足時執(zhí)行某一操作，可以用 IF 條件 THEN 程序塊 ELSE 程序塊 的結(jié)構(gòu)，其中 Else語句及其后的程序塊可以省略。需要注意的是，如果程序塊只有一句，則可直接寫出，否則應以 DO開始 ，以 END結(jié)束 。 比如， if max(a)20 then p=0。 else p=1。 上述程序當 max(a)20時 p=0，否則 p=1。 下述程序則是復合的 IF結(jié)構(gòu)。 IF x0 THEN DO。 PUT 39。 X為正數(shù) 39。 x = 2*x。 PUT x=。 END。 注意 SAS的分支結(jié)構(gòu)的寫法與其它語言有些不同，它不用 ENDIF結(jié)束。 ? 循環(huán)結(jié)構(gòu) SAS數(shù)據(jù)步可以使用循環(huán)結(jié)構(gòu)，主要的是兩種：計數(shù) DO循環(huán)和當型、直到型循環(huán)。 計數(shù) DO循環(huán)的寫法是： DO 計數(shù)變量 ＝ 起始值 TO 結(jié)束值 BY 步長 。 循環(huán)體語句 …… END。 在循環(huán)體中可以用 LEAVE語句跳出循環(huán)，相當于 C語言的 break語句 。用 CONTINUE語句可以立即結(jié)束本輪循環(huán)并轉(zhuǎn)入下一輪循環(huán)的判斷與執(zhí)行 . 當型循環(huán) 的語法是： DO WHILE(循環(huán)繼續(xù)條件 )。 循環(huán)體語句 …… END。 程序先判斷循環(huán)繼續(xù)條件是否成立，成立時執(zhí)行循環(huán)體語句，再判斷循環(huán)繼續(xù)條件，如此重復，直到循環(huán)繼續(xù)條件不再成立。 直到型循環(huán) 的寫法是： DO UNTIL (循環(huán)退出條件 )。 循環(huán)體語句 …… END。 程序先執(zhí)行循環(huán)體，然后判斷循環(huán)退出條件是否成立，成立則結(jié)束循環(huán)，否則繼續(xù)。注意每輪循環(huán)都是先執(zhí)行循環(huán)體再判斷是否退出。 如： data。 DO i = 1 TO 20 BY 2。 j = i**3。 put i 3. j 5.。 END。 run。 可以輸出一個 1， 3， 5， 7， … ， 19的立方表。 七、 SAS的基本統(tǒng)計分析功能 －－－ SAS Analysis 模塊 ? 前面我們已經(jīng)對 SAS系統(tǒng)進行了簡要的介紹，下面我們講如何利用 SAS Analysis 模塊計算描述統(tǒng)計量、進行假設檢驗、擬合曲線以及進行回歸分析 啟動 Analyst模塊 ? 選擇 Solutions菜單下的 Analysis，然后再選擇下面的 Analyst，如圖： ? 得到 Analyst界面 ? 數(shù)據(jù)集相關操作 ? 一、依分布產(chǎn)生數(shù)據(jù) ? 利用 Data菜單中 Random Variates可以產(chǎn)生服從特定分布的隨機數(shù)，可以選擇的分布類型包括： ? 正態(tài)分布 ? 均勻分布 ? 二項分布 ? 卡方分布 ? 泊松分布 ? 指數(shù)分布 ? Beta分布 ? Gama分布 ? 幾何分布 ? 極值分布 ? 例 利用 Analyst模塊生成 100個服從標準正態(tài)分布的數(shù)據(jù)，并把生成的數(shù)據(jù)命名為 Normal保存在Sasuser下。 ? 操作步驟： （ 1）在 Data菜單中選擇 Random Variates，然后選擇分布類型 —— Normal（正態(tài)分布），則彈出現(xiàn)一對話框： （ 2）在 Number of values to generate中填 100，在 New column name (變量名稱）中填 X，在Mean中填 0，在 Standard deviation 中填 1，單擊OK ? （ 3）可以看到數(shù)據(jù)集中產(chǎn)生了新的一列，在 File菜單中選 Save As By SAS Name，保存數(shù)據(jù)集 （ 4）在彈出的對話框中雙擊 Sasuser，并在 Member Name中填數(shù)據(jù)集名稱 Normal，單擊 Save即可 二、隨機抽樣 使用 Data中 Random Sample可以從數(shù)據(jù)集中隨機抽取記錄。 例：在剛才生成的正態(tài)數(shù)據(jù)集 normal中隨機抽取 50個數(shù)據(jù)。 步驟如下： （ 1）首先選擇 Data菜單下 Random Sample命令，彈出對話框： （ 2）在 Rows中填 50，單擊 OK即可得到一個表格含有隨機抽取的 50個樣本，如圖： 繪制統(tǒng)計圖 一、繪制變量直方圖并在圖上擬合分布曲線 利用 Grphs中 Histogram可以繪制變量直方圖并擬合分布曲線。 例：以上節(jié) normal數(shù)據(jù)集中 X為對象繪制直方圖，并在直方圖上擬合正態(tài)分布曲線。 操作步驟： （ 1）打開 Normal文件： File→Open By SAS Name→Sasuser→Normal→OK （ 2）在 Grphs下選 Histogram，單擊繪圖變量 X，再單擊 Analysis將其加入 （ 3）選擇 Fit在直方圖上擬合正態(tài)分布曲線，單擊 Fit彈出對話框 ? （ 3）選擇 Normal，單擊 OK回到原對話框，再單擊OK即得： 二、概率圖 在直方圖基礎上我們可以猜想變量是服從正態(tài)分布的，通過繪制概率圖可以來證實我們的想法。概率圖將數(shù)值排序，給出每個數(shù)值對應的分位數(shù)，然后打點作圖。如果這些點呈現(xiàn)線性特征，說明他們與理論分布相符，同時又在圖像上加上一條給定分布的曲線，并給出分布的系數(shù)。 例：繪制 Normal數(shù)據(jù)集中 X的概率圖，步驟如下： （ 1）首先選擇 Graphs下 Probality Plot，彈出對話框 （ 2）單擊 X，再單擊 Analysis，然后選擇分布為Normal，單擊 Ok得到概率圖： 可以看到 X的分布與正態(tài)分布基本吻合，圖中還給出了正態(tài)分布的參數(shù)：均值為－ ，標準差為 。 一、繪制散點圖 使用 Graphs下的 Scatter Plot可以繪制散點圖 例：病人在服用某種放射性藥物后，體內(nèi)的放射性物質(zhì)的含量 y與時間 x具有一定的關系，繪制二者的散點圖，然后找出他們的關系，數(shù)據(jù)如下： x 1 2 3 4 5 6 7 8 9 10 11 y 37 14 5 首先將數(shù)據(jù)輸入建立一個名為 Test的數(shù)據(jù)集，然后在 ANALYST模塊下打開 然后打開 Scatter Plot對話框 在繪制散點圖對話框中將 x放在 X軸， y放在 Y軸，單擊 OK可以得到二者的散點圖，為容易觀察，可在Display中選擇用線連接散點圖 繪制的圖形如下，可以看到放射性物質(zhì)的含量和時間具有負指數(shù)衰減關系 統(tǒng)計分析與計算 一、計算描述性統(tǒng)計量 利用 Statistics 的 Descriptive菜單可以計算 描述性統(tǒng)計量， Descriptive下有 4個選項： Summary Statistics計算簡單統(tǒng)計量 Distribution 計算分布信息 Correlations 計算變量間的相關系數(shù) Frequency Counts統(tǒng)計變量的頻數(shù) 例 1：計算