【正文】 list。 例 指定缺省邏輯庫(kù) 。 setin class point 10。 EDIT語(yǔ)句 EDTI語(yǔ)句用于讀寫(xiě) SAS數(shù)據(jù)集 。 pause *。 else do。 print y。 y=sqrt(x)。 end。 print count。 例中，將外部文件中的數(shù)據(jù)導(dǎo)入到一個(gè)向量。 例 。 z1=abs(x+y)。例 用 ANY函數(shù)的表達(dá)式。 else w=0。 /* 矩陣相加 */ C=XY。 names={Jenny, Linda, Jim, Samuel}。 T 3 2 0 6 0 3 2 2 2 簡(jiǎn)化下標(biāo)算符 簡(jiǎn)化下標(biāo)算符列表 例 簡(jiǎn)化下標(biāo)算符的應(yīng)用 。 t={ 3 2 1, 6 4 3, 2 2 2 }。 /*事先創(chuàng)建向量 cols*/ submat1=coffee[rows,cols]。 /*第 3行元素 */ print jim。 print coffee[rowname=names]。 a=(x)(yz)。 r = { 1 0. 5 0 0. 5 } 。 S = { 10 9 8 7 6} 。 aa=shape({99 33,99 33},3,3)。 DESIGNF( columnvector )。 c=block(a,b)。 grandtot=sum(coffee)。 print coffee [rowname=names]。Hello39。 corrb=diag(s)*covb*diag(s)。 /* 矩陣 X的第二列 */ x=orpol(x1,2)。 finish regress。 /* MSE */ cssy=ssq(ysum(y)/n)。 /* 矩陣 X39。 b=inv(t(x)*x)*t(x)*y。 x=inv(a)*c。 注意 ， 矩陣中的缺失值不能視為 0， 矩陣標(biāo)識(shí)中 ， 數(shù)值缺失值用點(diǎn) ” .”表示 。 proc iml。 B={2 2 2}。 CALL的其它一用法： call gstart。 例 自建子程序 。 run。 retain _seed_ 0。 Rv =Rancau(Repeat(0,10,1))。 其中： Argument可以是矩陣名 ， 矩陣本身或表達(dá)式 。}。 ANS={39。}。 /*定義一個(gè)列向量 Y*/ Z={1 2, 3 4, 5 6}。 a=100。 M 1矩陣稱(chēng)作列向量 。 /* 返回近似值 */ finish approx。 /*指定 ResDat為缺省的 SAS邏輯庫(kù) */ start approx(x)。 IML最基本的數(shù)據(jù)對(duì)象是矩陣 。 SAS程序中的語(yǔ)句可分為兩類(lèi)步驟： ? 數(shù)據(jù)步，用來(lái)處理輸入數(shù)據(jù) ； ? 過(guò)程步，通過(guò)一系列 SAS語(yǔ)句來(lái)調(diào)用和執(zhí)行程序，通常以一個(gè) SAS數(shù)據(jù)集作為輸入數(shù)據(jù)。 12 cards。 data A。 data A。 var data clpr。 data a input x y。 /*將 data錯(cuò)拼為 daat */ x=2。 data one。 data a (keep=date x)。/*自動(dòng)將數(shù)值變量 Y轉(zhuǎn)換為字符變量 */ put Z=。 NOTE: “DATA 語(yǔ)句”所用時(shí)間（總處理時(shí)間） : 實(shí)際時(shí)間 秒 CPU 時(shí)間 秒 data。 日志信息： 10 data。 code=input(stkcd, 6.)。 /*需要數(shù)值變量的函數(shù)中引用字符變量時(shí)， 自動(dòng)將字符變量 Y轉(zhuǎn)換為數(shù)值變量 */ run。 data。 /*X為數(shù)值變量 */ Y=?10?。 ?使用 FORMAT 或 INFORMAT語(yǔ)句 Data b。 創(chuàng)建變量 ?使用賦值語(yǔ)句 Data a。字符變量名后跟一個(gè)美元號(hào)（ $）表示該變量是字符型而不是數(shù)值型。?。 輸出結(jié)果： X= Y= Y=. 比較算符 例 比較算符應(yīng)用。 X=**。 */ 要注意的是，字符常數(shù)是由引號(hào)括起來(lái)的，但是字符變量的名稱(chēng)則沒(méi)有引號(hào)，也就是說(shuō)字符常數(shù)不能作為字符變量的名稱(chēng)。 數(shù)值常數(shù) 例 數(shù)值常數(shù)。/*取數(shù)據(jù)集前五個(gè)觀測(cè) */ Data one。 do I=1 to 100。 參數(shù) PARAMETER不是任選項(xiàng)，不用括號(hào)，OPTIONS是選項(xiàng)關(guān)鍵詞。 CALL子程序 用來(lái)改變變量值，或執(zhí)行一些系統(tǒng)功能。 data scores(keep=team game1 game2 game3)。 6 quit。 日志中信息 PROC SQL視圖 PROC SQL視圖為一個(gè)帶有名稱(chēng)的 PROC SQL查詢(xún)。 Data 。 ? 數(shù)據(jù)視圖可以保證讀取的數(shù)據(jù)集永遠(yuǎn)都是最新的。而 SAS視圖不能被索引。 一般來(lái)說(shuō)， SAS數(shù)據(jù)文件又分為： ? 普通的 SAS數(shù)據(jù)文件： SAS格式的數(shù)據(jù)文件。 庫(kù)引擎功能包括： ? 讀取和寫(xiě)入數(shù)據(jù)； ? 列出庫(kù)中的文件； ? 刪除和重命名文件。)。 建立 SAS邏輯庫(kù)的方法： ? 用菜單操作； ? 用 LIBNAME語(yǔ)句。外部文件一般用于儲(chǔ)存數(shù)據(jù)。 外部文件一般用來(lái)儲(chǔ)存： ? 要讀入 SAS數(shù)據(jù)文件的原始數(shù)據(jù) ? SAS程序語(yǔ)句 ? 過(guò)程步輸出 DBMS文件 SAS軟件可以和其它數(shù)據(jù)庫(kù)產(chǎn)品進(jìn)行數(shù)據(jù)文件轉(zhuǎn)換。 語(yǔ)句格式： LIBNAME libref engine 39。 例 引用非臨時(shí)庫(kù)的 SAS文件時(shí)必須使用兩級(jí)命名方式，而引用臨時(shí)庫(kù)的 SAS文件時(shí)，可以直接使用文件名，效果等同于。 SAS通過(guò)不同庫(kù)引擎讀寫(xiě)不同格式的文件。 ? 接口數(shù)據(jù)文件：以其它數(shù)據(jù)軟件格式儲(chǔ)存數(shù)據(jù)的數(shù)據(jù)文件。 ? …… 數(shù)據(jù)視圖 SAS數(shù)據(jù)視圖就是用來(lái)從其它文件中讀取數(shù)據(jù)的一種 SAS文件。 ? 更改一個(gè)數(shù)據(jù)視圖只需要改變這個(gè)視圖的查詢(xún)語(yǔ)句。 set 。它可以讀寫(xiě)的數(shù)據(jù)源同數(shù)據(jù)步視圖。 存儲(chǔ)編譯的 DATA步程序 語(yǔ)句格式 DATA datasetname(s) / PGM=storedprogramname (passwordoptionSOURCE=sourceoption)。 /*scores中只保留 team,game1game2 game3這 4個(gè)變量） proc print data=new(drop=year)。CALL子程序類(lèi)似于函數(shù)，但是不能在賦值語(yǔ)句中使用，所有的 SAS CALL子程序都只能被 CALL語(yǔ)句使用。 例 有效的 SAS語(yǔ)句。 x=y+1。 Set 。 1, –5, , , 2E4 , 20200 數(shù)值常數(shù)可以有多種格式展示。 日期時(shí)間常數(shù) 例 引用日期時(shí)間常數(shù)。 put X=。 if xy then c＝ 5； else c＝ 12。 put result。 變量屬性 SAS變量的屬性包括：長(zhǎng)度、輸人輸出格式和標(biāo)簽。 x=1。 Set b。 /*Y為字符變量 */ X=Y。 X=1。 函數(shù) INPUT將字符轉(zhuǎn)換為數(shù)值 INPUT函數(shù)的形式： INPUT(數(shù)據(jù)源 , 輸入格式 )。 run。 11 X=39。 Y=10。 run。 set 。 set 。 run。 z=x+y。 /* 將變量名 Date拼成 Data */ run。 input Stock $ Money Price。 input Name $ Age。 NOTE: 對(duì) Age 無(wú)效的數(shù)據(jù)（在行 14 1314 中）。 。 在 SAS/IML模式下 ， 可以立刻看到語(yǔ)句運(yùn)行的結(jié)果 ， 也可以把語(yǔ)句貯存在一個(gè)模塊中為以后使用 。 /* 定義模塊開(kāi)始 */ y=1。 /* 模塊結(jié)束 */ t=approx({3,5,7,9})。 1 1矩陣稱(chēng)作標(biāo)量 . 矩陣名字與標(biāo)識(shí) 矩陣名字是一個(gè)合法的 SAS名字 ， 不能超過(guò) 8字節(jié) 。 b=. 。 /*定義一個(gè) 3 2矩陣 Z*/ W=3Z。 重復(fù)因子 重復(fù)因子用括弧括起來(lái)放在要重復(fù)的元素前面 。Yes39。 賦值語(yǔ)句 結(jié)果＝表達(dá)式： w=3z。 SAS/IML軟件中的函數(shù)類(lèi)型： ? 矩陣查詢(xún)函數(shù)； ? 標(biāo)量函數(shù)； ? 概括函數(shù)； ? 矩陣變形函數(shù)； ? 矩陣算術(shù)函數(shù)； ? 線性代數(shù)和統(tǒng)計(jì)函數(shù) 。 Rv =Ranexp(Repeat(0,10,1))。 /*retain 賦初值 0給 _seed_*/ mu=0。 SAS/IML不支持 OF語(yǔ)句： A=mean(of x1x10)。 start mymod(a,b,c)。 /*激活 graphics模塊 */ call gopen。 show names。 A={1 2 3, 4 5 6, 7 8 9}。 ? SAS/IML軟件支持缺失值的方式是有限的 。 print a c x。 yhat=x*b。X的逆 */ beta=xpxi*(t(x)*y)。 /* 校正平方和 */ rsquare=(cssysse)/cssy。 /* 模塊結(jié)束 */ reset noprint。 /* 產(chǎn)生正交多項(xiàng)式 */ reset noprint。 矩陣操作 ? 輸入數(shù)據(jù)創(chuàng)建矩陣標(biāo)識(shí) ? 使用賦值語(yǔ)句創(chuàng)建矩陣 ? 使用矩陣表達(dá)式 ? 利用行列標(biāo)展現(xiàn)矩陣 ? 缺失值運(yùn)算舉例 輸入數(shù)據(jù)創(chuàng)建矩陣標(biāo)識(shí) SAS/IML中可以使用多種方式輸入數(shù)據(jù)創(chuàng)建矩陣標(biāo)識(shí)。 b=Hi there。 COFFEE 4 rows 5 cols (numeric) 4 2 2 3 2 3 3 1 2 1 2 1 0 2 1 5 4 4 3 4 NAMES 4 rows 1 col (character, size 6) JENNY LINDA JIM SAMUEL COFFEE JENNY 4 2 2 3 2 LINDA 3 3 1 2 1 JIM 2 1 0 2 1 SAMUEL 5 4 4 3 4 使用賦值語(yǔ)句創(chuàng)建矩陣 簡(jiǎn)單賦值語(yǔ)句 例 。 average=grandtot/ncol(coffee)。 例中 ， 創(chuàng)建分塊對(duì)角陣 C， 其主對(duì)