方差分析SAS課件

方差分析主要過(guò)程說(shuō)明

單因素方差分析雙因素?zé)o交互作用方差分析

方差分析主要過(guò)程說(shuō)明ANOVA過(guò)程格式及使用說(shuō)明GLM過(guò)程格式及使用說(shuō)明PROCANOVA語(yǔ)句的選項(xiàng)主要有：

DATA=數(shù)據(jù)集名指明要分析的SAS數(shù)據(jù)集，缺省時(shí)SAS將使用最近建立的數(shù)據(jù)集.OUTSTAT=輸出數(shù)據(jù)集指定分析計(jì)算結(jié)果輸出的數(shù)據(jù)集名.CLASS語(yǔ)句指明分類變量，是ANOVA過(guò)程的必需語(yǔ)句，并且必須出現(xiàn)在MODEL語(yǔ)句之前.分類變量可以為數(shù)值型或字符型，分類變量的個(gè)數(shù)表示方差分析的因素個(gè)數(shù).MODEL語(yǔ)句定義分析所用的效應(yīng)模型，即方差分析的因變量和效應(yīng)變量.在方差分析過(guò)程中，關(guān)鍵在于定義線性數(shù)學(xué)模型，常用的模型定義語(yǔ)句有：MODELy=a單因素一元方差分析MODELy=ab雙因素?zé)o交互作用一元方差分析MODELy=aba*b雙因素有交互作用一元方差分析MEANS語(yǔ)句用來(lái)計(jì)算該語(yǔ)句所列的每個(gè)效應(yīng)所對(duì)應(yīng)的因變量均值，其選項(xiàng)用于設(shè)定多重比較的方法.過(guò)程說(shuō)明：GLM過(guò)程格式及使用說(shuō)明

GLM即廣義線性模型（GeneralLinerModel）過(guò)程，對(duì)于非平衡數(shù)據(jù)，應(yīng)采用GLM過(guò)程.它使用最小二乘法對(duì)數(shù)據(jù)擬合廣義線性模型.該過(guò)程功能強(qiáng)大，可用于多種不同的統(tǒng)計(jì)分析中.GLM過(guò)程用于方差分析時(shí)，主要語(yǔ)句和使用格式與上述ANOVA過(guò)程類似.單因素方差分析單因素等重復(fù)方差分析(ANOVA過(guò)程)單因素不等重復(fù)的方差分析（GLM過(guò)程）

應(yīng)用實(shí)例

一個(gè)工廠用三種不同的工藝生產(chǎn)某種電池.從三種工藝生產(chǎn)的電池中分別抽取5個(gè)樣品，測(cè)得樣品壽命的數(shù)據(jù)如下（單位小時(shí)）：A重復(fù)工藝1工藝2工藝314026392463440338304344228485443250平均423044我們要研究的指標(biāo)是電池的壽命，工藝是影響壽命的一個(gè)因素，三種工藝分別是該因素的三個(gè)水平.在試驗(yàn)中我們假設(shè)其它因素都處于相同的狀態(tài).這里我們希望利用上面得到的數(shù)據(jù)來(lái)考察“工藝”的不同是否對(duì)“壽命”這個(gè)指標(biāo)有影響？源程序

Dataexam;/*建立數(shù)據(jù)集*/Dotrt=1to3;/*3個(gè)處理（trt）分別為1、2、3*/DoI=1to5;/*每個(gè)處理下5次重復(fù)*/Inputx@@;Output;End;End;Cards;404638424426343028323940434850;Procanova;/*調(diào)用方差分析過(guò)程*/Classtrt;/*定義處理為分類變量*/Modelx=trt;/*定義效應(yīng)模型*/Title'方差分析';Run;運(yùn)行結(jié)果

AnalysisofVarianceProcedureDependentVariable:XSourceDFSumofSquaresMeanSquareFValuePr>FModel2573.33333333286.6666666719.770.0002Error12174.0000000014.50000000CT14747.33333333R-SquareC.V.RootMSEXMean0.7671729.8479823.8078865538.66666667

以上結(jié)果相當(dāng)于方差分析表，F(xiàn)值為19.77，顯著性水平為0.0002，小于0.01，說(shuō)明各處理間的均值差異極顯著.某食品公司對(duì)一種食品設(shè)計(jì)了四種新包裝.為了考察哪種包裝受歡迎，選取了十個(gè)有近似相同銷售量的商店作試驗(yàn)，其中兩種包裝各指定兩個(gè)商店銷售，另兩種包裝各指定三個(gè)商店銷售.在試驗(yàn)期間，各商店的貨架排放位置，空間都盡量一致，營(yíng)業(yè)員也采用相同的促銷方法.一段時(shí)間的銷售量記錄如表所示應(yīng)用實(shí)例

包裝類型商店1121419242181217303132130395754源程序Dm"LOG;CLEAR;OUTPUT;CLEAR;";Datanew;Inputstr$x@@;/*x表示含銷售量，str表示包裝處理*/Cards;A112A118A214A212A213A319A317A321A424A430;Procprint;Title"單因素不等重復(fù)方差分析";Procglm;Classstr;Modelx=str;Run;運(yùn)行結(jié)果

GeneralLinearModelsProcedureDependentVariable:XSourceDFSumofSquaresMeanSquareFValuePr>FModel3258.0000000086.0000000011.220.0071Error646.000000007.66666667CT9304.00000000R-SquareC.V.RootMSEXMean0.84868415.382642.7688746218.00000000

由方差分析表中看到，組間平方和為258，組內(nèi)平方和為46，總的平方和為304,F值為11.22，顯著性水平達(dá)到0.007，包裝不同對(duì)銷售量的影響極顯著.

在一個(gè)小麥種植試驗(yàn)中，考察4種不同肥料（因素A）與3種不同品種（因素B），選擇12塊形狀、大小等條件盡量一致的地塊，每塊以做一個(gè)處理，得如下數(shù)據(jù)，見(jiàn)表（單位：公斤/畝）.在水平下檢驗(yàn)假設(shè)：（1）使用不同肥料小麥平均產(chǎn)量無(wú)差異；（2）不同品種的小麥平均產(chǎn)量無(wú)差異.應(yīng)用實(shí)例

產(chǎn)量B小麥品種123肥A料種類1164172174215515714731591661584158157153源程序

Datanew;Doa='a1','a2','a3','a4';Dob='b1','b2','b3';Inputy@@;Output;End;End;Cards;164172174155157147159166158158157153;Procprintdata=new;Run;Procanova;Classab;Modely=ab;Run;運(yùn)行結(jié)果

AnalysisofVarianceProcedureDependentVariable:YSourceDFSumofSquaresMeanSquareFValuePr>FModel5554.00000000110.800000006.160.0234Error6108.0000000018.00000000CT11662.00000000R-SquareC.V.RootMSEYMean0.8368582.6516504.24264069160.00000000SourceDFAnovaSSMeanSquareFValuePr>FA3498.00000000166.000000009.220.0115B256.