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