sas編程技術(shù)全局通用語句-資料下載頁
2025-08-11 14:17本頁面
【導(dǎo)讀】全局通用語句是可以用在任何地方的SAS語。說明，或者介紹整個程序的步驟或算法等。中間的信息可以包含分號，但這種形式的注釋不能嵌套。*對數(shù)據(jù)集means過程;例標(biāo)準(zhǔn)SAS程序開頭，記錄SAS程序信息的注釋形式。運(yùn)行SAS系統(tǒng)時，發(fā)布主機(jī)操作系統(tǒng)命令。每一個TITLE語規(guī)定一級標(biāo)題，最多可規(guī)定10級標(biāo)題。n緊跟在詞TIILE后面的數(shù)字，用來規(guī)定標(biāo)題的級別;例取消第3級及以后的所有標(biāo)題內(nèi)容。FOOTNOTE語句在每一頁的底部輸出一些腳注行。footnote‘清華大學(xué)金融系’;例取消所有已規(guī)定的腳注行。種情況下，第一個RUN語不能省略。項CANCEL結(jié)束當(dāng)前步的執(zhí)行。LIBNAME語句把一個libref和一個目錄名聯(lián)系起。三種格式反映了LIBNAME語句的三種用法。
【正文】 show names。 store a b。 show storage。 remove a b。 /*移走矩陣 a,b*/ show storage。 load a b。 quit。 proc iml。 A={1 2 3, 4 5 6, 7 8 9}。 B={2 2 2}。 print a b。 free a b。 /*清空矩陣 a, b*/ print a b。 quit。 缺失值 ? SAS/IML軟件也有缺失值的問題 。 注意 ， 矩陣中的缺失值不能視為 0， 矩陣標(biāo)識中 ， 數(shù)值缺失值用點 ” .” 表示 。 ? SAS/IML軟件支持缺失值的方式是有限的 。 許多矩陣運(yùn)算不允許有缺失值 ， 比如 ， 矩陣乘法 、 矩陣求逆等 。 另外 ，支持缺失值的矩陣運(yùn)算也會產(chǎn)生不一致的結(jié)果 。 線性回歸 IML模塊 解方程組 例 線性解方程組 。 用符號表示為： Ax = c 解方程組程序： proc iml。 a={3 1 2, 2 2 3, 4 1 4}。 c={8, 2, 9}。 x=inv(a)*c。 print a c x。 quit。 run。 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模塊程序 。 proc iml。 x={1 1 1, 1 2 4, 1 3 9, 1 4 16, 1 5 25}。 y={1,5,9,23,36}。 b=inv(t(x)*x)*t(x)*y。 yhat=x*b。 r=yyhat。 sse=ssq(r)。 dfe=nrow(x)ncol(x)。 mse=sse/dfe。 start regress。 /* 定義模塊開始 */ xpxi=inv(t(x)*x)。 /* 矩陣 X39。X的逆 */ beta=xpxi*(t(x)*y)。 /* 參數(shù)估計 */ yhat=x*beta。 /* 預(yù)測值 */ resid=yyhat。 /* 殘差 */ sse=ssq(resid)。 /* SSE */ n=nrow(x)。 /* 觀測值數(shù) */ dfe=nrow(x)ncol(x)。 /* 誤差自由度 DF */ mse=sse/dfe。 /* MSE */ cssy=ssq(ysum(y)/n)。 /* 校正平方和 */ rsquare=(cssysse)/cssy。 /* RSQUARE */ print,Regression Results, sse dfe mse rsquare。 stdb=sqrt(vecdiag(xpxi)*mse)。 /* 參數(shù)估計的標(biāo)準(zhǔn)差 */ t=beta/stdb。 /* 參數(shù)的 t檢驗 */ prob=1probf(tt,1,dfe)。 /* p值 */ print,Parameter Estimates, beta stdb t prob。 print,y yhat resid。 finish regress。 /* 模塊結(jié)束 */ reset noprint。 run regress。 /* 執(zhí)行 REGRESS 模塊 */ reset print。 /* 打開自動打印 */ covb=xpxi*mse。 /* 參數(shù)估計協(xié)方差陣 */ s=1/sqrt(vecdiag(covb))。 corrb=diag(s)*covb*diag(s)。 /* 參數(shù)估計的相關(guān)系數(shù) */ x1={1,2,3,4,5}。 /* 矩陣 X的第二列 */ x=orpol(x1,2)。 /* 產(chǎn)生正交多項式 */ reset noprint。 /* 關(guān)閉自動打印 */ run regress。 /* 運(yùn)行 REGRESS */ reset print。 covb=xpxi*mse。 s=1/sqrt(vecdiag(covb))。 corrb=diag(s)*covb*diag(s)。 reset fuzz。 corrb=diag(s)*covb*diag(s)。 矩陣操作 ? 輸入數(shù)據(jù)創(chuàng)建矩陣標(biāo)識 ? 使用賦值語句創(chuàng)建矩陣 ? 使用矩陣表達(dá)式 ? 利用行列標(biāo)展現(xiàn)矩陣 ? 缺失值運(yùn)算舉例 輸入數(shù)據(jù)創(chuàng)建矩陣標(biāo)識 SAS/IML中可以使用多種方式輸入數(shù)據(jù)創(chuàng)建矩陣標(biāo)識。 創(chuàng)建標(biāo)量矩陣標(biāo)識 例 創(chuàng)建標(biāo)量矩陣?yán)?。 proc iml。 x=12。 y=。 z=.。 a=39。Hello39。 b=Hi there。 print x y z a b。 X Y Z A B 12 . Hello Hi there 創(chuàng)建多元素矩陣 創(chuàng)建多元素矩陣時 ， 用括號把元素括起來 ， 用逗號 （ ， ） 把行分開 。 例 創(chuàng)建多元素矩陣?yán)?。 reset print。 coffee={4 2 2 3 2, 3 3 1 2 1, 2 1 0 2 1, 5 4 4 3 4}。 names={Jenny, Linda, Jim, Samuel}。 print coffee [rowname=names]。
畢業(yè)設(shè)計相關(guān)推薦
鄂ICP備17016276號-1