【正文】 do a=1 to 3。 means i/duncan alpha=.05。 。 end。 do l=1 to 5。 試對(duì)該組數(shù)據(jù)進(jìn)行方差分析 。 model y=blk trt。 cards。 run。 proc anova。 end。 do i=1 to 3。 [eganova2] 二因素有重復(fù)資料的方差分析 A因素有 a個(gè)水平 ，B因素有 b個(gè)水平 ， 共有 ab個(gè)處理組合 ， 每一組合有n個(gè)觀察值 ， 則該資料有 abn個(gè)觀察值 。 Class a b。 Cards。 Input y 。 [eganova12] 設(shè)有 A和 B兩個(gè)因素 ， A因素有 a個(gè)水平 ， B因素有 b個(gè)水平 ， 每一處理有單個(gè)觀察值 ， 全試驗(yàn)共有a*b個(gè)觀察值 ， 如下表： 兩向分組資料的方差分析 A因素 B因素 Ti. B1 B2 … Bb A1 Y11 Y12 … Y1b T1. A2 Y21 Y22 … Y2b T2. … … … … … … Aa Ya1 Ya2 … Yab Ta. … T 兩向分組資料方差分析表 22. //iT b T a b??22. //jT a T a b??22 /Y T a b??變異來(lái)源 df SS MS F A因素 a1 MSA MSA/MSe B因素 b1 MSB MSB/MSe 誤差 (a1)(b1) SSTSSASSB MSe 總變異 ab1 jH iM 1H 2H 3H 1M 13 14 14 2M 12 12 13 3M 3 3 3 4M 10 9 10 5M 2 5 4 將一種生長(zhǎng)激素配成 M1M5 5種濃度 ， 并用 H1H3三種時(shí)間浸漬某大豆品種的種子 ， 45天后得各處理每一植株的平均干物重 (克 )于下表 ， 試作方差分析 。 class i。 cards。 output。 石棉肺患者： 可疑患者： 非患者： GLM過(guò)程應(yīng)用實(shí)例 data new。 語(yǔ)句說(shuō)明： Means語(yǔ)句 計(jì)算平均數(shù)，并可選用多種方法進(jìn)行多重比較。 Run。 Class 變數(shù)表 。 [eganova1] GLM過(guò)程應(yīng)用范圍較廣 ， 除了用于方差分析外 ， 還可用于回歸分析 、 協(xié)方差分析等 。 Class i。 drop j。 Input y 。39。output。 兩因素帶互作模型。 單因素方差分析。 Model 線性數(shù)學(xué)模型 上述語(yǔ)句在于定義線性數(shù)學(xué)模型。 Run。 Class 變數(shù)表 。 概述 在 SAS中用于方差分析的主要過(guò)程有ANOVA（ 方差分析 ） 和 GLM（ 廣義線性模型 ） 。 方差分析的基本思想是將所有觀察值的總變異分解成不同變異來(lái)源的平方和和自由度 ， 進(jìn)而獲得不同變異來(lái)源的方差估計(jì) ， 以 F測(cè)驗(yàn)鑒定多個(gè)樣本平均數(shù)之間的差異顯著性 。 TTEST過(guò)程用于兩個(gè)獨(dú)立樣本的平均數(shù)差異顯著性測(cè)驗(yàn) 。 var y。 cards。 應(yīng)用舉例 3－ 兩個(gè)樣本平均數(shù)的組群比較 測(cè)定兩小麥品種千粒重 (克 )， A品種取樣 6次： ； B品種取樣 5次： 。 *如 data=new。 proc means n mean stderr t prt。 input a b 。 。 應(yīng)用舉例 1－ 單個(gè)平均數(shù)的測(cè)驗(yàn) 測(cè)定某稻田的地表光強(qiáng) 4次 ， 得結(jié)果為 , , , （ 千勒克斯 ） ， 試測(cè)驗(yàn)該結(jié)果與根據(jù) BeerLambert定律推出的理論值 μ0= ？ data new。 means過(guò)程 MEANS過(guò)程格式： PROC MEANS N MEAN STD STDERR T PRT。 5 統(tǒng)計(jì)推斷 統(tǒng)計(jì)推斷是統(tǒng)計(jì)分析的重要內(nèi)容 ， 用于顯著性測(cè)驗(yàn)的 SAS過(guò)程主要有 MEANS和TTEST。 clear。 語(yǔ)句行例： proc univariate data=aaa plot。 在默認(rèn)情況下 ， 即不使用 VAR語(yǔ)句指定分析變數(shù)時(shí) ， SUMMARY過(guò)程僅進(jìn)行對(duì)觀測(cè)值的計(jì)數(shù)工作 ，其他各種統(tǒng)計(jì)數(shù)的計(jì)算都將被忽略 。 MAXDEC＝ n， 設(shè)定結(jié)果輸出時(shí)的小數(shù)位數(shù) ， 缺省為 [eg9]。 指定統(tǒng)計(jì)量的輸出數(shù)據(jù)集名 關(guān)鍵字 =新變量名 。 指明分組變量 CLASS 變量表 。 這些過(guò)程主要有： means, summary, univariate。 若原數(shù)據(jù)未經(jīng)排序 ， 應(yīng)使用sort過(guò)程對(duì)相應(yīng)的分組變數(shù) （ x） 進(jìn)行排序 。 如上例中我們只想顯示大于 50的 x變數(shù)的值 ， 則在 print過(guò)程中加入 where語(yǔ)句如下： Where x50。 Data=數(shù)據(jù)集 該選項(xiàng)用于指明所需處理的數(shù)據(jù)集名 ， 無(wú)此選項(xiàng)時(shí)將處理最后打開(kāi)的數(shù)據(jù)集 。 Where 條件表達(dá)式 ...。 文件 導(dǎo)入數(shù)據(jù) 文件類型 （ 有些選項(xiàng) ） 導(dǎo)入至哪一個(gè)數(shù)據(jù)集 （ 臨時(shí) ？ 永久 ？ ） 應(yīng)用 。 Input x1x4 。 它們均可用于標(biāo)志數(shù)據(jù)塊的開(kāi)始 ， 隨后緊跟需要讀入的數(shù)據(jù) 。 Zhang Shan F 18 176 75 Li Si M 19 163 55 Wang Wu F 17 169 70 Zhao Liu M 18 173 72 。 cards。 基本格式為： input 變數(shù)名 變數(shù)類型起止列數(shù) ； 例如： Input x y z； *確定 x、 y、 z 三個(gè)變數(shù)； Input x1x10。 則系統(tǒng)邏輯庫(kù)中建立 mylib數(shù)據(jù)庫(kù) ， 其中有 mydata數(shù)據(jù)集 (e:\work\)，其中的數(shù)據(jù)可重復(fù)使用 。 此時(shí)在臨時(shí)庫(kù) work中生成了一個(gè)名為new的數(shù)據(jù)集 。 數(shù)據(jù)庫(kù)可以分為永久庫(kù)和臨時(shí)庫(kù)兩種。 有了 SAS軟件 ， 你不必清楚怎么做 ， 但應(yīng)清楚做什么 ！ 只需將數(shù)據(jù)按 SAS過(guò)程步的格式準(zhǔn)備好 ， 交由其完成相應(yīng)的統(tǒng)計(jì)分析 。 注釋段落：用字符組 “ /*”和 “ */”包括起來(lái)的任何字符內(nèi)容 ， 可占多行 。 … otherwise 執(zhí)行語(yǔ)句 Z。 Do until … end。 如： if … then statement 1。 class 變數(shù)； 等 。 PROC 過(guò)程名 DATA＝數(shù)據(jù)集名 選項(xiàng) 。 過(guò)程步 (Proc Step)將數(shù)據(jù)集完成相應(yīng)統(tǒng)計(jì)分析 。 命令欄 2 SAS程序語(yǔ)句基本結(jié)構(gòu)簡(jiǎn)介 數(shù)據(jù)步與過(guò)程步 SAS軟件就是將數(shù)據(jù)交給 SAS過(guò)程 (proc)計(jì)算分析并顯示結(jié)果的軟件 。 菜單和工具欄 工具欄圖標(biāo)提供了常見(jiàn)任務(wù)的快捷操作方式 新建 打開(kāi) 保存 打印 打印預(yù)覽 剪切 復(fù)制 粘貼 撤銷 建立新庫(kù) 資源瀏覽器窗口 執(zhí)行 清除 暫停 幫助 在命令欄可輸入一些命令 ， 可方便地實(shí)現(xiàn)某些功能 。 View(瀏覽 )－支持 SAS系統(tǒng)用戶在多窗口之間切換 。 6) 其它子窗口 除了增強(qiáng)型編輯器外 ， SAS還提供了普通的程序編輯 （ program editor） 窗口； 當(dāng)使用 SAS作圖時(shí) ， 相應(yīng)的統(tǒng)計(jì)圖會(huì)在專門(mén)的Graph窗口中輸出； SAS的數(shù)據(jù)集顯示窗口 Viewtable等 。 結(jié)果輸出窗口中的內(nèi)容可以保存為擴(kuò)展名為 .lst的純文本文件 。 窗口內(nèi)容可保存為擴(kuò)展名為 “ .log”的文本文件 。 1) 增強(qiáng)型程序編輯 (Enhanced Editor)窗口 Log窗口顯示程序運(yùn)行有關(guān)信息： 程序行 ， 黑色 ， 記錄執(zhí)行過(guò)的每一條語(yǔ)句 。 SAS的操作方式 SAS的工作界面 菜單欄 命令欄 狀態(tài)欄 工具欄 窗體 窗口欄 1 SAS系統(tǒng)的工作環(huán)境 SAS啟動(dòng)時(shí)，默認(rèn)會(huì)打開(kāi)以下五個(gè)窗口： 1）增強(qiáng)型程序編輯器 (Enhanced Editor)窗口 2）日志（ Log）窗口 3）結(jié)果輸出（ Output）窗口 4）資源管理器（ Explorer）窗口 5）結(jié)果（ Results）窗口 主要窗口 系統(tǒng)默認(rèn)提供的程序編輯窗口，編輯 SAS應(yīng)用程序。 SAS的基本部分是 SAS/BASE模塊 ， 該模塊是SAS系統(tǒng)的核心 ， 承擔(dān)著主要的數(shù)據(jù)管理任務(wù) ， 并管理 SAS的用戶使用環(huán)境 ， 進(jìn)行用戶語(yǔ)言的處理 ， 調(diào)用其它 SAS模塊和產(chǎn)品 。 有些機(jī)構(gòu)和雜志只認(rèn) SAS分析的結(jié)果 。 面向 Windows9 Windows NT 和 UNIX 的 Orlando（ SAS System ） 上市； 1996：支持 Web， JMP 軟件和 SAS 軟件通過(guò) Windows 95 標(biāo)志認(rèn)證； 1997： SAS的全球員工人數(shù)達(dá) 5,000人； S座辦公樓竣工； SAS 收購(gòu)生命科學(xué)統(tǒng)計(jì)軟件 Stat View； SAS 總裁 Jim Goodnight 榮獲 1997年 American Academy of Achievements頒發(fā)的榮譽(yù)； 《 Fortune》 雜志的 “ 美國(guó) 100強(qiáng)公司 ” 評(píng)比中排名第 3； SAS/Warehouse Administrator、 SAS/IntrNet 軟件投入生產(chǎn)； 1998： SAS推出 CRM、 SAS/ACCESS、 OLEDB for OLAP 和 SAS/Enterprise Reporter 和 HR Vision 上市； 1999： SAS收入超過(guò) 10億美元； T座辦公樓開(kāi)放； SAS支持報(bào)刊調(diào)查工作；美國(guó)食品和藥品管理局選擇 SAS技術(shù) ， 作為分析數(shù)據(jù)的標(biāo)準(zhǔn) 。 公司大事記 1985：與 HP合作； SAS/AF、 SAS/DMI和 （ SAS/IML） 軟件；首個(gè) PC DOS SAS System 版本 （ Base SAS 和 SAS/RTERM 軟件 ） 推出；PC的微型機(jī)－大型機(jī)鏈接 （ SAS/C 編譯器 ） 上市； 1986： SAS軟件研究所收購(gòu) Lattice有限公司；與微軟公司建立合作關(guān)系；SAS/QC 、 SAS/IML 和 SAS/STAT 軟件上市； 1987： J座辦公樓竣工；與蘋(píng)果電腦合作； SAS/SHARE、 DB2和 SQLDS；用 C語(yǔ)言重新編寫(xiě)了 SAS系統(tǒng)； ； 1988：?jiǎn)T工人數(shù)達(dá)千名；收購(gòu) NeoVisuals公司；與 Sun Microsystems 合作 。 經(jīng)過(guò) 30多年發(fā)展 ， SAS已被全世界 120多個(gè)國(guó)家和地區(qū)的科研機(jī)構(gòu)和人員普遍采用 ， 涉及教育 、 科研 、金融等各個(gè)領(lǐng)域 。 其創(chuàng)業(yè)產(chǎn)品 — 統(tǒng)計(jì)分析系統(tǒng)部分 ， 以應(yīng)有盡有 、 包羅萬(wàn)象和強(qiáng)大精準(zhǔn)的數(shù)據(jù)分析能力 ， 一直為業(yè)界推崇 ， 被視為最權(quán)威的統(tǒng)計(jì)分析標(biāo)準(zhǔn)軟件 。 發(fā)生發(fā)展（歷史） 1976： SAS軟件研究院 (北卡羅萊那州 )成立；推出首個(gè)產(chǎn)品 － Base SAS軟件； 1977： SAS由創(chuàng)業(yè)初期的 7名員工組成； 1978： SAS的客戶數(shù)達(dá)到了 600； 1979： SAS公司軟件第一次外銷； SAS軟件取得技術(shù)進(jìn)步 ， 可以運(yùn)行在IBM的 VM/CMS主機(jī)環(huán)境下； 1980： SAS歐洲總部在英國(guó)開(kāi)張；公司遷于美國(guó)北卡州卡莉市的 A座新址； SAS/GRAPH軟件和 SAS/ETS計(jì)量經(jīng)濟(jì)學(xué)軟件以及 SAS時(shí)間序列分析軟件的圖形展現(xiàn)模塊發(fā)布； 1981： SAS軟件銷售速度加快； 3000多個(gè)軟件安裝點(diǎn)上安裝 SAS 軟件； 1982：歐洲總部搬遷到海德堡；亞太 （ 新西蘭辦事處 、 澳大利亞分公司 ）成立；與 DEC公司建立合作關(guān)系； 1983： E座辦公樓竣工使用； SAS軟件全球銷量增加； 1984：項(xiàng)目管理軟件上市； SAS 收購(gòu) SYSTEM 2022 數(shù)據(jù)管理軟件； SAS 軟件擴(kuò)展 ， 從大型機(jī)到微型計(jì)算機(jī) 。 公司大事記（續(xù)） 1994： SAS全球員工人數(shù)達(dá) 3,000人； SAS/SPECTRAVIEW, SAS/SHARE*NET、JMP 上市； 1995： Rapid Warehousing Program； EIS市場(chǎng)領(lǐng)先； SAS for Macintosh（ 本 ） 上市 。 公司大事記（續(xù)） 1） 功能強(qiáng)大 ， 統(tǒng)計(jì)方法齊 、 新 、 優(yōu) SAS提供了從基本統(tǒng)計(jì)數(shù)計(jì)算到各種試驗(yàn)設(shè)計(jì)的方差分析 ， 相關(guān)回歸分析以及多變量分析的各種統(tǒng)計(jì)分析過(guò)程 ， 幾乎囊括了所有的最新統(tǒng)計(jì)分析方法 ，其分析技術(shù)先進(jìn) ， 可靠 。 3） 運(yùn)算速度快 ， 聯(lián)機(jī)幫助功能強(qiáng)大