【正文】 Means a b/duncan。 Class a b。 proc print。 Cards。 End。 Input y 。 Do a=1 to 5。 [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??變異來源 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 將一種生長激素配成 M1M5 5種濃度 ， 并用 H1H3三種時(shí)間浸漬某大豆品種的種子 ， 45天后得各處理每一植株的平均干物重 (克 )于下表 ， 試作方差分析 。 means i/lsd alpha=.05。 class i。 proc print。 cards。 end。 output。 do j=1 to 11。 石棉肺患者： 可疑患者： 非患者： GLM過程應(yīng)用實(shí)例 data new。 在 SAS中，數(shù)據(jù)行中的小數(shù)點(diǎn)“ .”代表缺失數(shù)據(jù)。 語句說明： Means語句 計(jì)算平均數(shù)，并可選用多種方法進(jìn)行多重比較。 Class語句 指定分類變數(shù)，此語句必須設(shè)定，并且應(yīng)出現(xiàn)在 model語句之前。 Run。 Means 效應(yīng)表 /選項(xiàng) 。 Class 變數(shù)表 。 與 ANOVA過程相比 ， GLM過程運(yùn)行時(shí)一般占用較多的計(jì)算機(jī)內(nèi)存和時(shí)間 。 [eganova1] GLM過程應(yīng)用范圍較廣 ， 除了用于方差分析外 ， 還可用于回歸分析 、 協(xié)方差分析等 。 Means i/lsd alpha=.05。 Class i。 24 30 28 26 27 24 21 26 31 28 25 30 32 33 33 28 21 22 16 21 。 drop j。 end。 Input y 。 do i=1 to 5。39。 log。output。 模型定義中可用 |和 n簡化表達(dá)方法 水稻施肥試驗(yàn) ， 5個(gè)處理 ABCDE分別代表氨水 、廢氨水 、 碳銨 、 尿素和對照 （ 不施肥 ） ， 每處理 4盆（ 完全隨機(jī)試驗(yàn) ） ， 產(chǎn)量 （ 克 /盆 ） 結(jié)果如下： A: 24 30 28 26 B: 27 24 21 26 C: 31 28 25 30 D: 32 33 33 28 E: 21 22 16 21 請作方差分析。 兩因素帶互作模型。 兩因素主效模型。 單因素方差分析。因而應(yīng)依據(jù)試驗(yàn)設(shè)計(jì)，給出正確的線性數(shù)學(xué)模型。 Model 線性數(shù)學(xué)模型 上述語句在于定義線性數(shù)學(xué)模型。 多重比較方法選擇有 T或 LSD（ 配對 t 測驗(yàn)或 Fisher氏最小顯著差數(shù)法 ） ， DUNCAN（ Duncan氏新復(fù)極差測驗(yàn)法 ） ， TUCKY（ Tukey氏固定極差測驗(yàn)法 ） 和 DUNNETT（ Dunt氏新復(fù)極差測驗(yàn)法 ） 等 。 Run。 Means 效應(yīng)表 /選項(xiàng) 。 Class 變數(shù)表 。 對于非平衡資料 ， 應(yīng)采用 GLM過程 。 概述 在 SAS中用于方差分析的主要過程有ANOVA（ 方差分析 ） 和 GLM（ 廣義線性模型 ） 。 當(dāng)處理效應(yīng)為隨機(jī)模型時(shí) ， 估計(jì)方差分量 。 方差分析的基本思想是將所有觀察值的總變異分解成不同變異來源的平方和和自由度 ， 進(jìn)而獲得不同變異來源的方差估計(jì) ， 以 F測驗(yàn)鑒定多個(gè)樣本平均數(shù)之間的差異顯著性 。 推斷過程的比較和選擇 6 方差分析 兩個(gè)平均數(shù)的假設(shè)測驗(yàn) ， 一般可以采用 t 測驗(yàn)的方法完成 。 TTEST過程用于兩個(gè)獨(dú)立樣本的平均數(shù)差異顯著性測驗(yàn) 。 [eg13] NEANS過程可用于推斷單個(gè)樣本是否抽自平均數(shù)為 0的總體 。 var y。 proc ttest。 cards。 data new。 應(yīng)用舉例 3－ 兩個(gè)樣本平均數(shù)的組群比較 測定兩小麥品種千粒重 (克 )， A品種取樣 6次： ； B品種取樣 5次： 。 *指明分類變量，這里分類變量只能有兩個(gè)水平，該語句是必需的； VAR 變量 。 *如 data=new。 run。 proc means n mean stderr t prt。 cards。 input a b 。 run。 。 y=y3。 應(yīng)用舉例 1－ 單個(gè)平均數(shù)的測驗(yàn) 測定某稻田的地表光強(qiáng) 4次 ， 得結(jié)果為 , , , （ 千勒克斯 ） ， 試測驗(yàn)該結(jié)果與根據(jù) BeerLambert定律推出的理論值 μ0= ？ data new。 結(jié)果輸出包括測驗(yàn)平均數(shù)與 0差異顯著性的 t 值和概率 p。 means過程 MEANS過程格式： PROC MEANS N MEAN STD STDERR T PRT。 在成對比較時(shí) ， 一般先計(jì)算對比差值 d，然后利用 MEANS過程進(jìn)行分析 。 5 統(tǒng)計(jì)推斷 統(tǒng)計(jì)推斷是統(tǒng)計(jì)分析的重要內(nèi)容 ， 用于顯著性測驗(yàn)的 SAS過程主要有 MEANS和TTEST。 clear。 clear。 在運(yùn)行程序后 ， 輸出 、 日志等窗口的內(nèi)容將在其中累積 ， 有時(shí)影響相關(guān)信息的查看 ， 有必要將其清除 。 語句行例： proc univariate data=aaa plot。 Univariate 過程的重點(diǎn)在于描述變量的分布 。 在默認(rèn)情況下 ， 即不使用 VAR語句指定分析變數(shù)時(shí) ， SUMMARY過程僅進(jìn)行對觀測值的計(jì)數(shù)工作 ，其他各種統(tǒng)計(jì)數(shù)的計(jì)算都將被忽略 。 Output out=數(shù)據(jù)集名 （ 若無指定數(shù)據(jù)集 ， 自動(dòng)以 data*命名 ） ， 其后可將下列之一統(tǒng)計(jì)數(shù)給以新的變數(shù)命名 ， 如 CSS=aa。 MAXDEC＝ n， 設(shè)定結(jié)果輸出時(shí)的小數(shù)位數(shù) ， 缺省為 [eg9]。 未指定時(shí) ， 輸出的統(tǒng)計(jì)數(shù)只有 N, MEAN, STD, MAX和 MIN。 指定統(tǒng)計(jì)量的輸出數(shù)據(jù)集名 關(guān)鍵字 =新變量名 。 指明頻數(shù)變量 WEIGHT 變量表 。 指明分組變量 CLASS 變量表 。 指定要分析的數(shù)據(jù)集名及 一些選項(xiàng) VAR 變量表 。 這些過程主要有： means, summary, univariate。 By x。 若原數(shù)據(jù)未經(jīng)排序 ， 應(yīng)使用sort過程對相應(yīng)的分組變數(shù) （ x） 進(jìn)行排序 。 By語句用于指定分組變數(shù) ， 如果你需要分組處理數(shù)據(jù) ， 例如要按 （ x） 類別分組輸出統(tǒng)計(jì)結(jié)果 ， 采用 by語句即可 。 如上例中我們只想顯示大于 50的 x變數(shù)的值 ， 則在 print過程中加入 where語句如下： Where x50。 如在數(shù)據(jù)集中有 x、 y、 z三個(gè)變數(shù) ， 若只想對個(gè)別變數(shù)統(tǒng)計(jì)處理 ， 則 print過程如下： Proc print； *對變數(shù)列出顯示； var x； *僅對 x變數(shù)進(jìn)行處理； Var x, y。 Data=數(shù)據(jù)集 該選項(xiàng)用于指明所需處理的數(shù)據(jù)集名 ， 無此選項(xiàng)時(shí)將處理最后打開的數(shù)據(jù)集 。 Run。 Where 條件表達(dá)式 ...。 導(dǎo)入數(shù)據(jù) 4 SAS過程步之描述性統(tǒng)計(jì) SAS過程步基礎(chǔ) PROC 過程名 data=數(shù)據(jù)集名 其他選項(xiàng) 。 文件 導(dǎo)入數(shù)據(jù) 文件類型 （ 有些選項(xiàng) ） 導(dǎo)入至哪一個(gè)數(shù)據(jù)集 （ 臨時(shí) ？ 永久 ？ ） 應(yīng)用 。 幾點(diǎn)注意： 1） 該命令只能讀入 txt文本或 dbf文件； 2） 文本文件必須是純文本格式 （ 不帶逗號(hào) 、 不帶制表符或 TAB符 ） ； 3） 也不宜有表頭等格式 。 Input x1x4 。 Cards/datalines語句 Data aaa。 它們均可用于標(biāo)志數(shù)據(jù)塊的開始 ， 隨后緊跟需要讀入的數(shù)據(jù) 。 Run。 Zhang Shan F 18 176 75 Li Si M 19 163 55 Wang Wu F 17 169 70 Zhao Liu M 18 173 72 。 Input name $110 sex $ age height weight。 cards。 *數(shù)據(jù)集為 ； input x y 。 基本格式為： input 變數(shù)名 變數(shù)類型起止列數(shù) ； 例如： Input x y z； *確定 x、 y、 z 三個(gè)變數(shù)； Input x1x10。 如：data=。 則系統(tǒng)邏輯庫中建立 mylib數(shù)據(jù)庫 ， 其中有 mydata數(shù)據(jù)集 (e:\work\)，其中的數(shù)據(jù)可重復(fù)使用 。 data 語句 若 libname mylib ‘e:\work’。 此時(shí)在臨時(shí)庫 work中生成了一個(gè)名為new的數(shù)據(jù)集 。 永久庫通過 Libname語句由用戶自行定義 ，永久庫的數(shù)據(jù)文件將被保留 ， 可以重復(fù)使用 。 數(shù)據(jù)庫可以分為永久庫和臨時(shí)庫兩種。 對于復(fù)雜的數(shù)據(jù)準(zhǔn)備 ， 擬在稍后再作交代 ，本節(jié)只討論相對簡單的數(shù)據(jù)準(zhǔn)備 。 有了 SAS軟件 ， 你不必清楚怎么做 ， 但應(yīng)清楚做什么 ！ 只需將數(shù)據(jù)按 SAS過程步的格式準(zhǔn)備好 ， 交由其完成相應(yīng)的統(tǒng)計(jì)分析 。 注釋 3 SAS的數(shù)據(jù) (data)步 通俗的講 ， SAS軟件就是將數(shù)據(jù)交給 SAS過程（ proc） 計(jì)算分析并顯示結(jié)果的軟件 。 注釋段落：用字符組 “ /*”和 “ */”包括起來的任何字符內(nèi)容 ， 可占多行 。 其基本格式為： 注釋語句：使用星號(hào) “ *” 開始 ， 可占多行 ， 以分號(hào) “ 。 … otherwise 執(zhí)行語句 Z。 when (條件 1) 執(zhí)行語句 A。 Do until … end。 Do … end。 如： if … then statement 1。 注：尖括號(hào)內(nèi)的選項(xiàng)可以省略 ， 過程將按最通常的情況來處理 ， 即：處理最新建立的 SAS數(shù)據(jù)集及其中所有的變數(shù) 。 class 變數(shù)； 等 。 WHERE 條件表達(dá)式 。 PROC 過程名 DATA＝數(shù)據(jù)集名 選項(xiàng) 。 最終以分號(hào) “ 。 過程步 (Proc Step)將數(shù)據(jù)集完成相應(yīng)統(tǒng)計(jì)分析 。 用 SAS軟件 ， 你不必清楚怎么做 ， 但需清楚做什么 ！ 只需將數(shù)據(jù)按一定格式準(zhǔn)備 ， 交由過程步即可實(shí)現(xiàn)相應(yīng)的統(tǒng)計(jì)分析 。 命令欄 2 SAS程序語句基本結(jié)構(gòu)簡介 數(shù)據(jù)步與過程步 SAS軟件就是將數(shù)據(jù)交給 SAS過程 (proc)計(jì)算分析并顯示結(jié)果的軟件 。 在命令欄內(nèi)還可輸入 help， 獲得即時(shí)幫助 。 菜單和工具欄 工具欄圖標(biāo)提供了常見任務(wù)的快捷操作方式 新建 打開 保存 打印 打印預(yù)覽 剪切 復(fù)制 粘貼 撤銷 建立新庫 資源瀏覽器窗口 執(zhí)行 清除 暫停 幫助 在命令欄可輸入一些命令 ， 可方便地實(shí)現(xiàn)某些功能 。 Run(運(yùn)行 )－用于執(zhí)行或調(diào)用程序。 View(瀏覽 )－支持 SAS系統(tǒng)用戶在多窗口之間切換 。 其中的主要命令及功能： File(文件 )－支持 SAS文件的調(diào)入 、 保存及打印 。 6) 其它子窗口 除了增強(qiáng)型編輯器外 ， SAS還提供了普通的程序編輯 （ program editor） 窗口； 當(dāng)使用 SAS作圖時(shí) ， 相應(yīng)的統(tǒng)計(jì)圖會(huì)在專門的Graph窗口中輸出； SAS的數(shù)據(jù)集顯示窗口 Viewtable等 。 在該窗口中將 SAS系統(tǒng)的所有輸出結(jié)果依次按照目錄樹的結(jié)構(gòu)加以排列 ， 每一個(gè)過程步的結(jié)果被表示為一個(gè)節(jié)點(diǎn) ， 展開該節(jié)點(diǎn)就可進(jìn)一步看到表示不同輸出內(nèi)容的子節(jié)點(diǎn) ， 使用鼠標(biāo)右鍵單擊每個(gè)節(jié)點(diǎn) ， 就可對輸出結(jié)果的各部分進(jìn)行察看 、 存儲(chǔ) 、 打印 、 刪除等操作 。 結(jié)果輸出窗口中的內(nèi)容可以保存為擴(kuò)展名為 .lst的純文本文件 。 結(jié)果輸出窗口中的內(nèi)容是分頁顯示的 ， 頁上方