sas編程技術(shù)全局通用語(yǔ)句(參考版)
2024-08-24 14:17本頁(yè)面
【正文】 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}。 例 創(chuàng)建多元素矩陣?yán)?。 print x y z a b。Hello39。 z=.。 x=12。 創(chuàng)建標(biāo)量矩陣標(biāo)識(shí) 例 創(chuàng)建標(biāo)量矩陣?yán)?。 corrb=diag(s)*covb*diag(s)。 corrb=diag(s)*covb*diag(s)。 covb=xpxi*mse。 /* 關(guān)閉自動(dòng)打印 */ run regress。 /* 矩陣 X的第二列 */ x=orpol(x1,2)。 corrb=diag(s)*covb*diag(s)。 /* 打開(kāi)自動(dòng)打印 */ covb=xpxi*mse。 run regress。 finish regress。 /* p值 */ print,Parameter Estimates, beta stdb t prob。 /* 參數(shù)估計(jì)的標(biāo)準(zhǔn)差 */ t=beta/stdb。 /* RSQUARE */ print,Regression Results, sse dfe mse rsquare。 /* MSE */ cssy=ssq(ysum(y)/n)。 /* 觀測(cè)值數(shù) */ dfe=nrow(x)ncol(x)。 /* 殘差 */ sse=ssq(resid)。 /* 參數(shù)估計(jì) */ yhat=x*beta。 /* 矩陣 X39。 start regress。 dfe=nrow(x)ncol(x)。 r=yyhat。 b=inv(t(x)*x)*t(x)*y。 x={1 1 1, 1 2 4, 1 3 9, 1 4 16, 1 5 25}。 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 結(jié)果： 線性回歸 IML模塊程序 例 IML模塊程序 。 quit。 x=inv(a)*c。 a={3 1 2, 2 2 3, 4 1 4}。 線性回歸 IML模塊 解方程組 例 線性解方程組 。 許多矩陣運(yùn)算不允許有缺失值 ， 比如 ， 矩陣乘法 、 矩陣求逆等 。 注意 ， 矩陣中的缺失值不能視為 0， 矩陣標(biāo)識(shí)中 ， 數(shù)值缺失值用點(diǎn) ” .” 表示 。 quit。 free a b。 B={2 2 2}。 proc iml。 load a b。 remove a b。 store a b。 B={2 2 2}。 proc iml。 /*顯示散點(diǎn)圖 */ 命令 利用 SAS/IML命令用于實(shí)現(xiàn)一些具體的系統(tǒng)功能 ， 例如 ， 存貯和加載矩陣 、 或進(jìn)行特定的數(shù)據(jù)處理等 。 /*打開(kāi)一個(gè)新的 graphics片段 */ call gpoint(x,y)。 CALL的其它一用法： call gstart。 /*運(yùn)行模型 MYMOD*/ run mymod(x, y, 10)。 finish。 a=sqrt(c)。 例 自建子程序 。 第一個(gè)例句 ， 用 CALL語(yǔ)句創(chuàng)建矩陣 , 得到矩陣 t的特征值和特征向量 。 call eigen(val,vec,t)。 在 SAS/IML中 x1x10被當(dāng)成 x1減 x10， 而不是 x1到 x10. CALL語(yǔ)句與子程序 CALL SUBROUTINE arguments 。 run。 end。 Normal1 = mu + sigma * rannor(_seed_)。 sigma=1。 retain _seed_ 0。 例 DATA步下隨機(jī)數(shù)產(chǎn)生 SAS程序例子 。 Rv =Rantri(Repeat(0,10,1),3)。 Rv =Rangam(Repeat(0,10,1),1)。 Rv =Rancau(Repeat(0,10,1))。 Rv =Rannor(Repeat(0,10,1))。 Rv =Ranuni(Repeat(0,10,1))。 函數(shù) 例 IML隨機(jī)數(shù)函數(shù)應(yīng)用舉例 。 其中： Argument可以是矩陣名 ， 矩陣本身或表達(dá)式 。 )語(yǔ)句類型 SAS/IML軟件語(yǔ)句大致分為三類： ? 控制語(yǔ)句； ? 函數(shù)和 CALL語(yǔ)句； ? 命令語(yǔ)句 。 r=rank(x)。 函數(shù)作為表達(dá)式： a=sqrt(b)。}。 39。, 39。 39。 ANS={39。No39。Yes39。 例 用與不用重復(fù)因子的等價(jià)語(yǔ)句 。}。 39。 /*字符長(zhǎng)度以較長(zhǎng)的字符為準(zhǔn) ， 這里為 6*/ B={39。 /*創(chuàng)建一個(gè)新矩陣 W， 它的每個(gè)元素是矩陣 Z相應(yīng)元素的 3倍 */ 創(chuàng)建字符矩陣 例 ， 字符長(zhǎng)度取最長(zhǎng)字符的長(zhǎng)度 。 /*定義一個(gè)列向量 Y*/ Z={1 2, 3 4, 5 6}。 X={1 2 3 4 5 6}。 用逗號(hào)將矩陣的行分開(kāi) 。 c=IML。 a=100。 proc iml。 A, XX, MATRIX1. 例 。 矩陣標(biāo)識(shí)是由它的值表示的矩陣 。 M 1矩陣稱作列向量 。 SAS/IM軟件可以定義數(shù)值矩陣也可以定義字符矩陣 。 run。 /* 調(diào)用函數(shù) APPROX */ print t。 /* 返回近似值 */ finish approx。 /* 計(jì)算估計(jì)值之間的差異 */ end。
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1