【正文】 將 date改名 */ 輸入和輸出格式 輸出格式 format是 SAS用來確定如何輸出或?qū)懗鰯?shù)據(jù)值的指令，用 FORMAT控制數(shù)據(jù)值在輸出時所使用的格式。所有的 SAS語句以分號結(jié)尾。 /* data為關(guān)鍵詞， data1為 SAS名 */ data _null_。 run。 表達(dá)式 表達(dá)式由一系列操作符和操作對象構(gòu)成，產(chǎn)生一個目標(biāo)值。 科學(xué)計數(shù)法： ， 。 ?01jan00?d。 put Y=。 比較準(zhǔn)則： ? 數(shù)值和字符都可以比； ? 結(jié)果為真賦值 1，假賦值 0； ? 字符值從左到右逐個按 ASCII碼排列序列進(jìn)行比較； ? 缺失值參加比較時，它比任何有效值都小。 SAS變量 SAS變量分為數(shù)值變量和字符變量。 變量的長度是指存貯變量值的字節(jié)數(shù)。/*數(shù)據(jù)集 a中有個變量 x，值為 1*/ ?使用 INPUT語句 Data b。 Run。 data。 /*Y為字符變量 */ if XY 。數(shù)據(jù)源既可以是數(shù)值型，也可以是字符型。 data。 /*X為字符變量 */ 12 Y=10。 /*自動將數(shù)值變量 Y轉(zhuǎn)換為字符變量 */ put Z=。 函數(shù) PUT將數(shù)值轉(zhuǎn)換為字符 PUT函數(shù)的形式： PUT(數(shù)據(jù)源 , 輸出格式 )。 run。 /* 只保留數(shù)據(jù)集中的數(shù)值變量 */ run。 /*將 data錯拼為 daat */ 14 WARNING 14169: 假定符號 DATA 錯拼為 daat。 cards。 NOTE: 正在寫入 HTML Body（主體）文件 : 8 var data clpr。 cards。 Wang 35 Zhang A2 Hu 22 。 NOTE: “DATA 語句 ” 所用時間（總處理時間） : 實際時間 秒 CPU 時間 秒 16 。 ?比較規(guī)范的書寫有助于閱讀和檢查，還可以減少書寫錯誤。因為 SAS/IML軟件是 SAS系統(tǒng)的一個模塊 ， 所以 ， 可以使用 SAS數(shù)據(jù)集和外部文件 ， 也可以對一個已有的 SAS數(shù)據(jù)集進(jìn)行編輯或者創(chuàng)建一個新的 SAS數(shù)據(jù)集 。 /* 循環(huán)開始 , le3= */ z=y。 /* 打印矩陣 */ quit。 例 矩陣名舉例 。 創(chuàng)建數(shù)值矩陣 含有多個元素的矩陣標(biāo)識要用括號把元素括起來 。 a={IML Modual}。 ANS={[2] 39。Yes39。 y=inv(x)。 Rv=Uniform(Repeat(0, 10, 1))。 Rv =Ranpoi(Repeat(0,10,1),2)。 do _i_ = 1 to 1000。 例 使用 CALL子程序 。 b=log(c)。 /*作散點圖 */ call gshow。 show storage。 print a b。 另外 ，支持缺失值的矩陣運算也會產(chǎn)生不一致的結(jié)果 。 run。 sse=ssq(r)。 /* 預(yù)測值 */ resid=yyhat。 stdb=sqrt(vecdiag(xpxi)*mse)。 /* 執(zhí)行 REGRESS 模塊 */ reset print。 /* 運行 REGRESS */ reset print。 proc iml。 X Y Z A B 12 . Hello Hi there 創(chuàng)建多元素矩陣 創(chuàng)建多元素矩陣時 ， 用括號把元素括起來 ， 用逗號 （ ， ） 把行分開 。 daycost=coffee。 創(chuàng)建特殊矩陣函數(shù) SAS/IML軟件有很多內(nèi)置函數(shù)來產(chǎn)生一些特殊矩陣。 one=j(2,5, 88)。 即 D={1 0, 1 0, 1 0, 0 1, 0 1, 1 1, 1 1} 例 由輸入矩陣創(chuàng)建一定形狀的矩陣。 如果步長不是 1， 可以用 DO函數(shù) 。 : 39。表達(dá)式可以出現(xiàn)在賦值語句的任何一邊。 二目元素算符 二目元素算符是對兩個矩陣的元素分別進(jìn)行操作。 也可以尋找超出行范圍的元素 ， 這時得到的是以行為順序第 6個元素 。 例 選擇矩陣的子陣 。 例 根據(jù)下標(biāo)位置進(jìn)行相應(yīng)元素的賦值 。 print i。 print A。 print coffee。 ELSE statement2 。 如果想要在 A和 B的對應(yīng)元素有一個不相等時為真 ， 要用 ANY函數(shù) 。 DO組語句 DO組的兩個主要用法： ? 將一組語句設(shè)成一個 DO組 ， 應(yīng)用時作為一個單位執(zhí)行； ? 將條件語句中的一組語句設(shè)成一個 DO組 。 else do。D:\ResDat\39。 例 循環(huán) DO例句 。 DO UNTIL語句 語句格式： DO UNTIL expression。 GOTO和 LINK語句都可以使得 IML立刻轉(zhuǎn)移到標(biāo)記語句 。 negative: print Sorry, X is negative。 x=abs(x)。 end。 pause msg。 CREATE SASdataset VAR operand 。 use 。 確定觀測范圍 確定觀測值范圍 range的關(guān)鍵字 Operand舉例列表 關(guān)鍵字 觀測范圍A L L 所有觀測C U R R E N T 當(dāng)前觀測N E X T numbe r 下一個或下 numbe r 個觀測A F T E R 當(dāng)前觀測后的所有觀測P O I N T ope r a nd 由 ope r a nd 指定觀測的記錄號O pe r a nd 舉例單個記錄號 p o i n t 5包括幾個記錄號的標(biāo)識 p o i n t { 2 5 1 0 }含有記錄號的矩陣名 p o i n t p括號中的表達(dá)式 p o i n t ( p + 1 )例 確定觀測范圍例句 。 p={2 4 9}。 reset deflib=ResDat。 /*顯示所有打開的 SAS數(shù)據(jù)集及其相應(yīng)的狀態(tài) */ show contents。 激活一個 SAS數(shù)據(jù)集 常用激活語句及功能 對一個已經(jīng)打開的數(shù)據(jù)集使用 USE， EDIT或 CREATE語句 ， 那么這個數(shù)據(jù)集變成了當(dāng)前數(shù)據(jù)集 。 ABORT語句 ABORT 。 PAUSE語句 PAUSE message * 。 end。 程序段二： do。 例 使用 GOTO或 LINK語句 。 do until(count5)。 DO While語句 語句格式： DO WHILE expression。 /* 開始讀循環(huán) */ input xx。 z2=abs(x+y)。 additional statements END 。 I f a l l ( x y ) th e n s ta te m e n t。 if max(a)20 then p=0。 X={1 2 ., . 5 6, 7 . 9}。 例 使用自動命名選項。 例如表達(dá)式可以為： t[loc(t0)]=0。 /*將矩陣 coffee的第 1行 2列處的元素值改為 4*/ print coffee。 /*第一個向量 ｛ 1 3｝ 選擇行 ， 第二個向量 ｛ 2 3 5｝ 選擇列 */ print submat1。 print c6。 邏輯 ” 和 ” 減 小于 元素相乘 = 小于等于 元素乘方 大于/ 除 = 大于等于 元素最大值 ^= 不等于 元素最小值 = 等于| 邏輯 ” 或 ” M O D ( m , n) 余數(shù)下標(biāo) 下標(biāo)的一般形式： operand[row,column] 其中： ? operand通常是一個矩陣名 ， 也可以是表達(dá)式或矩陣本身； ? row規(guī)定矩陣的行數(shù)，可以是表達(dá)式，也可以是標(biāo)量或向量； ? Column規(guī)定矩陣的列數(shù) ， 可以是表達(dá)式 ， 也可以是標(biāo)量或向量 。 算符優(yōu)先級 優(yōu)先級 算符^ 下標(biāo) ( 前綴 ) **II * / III + IV | | / / :V = = = ^ =VI amp。 。 例句 產(chǎn)生向量R = 1: 5。 ?其中， nrow, ncol, 和 padvalue是隨意的，但一般情況下要給他們賦值。 I( dimension )。 BLOCK( matrix1,matrix2,...,matrix15)。 ones={1,1,1,1,1}。 reset print。 y=。 s=1/sqrt(vecdiag(covb))。 /* 參數(shù)估計協(xié)方差陣 */ s=1/sqrt(vecdiag(covb))。 /* 參數(shù)的 t檢驗 */ prob=1probf(tt,1,dfe)。 /* SSE */ n=nrow(x)。 mse=sse/dfe。 proc iml。 用符號表示為： Ax = c 解方程組程序： proc iml。 /*清空矩陣 a, b*/ print a b。 /*移走矩陣 a,b*/ show storage。 例 SAS/IML命令用法舉例 。 run mymod(a, b, 10)。 call delete(DAT)。 /*均值為 mu,標(biāo)準(zhǔn)差為 sigma*/ output。 例中 ， 形式 REPEAT(X, Y, Z)中 ， X是隨機(jī)數(shù)種子 ， Y是產(chǎn)生的隨機(jī)數(shù)個數(shù) ， Z是產(chǎn)生的隨機(jī)數(shù)列數(shù) 。 Rv =Normal(Repeat(0,10,1))。 表達(dá)式中的三種運算符 表達(dá)式中的算符 算符 在表達(dá)式中的位置前綴運算符 放在運算對象前面 ( A )中綴運算符 放在運算對象之間 ( A * B )后綴運算符 放在運算對象后面 ( A 39。No39。, [2] 39。abc39。 例 由矩陣標(biāo)識創(chuàng)建的數(shù)值陣 。 100, ., HI, {1 2 2 8 8 }, {3, 3, 6,6}. 由矩陣標(biāo)識創(chuàng)建矩陣 創(chuàng)建標(biāo)量矩陣 例 。 T 3 OUTPUT窗口結(jié)果 牛頓迭代法求 f(z)=0一個實根的方法 牛頓迭代法公式： 1()()nnnnfzzzfz??? ?上例中函數(shù)為： xzzf ?? 2)(理解 IML語言 定義矩陣 SAS/IM軟件對矩陣作運算 ， 這里矩陣的定義同矩陣代數(shù)中矩陣的定義一樣 ， 是一個兩維的 （ 行 列 ） 數(shù)組 。 /* 估計開方根 ， 牛頓迭代 */ w=abs(yz)。 一個簡單 IML交互程序 例 求平方根 IML程序 。 第 18章 IML編程技術(shù) 清華大學(xué)經(jīng)管學(xué)院 朱世武 Resdat樣本數(shù)據(jù)： SAS論壇： 本章內(nèi)容包括： ? 概述； ? 理解 IML語言； ? 線型回歸 IML模塊 。 18 run。 run。 run。 9 run。 run。 6