【正文】 。 語句格式： OPTIONS option1 ...optionn。 例中 ， 運行調試好的大型程序時 ， 為了抑制 LOG和 OUTPUT窗口的輸出 ， 用上面的系統(tǒng)選項就可以解決問題 。 IML最基本的數(shù)據(jù)對象是矩陣 。 ? SAS/IML軟件是一種編程語言； ? SAS/IML軟件對矩陣作運算； ? SAS/IML軟件有強大的算符功能 ， 有許多函數(shù)和 CALL子程序； ? SAS/IML軟件將算符應用到整個矩陣； ? SAS/IML軟件是交互式的； ? SAS/IML軟件是動態(tài)的； ? SAS/IML軟件可以處理數(shù)據(jù)； ? SAS/IML軟件可以作圖 。 /*指定 ResDat為缺省的 SAS邏輯庫 */ start approx(x)。 y=.5(z+x/z)。 /* 返回近似值 */ finish approx。 run。 M 1矩陣稱作列向量 。 A, XX, MATRIX1. 例 。 a=100。 用逗號將矩陣的行分開 。 /*定義一個列向量 Y*/ Z={1 2, 3 4, 5 6}。 /*字符長度以較長的字符為準 ， 這里為 6*/ B={39。}。Yes39。 ANS={39。, 39。}。 r=rank(x)。 其中： Argument可以是矩陣名 ， 矩陣本身或表達式 。 Rv =Ranuni(Repeat(0,10,1))。 Rv =Rancau(Repeat(0,10,1))。 Rv =Rantri(Repeat(0,10,1),3)。 retain _seed_ 0。 Normal1 = mu + sigma * rannor(_seed_)。 run。 call eigen(val,vec,t)。 例 自建子程序 。 finish。 CALL的其它一用法： call gstart。 /*顯示散點圖 */ 命令 利用 SAS/IML命令用于實現(xiàn)一些具體的系統(tǒng)功能 ， 例如 ， 存貯和加載矩陣 、 或進行特定的數(shù)據(jù)處理等 。 B={2 2 2}。 remove a b。 proc iml。 free a b。 注意 ， 矩陣中的缺失值不能視為 0， 矩陣標識中 ， 數(shù)值缺失值用點 ” .” 表示 。 線性回歸 IML模塊 解方程組 例 線性解方程組 。 x=inv(a)*c。 1 2 31 2 31 2 33 2 82 2 3 24 4 9x x xx x xx x x? ? ??? ? ? ??? ? ? ?? A C X 3 1 2 8 3 2 2 3 2 5 4 1 4 9 2 結果： 線性回歸 IML模塊程序 例 IML模塊程序 。 b=inv(t(x)*x)*t(x)*y。 dfe=nrow(x)ncol(x)。 /* 矩陣 X39。 /* 殘差 */ sse=ssq(resid)。 /* MSE */ cssy=ssq(ysum(y)/n)。 /* 參數(shù)估計的標準差 */ t=beta/stdb。 finish regress。 /* 打開自動打印 */ covb=xpxi*mse。 /* 矩陣 X的第二列 */ x=orpol(x1,2)。 covb=xpxi*mse。 corrb=diag(s)*covb*diag(s)。 x=12。Hello39。 例 創(chuàng)建多元素矩陣例 。 print coffee [rowname=names]。 coffee={4 2 2 3 2, 3 3 1 2 1, 2 1 0 2 1, 5 4 4 3 4}。 print x y z a b。 z=.。 創(chuàng)建標量矩陣標識 例 創(chuàng)建標量矩陣例 。 corrb=diag(s)*covb*diag(s)。 /* 關閉自動打印 */ run regress。 corrb=diag(s)*covb*diag(s)。 run regress。 /* p值 */ print,Parameter Estimates, beta stdb t prob。 /* RSQUARE */ print,Regression Results, sse dfe mse rsquare。 /* 觀測值數(shù) */ dfe=nrow(x)ncol(x)。 /* 參數(shù)估計 */ yhat=x*beta。 start regress。 r=yyhat。 x={1 1 1, 1 2 4, 1 3 9, 1 4 16, 1 5 25}。 quit。 a={3 1 2, 2 2 3, 4 1 4}。 許多矩陣運算不允許有缺失值 ， 比如 ， 矩陣乘法 、 矩陣求逆等 。 quit。 B={2 2 2}。 load a b。 store a b。 proc iml。 /*打開一個新的 graphics片段 */ call gpoint(x,y)。 /*運行模型 MYMOD*/ run mymod(x, y, 10)。 a=sqrt(c)。 第一個例句 ， 用 CALL語句創(chuàng)建矩陣 , 得到矩陣 t的特征值和特征向量 。 在 SAS/IML中 x1x10被當成 x1減 x10， 而不是 x1到 x10. CALL語句與子程序 CALL SUBROUTINE arguments 。 end。 sigma=1。 例 DATA步下隨機數(shù)產(chǎn)生 SAS程序例子 。 Rv =Rangam(Repeat(0,10,1),1)。 Rv =Rannor(Repeat(0,10,1))。 函數(shù)