【正文】 1d26d6ca46471a74146acc8594878980 SAS 系統(tǒng)和數(shù)據(jù)分析 電子商務(wù)系列 上海財(cái)經(jīng)大學(xué)經(jīng)濟(jì)信息管理系 IS/SHUFE Created by Page 1 of 21 第十五課 用在 DATA 步的控制語句 DATA 步的基本概念、流程和有關(guān)文件的操作語句我們前面已介紹。但我們所介紹的DATA步中的 SAS 語句都是按語句出現(xiàn)的次序?qū)γ恳粋€(gè)觀測進(jìn)行處理。有時(shí)需要對一些確定的觀測跳過一些 SAS 處理語句，或者改變 SAS 語句的處理次序，就需要用到 DATA 步中的控制語句，實(shí)現(xiàn) SAS 程序的分支、轉(zhuǎn)移和循環(huán)等改變處理次序的功能。 SAS 系統(tǒng)提供的控制語句從實(shí)現(xiàn)功能的角度看主要有以下五大類： ? 實(shí)現(xiàn)循環(huán)（ DO 語句） ? 實(shí)現(xiàn)選擇（ SELECT語句） ? 實(shí)現(xiàn)分支（ IF語句） ? 實(shí)現(xiàn)轉(zhuǎn)移（ GOTO 語句） ? 實(shí)現(xiàn)連接（ LINK語句） 一、 實(shí)現(xiàn)循環(huán)（ DO 語句） 循環(huán)程序中使用 DO 語句的主要形式有四種，如下所示： ? DO 語句的程序格式之一： IF 條件表達(dá)式 THEN DO 。 一些 SAS語句 。 END 。 ? DO 語句的程序格式之二： DO 變量 =開始值 TO 終值 BY 步長值 。 一些 SAS語句 。 END 。 ? DO 語句的程序格式之三： DO WHILE （ 條件表 達(dá)式 ） 。 一些 SAS語句 。 END 。 ? DO 語句的程序格式之四： DO UNTIL （ 條件表達(dá)式 ） 。 一些 SAS語句 。 END 。 1d26d6ca46471a74146acc8594878980 SAS 系統(tǒng)和數(shù)據(jù)分析 電子商務(wù)系列 上海財(cái)經(jīng)大學(xué)經(jīng)濟(jì)信息管理系 IS/SHUFE Created by Page 2 of 21 DO WHILE 和 DO UNTIL語句中的表達(dá)式是用括號括起來的。兩種循環(huán)程序格式的區(qū)別是，對條件表達(dá)式的判斷位置。 DO WHILE 是在循環(huán)體的開頭，而 DO UNTIL 是在循環(huán)體的結(jié)束，也就是說 DO UNTIL 至少執(zhí)行循環(huán)體中 一些 SAS 語句一次。 下面我們舉例來說明 DO 語句的使用。 1. 使用循環(huán) DO 組產(chǎn)生隨機(jī)數(shù)數(shù)據(jù)集 例如，我們需要產(chǎn)生一組均勻分布的隨機(jī)數(shù)流的數(shù)據(jù)集，程序如下： Data DoRanuni 。 seed = 20202001 。 Do I = 1 to 10 by 2 。 X1=ranuni(seed ) 。 X2=ranuni(seed ) 。 Output 。 End 。 Proc print data=DoRanuni 。 Run 。 程序中的 X1 和 X2 都采用相同種子變量值 SEED=20202001 來產(chǎn)生的均勻分布的隨機(jī)數(shù)流。在數(shù)據(jù)步 DATA中使用 DO 循環(huán)語句時(shí)常常與 OUTPUT 語句配合來產(chǎn)生數(shù)據(jù)集。 OUTPUT語句作用是把當(dāng)前的觀測輸出到正在被創(chuàng)建的數(shù)據(jù)集 DoRanuni 中。第一次順序執(zhí)行產(chǎn)生Seed、 I、 X X2 四個(gè)變量， OUTPUT 輸出后，遇到 END語句回到 DO語句，產(chǎn)生 I、 XX2 變量的第二次值， Seed 變量因?yàn)闆]有遇到 DATA語句，繼續(xù)保持原來值， DO－ END 循環(huán)結(jié) 束后， DATA 步也就結(jié)束了。均勻分布隨機(jī)數(shù)是最基本也是最重要的隨機(jī)數(shù)，其 他 分布的隨機(jī)數(shù)都可以用均勻隨機(jī)數(shù)經(jīng)過變換得到。最常用的均勻分布隨機(jī)函數(shù)是 RANUNI（ seed） ，這個(gè)函數(shù)是一個(gè)模為 231－ 1，乘子為 397204094 的素?cái)?shù)模發(fā)生器。 Seed 必須是小于模 231－ 1任何數(shù)值 的 常數(shù)。相同的 Seed 值會(huì)產(chǎn)生相同的隨機(jī)數(shù)序列數(shù)，但不同次調(diào)用隨機(jī)函數(shù)所產(chǎn)生的值通常是不同的，因此計(jì)算機(jī)所產(chǎn)生的隨機(jī)數(shù)是一種偽隨機(jī)數(shù)。這個(gè)程序中的 X1 和 X2都采用相同種子變量值 SEED=20202001 所產(chǎn)生的均勻分布的隨機(jī)數(shù)流。 SAS 系統(tǒng)提供產(chǎn)生了11 種常見分布隨機(jī)數(shù)的函數(shù)， 如 表 所示，隨機(jī)數(shù)是我們實(shí)驗(yàn)和研究問題的重要的輸入數(shù)據(jù)。因此要能編寫程序，產(chǎn)生符合要求分布的隨機(jī)數(shù)數(shù)據(jù)集。 程序運(yùn)行結(jié)果 如 圖 所示 。 1d26d6ca46471a74146acc8594878980 SAS 系統(tǒng)和數(shù)據(jù)分析 電子商務(wù)系列 上海財(cái)經(jīng)大學(xué)經(jīng)濟(jì)信息管理系 IS/SHUFE Created by Page 3 of 21 表 SAS 系統(tǒng)的各種隨機(jī)函數(shù) 隨機(jī)數(shù)函數(shù)名 作 用 UNIFORM（ seed） 產(chǎn)生（ 0， 1）區(qū)域均勻分布隨機(jī)數(shù)，乘同余發(fā)生器 RANUNI(seed) 產(chǎn)生（ 0， 1）區(qū)域均勻分布隨機(jī)數(shù)，素?cái)?shù)模發(fā)生器 NORMAL(seed) 產(chǎn)生標(biāo)準(zhǔn)正態(tài)分布隨機(jī)數(shù)，利用中心極限定理近似公式 RANNOR(seed) 產(chǎn)生標(biāo)準(zhǔn)正態(tài)分布隨機(jī)數(shù)，利用變換抽樣法 RANEXP(seed) 產(chǎn)生λ =1 的指數(shù)分布隨機(jī)數(shù) RANGAM(seed,alpha) 產(chǎn)生伽馬分布隨機(jī)數(shù)， alpha0,seed 為任意數(shù)值 RANTRI(seed,h) 產(chǎn)生三角分布隨機(jī)數(shù)， 0h1,seed 為任意數(shù)值 RANCAU(seed) 產(chǎn)生標(biāo)準(zhǔn)柯西分布隨機(jī)數(shù) RANBIN(seed,n,p) 產(chǎn)生二項(xiàng)分布隨機(jī)數(shù)， n0 的整數(shù) ,0p1,seed 為 任意數(shù)值 RANPOI(seed,lambda) 產(chǎn)生泊松分布隨機(jī)數(shù)， lambda0,seed 為任意數(shù)值 RANTBL(seed,p1,…,p2,…pn) 產(chǎn)生離散分布隨機(jī)數(shù)， 0≤ pi≤ 1,seed 為任意數(shù)值 注：種子 seed 一般取 0，或 5 位， 6 位， 7 位的奇整數(shù)。 對于均值為 M，標(biāo)準(zhǔn)差為 S 的正態(tài)分布隨機(jī)數(shù)，可由標(biāo)準(zhǔn)正態(tài)分布隨機(jī)數(shù)的線性函數(shù)得到： X=M+S*NORMAL（ seed） 2. 在循環(huán) DO 組中使用下標(biāo)數(shù)組產(chǎn)生數(shù)據(jù)集 當(dāng)我們需要用同一種方法來處理很多變量時(shí) ,可以用數(shù)組語句定義這組變量為數(shù)組的一些 元素，這個(gè)數(shù)組中 的 一些元素就可以在 DATA步中較后面的 SAS 語句里以數(shù)組下標(biāo)的形式被引用。數(shù)組 ARRAY 語句的基本格式為： 圖 用循環(huán) DO 組產(chǎn)生隨機(jī)數(shù)數(shù)據(jù)集 1d26d6ca46471a74146acc8594878980 SAS 系統(tǒng)和數(shù)據(jù)分析 電子商務(wù)系列 上海財(cái)經(jīng)大學(xué)經(jīng)濟(jì)信息管理系 IS/SHUFE Created by Page 4 of 21 Array 數(shù)組名 {下標(biāo) } $長度 數(shù)組元素 （初始值） 。 例如，以下的幾種數(shù)組定義方式都是合法的： ? Array x{3} T1 T2 T3 。 ? Array x{5,3} T1T15 。 ? Array x{2:6,2:4} T1T15 。 ? Array x{3} T1 T2 T3 (100,99,98) 。 ? Array x{*} T1 T2 T3 。 第一種方式表示定義了一個(gè)一維名為 X 的數(shù)組，它有三個(gè)元素，對應(yīng)的變量為 T1， T2和 T3。第二種方式表示定義了一個(gè)二維名為 X 的數(shù)組，它共有 5 3=15 個(gè)元素，對應(yīng)的變量為 T1到 T15。第三種方式與第二種方式 的 區(qū)別是還規(guī)定每一維下標(biāo)的下界和上界，通常不特別指明下標(biāo)的下界從 1 開始。第四種方式給出了數(shù)組中相應(yīng)元素的初始值 T1=100， T2=99，T3=98。第五種方式下標(biāo)用星號 *來代替，表示 SAS 系統(tǒng)通過數(shù)組中的變量個(gè)數(shù)來確定下標(biāo)。 下面我們通過一個(gè)在循環(huán)中使用數(shù)組變量來產(chǎn)生一個(gè)新的數(shù) 據(jù)集的例子，來說明循環(huán)中使用數(shù)組變量的用法。假設(shè)要由一個(gè)老的數(shù)據(jù)集產(chǎn)生一個(gè)新的數(shù)據(jù)集，新的數(shù)據(jù)集要新增 n個(gè)變量，新增變量 Ti的值與原數(shù)據(jù)集的變量 Xj值和新增變量的位置值 i有關(guān)。為簡便起見，假設(shè)新增變量 Ti=∑ Xji。程序如下： Data DoArray 。 Input X1X3 。 Array a{4} T1T4 。 Do i = 1 to 4 。 a{i}=(X1+X2+X3)*i 。 End 。 Card 。 1 2 3 4 5 6 Proc print dat