【正文】 使用 CALL語句時(shí)隨時(shí)可以看到當(dāng)前的種子值 ， 用隨機(jī)數(shù)函數(shù)則看不到 。 例中 ， 產(chǎn)生兩個(gè)獨(dú)立的種子流和隨機(jī)數(shù)流 。 proc print。 output。 call ranuni(seed3, X3)。 retain seed3 seed4 135279821。 結(jié)果顯示： Obs seed1 seed2 I x1 x2 1 161321804 161321804 1 2 161321804 161321804 2 3 161321804 161321804 3 4 161321804 161321804 4 5 161321804 161321804 5 例 使用 CALL子程序產(chǎn)生兩個(gè)獨(dú)立的種子流和隨機(jī)數(shù)流 。因?yàn)椋?SEED2的值在這里不起作用， X2的第一個(gè)值并不是由SEED2產(chǎn)生的，而是產(chǎn)生第一個(gè) X1后一個(gè)新種子的結(jié)果。 run。 options nocenter。 output。 x1= ranuni(seed1)。 retain seed1 seed2 161321804。 R ou t i ne S A S 隨機(jī)子程序名S e e d 存放當(dāng)前種子值的變量A r gu me nt 特殊分布要求的參數(shù)V a r i a t e 存放生成隨機(jī)數(shù)的變量名例 使用隨機(jī)函數(shù)產(chǎn)生兩個(gè)隨機(jī)數(shù)變量屬于同一個(gè)隨機(jī)數(shù)流 。 同時(shí)產(chǎn)生幾個(gè)隨機(jī)數(shù)流時(shí)，用 CALL子程序比用隨機(jī)函數(shù)的效果更好。 CALL語句激活隨機(jī)數(shù)子程序的格式為： CALL routine (seed, argument, variate)。 隨機(jī)數(shù)子程序 用隨機(jī)數(shù)子程序可以更好地控制種子流和隨機(jī)數(shù)流 。如果要控制隨機(jī)數(shù)流，就要用隨機(jī)數(shù)函數(shù)的 CALL子程序。 表 SEED取值與初始化類型 自變量 S E E D 初始化類型≤ 0 等價(jià)于用計(jì)算機(jī)的日期時(shí)間值作為初始種子值來初始化隨機(jī)數(shù)流＞ 0 用自變量 S E E D 的值作為當(dāng)前的初始種子值來初始化隨機(jī)數(shù)流產(chǎn)生隨機(jī)數(shù)的過程中，自變量 SEED的值保持不變，種子則不斷變化。 隨機(jī)數(shù)函數(shù)自變量 SEED 隨機(jī)數(shù)函數(shù)使用一個(gè)自變量 SEED來選擇產(chǎn)生隨機(jī)數(shù)的初始種子值，由這個(gè)值開始產(chǎn)生隨機(jī)數(shù)流。 drop _seed_ _i_。 /*區(qū)間 [a, b]上的均勻分布 */ output。 do _i_ = 1 to 1000。 a=1。 data RV1。 run。 end。 uniform1 = a + (ba) * ranuni(_seed_)。 b=2。 retain _seed_ 0。 均勻分布 例 用函數(shù) RANUNI產(chǎn)生均勻分布隨機(jī)數(shù) 。 drop _seed_ _i_。/*均值為 mu,標(biāo)準(zhǔn)差為 sigma*/ output。 do _i_ = 1 to 1000。 mu=0。 data RV。 run。 end。 Normal1 = mu + sigma * rannor(_seed_)。 sigma=1。 retain _seed_ 0。 正態(tài)分布 例 用函數(shù) RANNOR產(chǎn)生正態(tài)分布隨機(jī)數(shù) 。 隨機(jī)數(shù)函數(shù) SAS系統(tǒng)提供 11種隨機(jī)數(shù)函數(shù)。 方差 VAR(of xlxn) 或 VAR (x1,x2, … ) 計(jì)算自變量中非缺失值的方差。 求和 SUM(of xlxn) 或 SUM (x1,x2, … ) 計(jì)算自變量的和。 均值 MEAN(of xlxn) 或 MEAN(x， y， z, … ) 計(jì)算非缺失自變量的算術(shù)平均。 一般情況下 ， 求變量 X的樣本統(tǒng)計(jì)函數(shù)值時(shí)要用別的方法 。 run。 P3=probnorm()。 P1=probnorm(0)。 例 計(jì)算三個(gè)特殊的正態(tài)概率值 。 run。 y=15263*24*3600+x。 dhms=dhms (15263,21,50,51)。 dhms=dhms (date (), 21,50,51)。 x=21*60*60+50*60+51。 hms=hms (21,50,51)。 put mdy=。 data。 run。 day=day (date)。 qtr=qtr (date)。 keep date year qtr month day。 data a。 run。 format x y yymmdd10.。 x=date ()。 當(dāng)然 ， 也可以用一定的日期格式表示它們 。 run。 y=today()。 data。 例中 ， 兩段程序的結(jié)果一樣 。 put actual= days360=。 actual=datdif(‘01jan2021’d, ‘01jan2021’d, ‘a(chǎn)ct/act’)。 run。)。 /*按每個(gè)月的實(shí)際天數(shù)算 */ days360=datdif(sdate, edate, 39。act/act39。d。 edate=39。01jan202139。 data _null_。 ? 由于日期函數(shù)是一類非常重要的函數(shù) ， 特別是對(duì)金融數(shù)據(jù)處理和金融計(jì)算 ， 借助日期函數(shù)可以極大提高效率 。 ?如 1960年 1月 1日 9時(shí) 0分 ， 按日記的數(shù)值就是 0， 按小時(shí)記的數(shù)值就是 9。 run。 /*自由度為 2的 t分布的 */ q2=tinv(.95, 2, 3)。 /* 標(biāo)準(zhǔn)正態(tài)分布小于 */ put Y。 data。 目標(biāo)變量長(zhǎng)度不用缺省規(guī)則的函數(shù)見課本表 。 PUT函數(shù)是一個(gè)例外 ， 它不管自變量是什么類型 ，其結(jié)果總是字符值 。 例中 ， 兩個(gè) SUM函數(shù)的表示法是等價(jià)的 。 /*定義數(shù)組 y{10} */ x=sum (of y{*})。 例 數(shù)組元素自變量的縮寫方法。 自變量表示法： 當(dāng)函數(shù)有多個(gè)自變量時(shí) ， 必須用逗號(hào)分隔開 。 ， 自變量個(gè)數(shù)不能超過 32767。 變量名 X ， YX = R e pe a t( ‘ ‘ , 20 ) 。 例 函數(shù)及其自變量類型舉例。 例中 ， 兩段程序的結(jié)果相同 ， 第二段只有一個(gè)語句 ， 非常簡(jiǎn)單 。 使用 SAS函數(shù) 。 if totxy then least=totx。 不用 SAS函數(shù) 。 /* Where語句中使用 Year函數(shù) * run。 set 。 例 條件語句中使用函數(shù) 。 /*函數(shù) weekday給出變量 date在一周內(nèi)的哪一天 */ min=min(sum(oppr, hipr, lopr,clpr)， 1000)。 y=year(date)。 data a(keep=date y d min)。 第 3章 SAS函數(shù)與 CALL子程序 清華大學(xué)經(jīng)管學(xué)院 朱世武 Resdat樣本數(shù)據(jù)： SAS論壇： SAS函數(shù)定義 SAS函數(shù)是一個(gè)子程序 ， 對(duì)自變量返回一個(gè)結(jié)果值 。 ?比較規(guī)范的書寫有助于閱讀和檢查，還可以減少書寫錯(cuò)誤。 。 NOTE: “PROCEDURE PRINT”所用時(shí)間（總處理時(shí)間） : 實(shí)際時(shí)間 秒 CPU 時(shí)間 秒 日志窗口信息： SAS輸出 1SAS日志 2程序結(jié)果 3SAS 控制臺(tái)日志 SAS程序 SAS程序由多個(gè) SAS語句構(gòu)成。 18 run。 NOTE: “DATA 語句 ” 所用時(shí)間（總處理時(shí)間） : 實(shí)際時(shí)間 秒 CPU 時(shí)間 秒 16 。 NOTE: 對(duì) Age 無效的數(shù)據(jù)（在行 14 1314 中）。 11 input Name $ Age。 run。 Wang 35 Zhang A2 Hu 22 。 input Name $ Age。 數(shù)據(jù)