【導(dǎo)讀】方差分析又稱變異數(shù)分析，可簡記為ANOVA，主要用于檢驗計量。資料中的兩個或兩個以上均值間差別顯著性的方法。當(dāng)欲比較幾組均值時，理論上抽得的幾。個樣本，都假定來自正態(tài)總體，且有一個相同的方差，僅僅均值可以不相同。的含義，稱之謂效應(yīng)的可加性。所謂的方差是離均差平方和除以自由度，在方差分析中常簡。各部分的均方，然后列出方差分析表算出F值，作出統(tǒng)計推斷。種效應(yīng)的作用就越了解，統(tǒng)計推斷就越準(zhǔn)確。方差分析表的一般形式如表所示。項的F值，F(xiàn)值實際上是兩個均方之比值，通常情況下，分母的均方是誤差項的均方?？赡軐傋儺愑袑嵸|(zhì)影響的。要求進行試驗，得出原始觀察值，按原來設(shè)計方案算出方差分析表中的各項。一次試驗時可以有多個觀察指標(biāo)，方差分析時也可以同時對多個因變量進。是地區(qū)這一因素的四個水平。水平?jīng)Q定了受試對象各自應(yīng)該歸屬的組別。①每一總體均服從正態(tài)分布；為了能更仔細地描述數(shù)據(jù)，常在方差分析模型中引入一般平均與效應(yīng)的概念。

　　

【正文】 Means with the same letter are not significantly different. Duncan Grouping Mean N BRAND A 4 TUFFY A B A 4 XTRA B A B A 4 CHAMP B B 4 ACME C 4 AJAX T Confidence Intervals for variable: WEAR Alpha= Confidence= df= 15 MSE= Critical Value of T= Half Width of Confidence Interval= Lower Upper BRAND N Confidence Mean Confidence Limit Limit TUFFY 4 XTRA 4 CHAMP 4 ACME 4 AJAX 4 上海財經(jīng)大學(xué)經(jīng)濟信息管理系 IS/SHUFE Page 20 of 39 表 （ c）中結(jié)果分析： 注意在顯著水平為 上，兩兩比較的最小顯著差為 ，如果顯著則被標(biāo)上“ ***”。例如， TUFFY 均值減 XTRA 均值 =－ =，顯著。綜合分析的結(jié)果表明， AJAX 品牌均值顯著與其他品牌均值不同，且為最小的均值；TUFFY 品牌均值也顯著與其他品牌均值不同，且為最大的均值； XTRA、 CHAMP、 ACME三個品牌均值之間無顯著差異。 3． 有計劃的均值比較和參數(shù)估計 例 繼續(xù)上例的分析。有時在實際情況中，多重比較要按某種分類標(biāo)準(zhǔn)來進行，例如，假設(shè)我們知道 5 種品牌的制造商情況，品牌 ACMX、 AXAX 和 CHAMP 來自美國 .制造商，而品牌 TUFFY 和 XTRA 來自非美國 。 假設(shè) 我們有興趣比較美國品牌的均值 與非美國品牌的均值是否有差異。 程序如下： proc glm data=。 class brand。 model wear=brand。 The SAS System T tests (LSD) for variable: WEAR NOTE: This test controls the type I parisonwise error rate not the experimentwise error rate. Alpha= Confidence= df= 15 MSE= Critical Value of T= Least Significant Difference= Comparisons significant at the level are indicated by 39。***39。. Lower Difference Upper BRAND Confidence Between Confidence Comparison Limit Means Limit TUFFY XTRA *** TUFFY CHAMP *** TUFFY ACME *** TUFFY AJAX *** XTRA TUFFY *** XTRA CHAMP XTRA ACME XTRA AJAX *** CHAMP TUFFY *** CHAMP XTRA CHAMP ACME CHAMP AJAX *** ACME TUFFY *** ACME XTRA ACME CHAMP ACME AJAX *** AJAX TUFFY *** AJAX XTRA *** AJAX CHAMP *** AJAX ACME *** 上海財經(jīng)大學(xué)經(jīng)濟信息管理系 IS/SHUFE Page 21 of 39 contrast 39。US vs .39。 brand 2 2 2 3 3。 estimate 39。US vs .39。 brand 2 2 2 3 3。 run。 程序說明：使用 contrast 語句來產(chǎn)生有計劃的均值比較分析和使用 estimate 語句進行參數(shù)估計。注意在 anova過程中沒有這兩條語句，必須使用 glm 過程。使用 contrast 語句前，應(yīng)該首先表達出所關(guān)心的均值線性組合的原假設(shè)，如 ： 0)(3)(2:)(2/1)(3/1:00?????????X T R AT U F F YC H A M PA J A XA C M EX T R AT U F F YC H A M PA J A XA C M EHH??????????等價于 contrast語句的三個 基本參數(shù)，一是標(biāo)簽（ 39。US vs .39。） ，二是效應(yīng)名（ brand），三是效應(yīng)的數(shù)字系數(shù)表（ 2 2 2 3 3）。應(yīng)特別注意的是，數(shù)字系數(shù)的次序是匹配分類變量按字母數(shù)字次序的水平值。事實上，均值線性組合的系數(shù)同樣是 model 語句中效應(yīng)參數(shù)組合的系數(shù)，這是因為， ii ??? ?? ，將它們分別代入均值線性組合 后，可得到 ： )(3)(2)(3)(2)(3)(2X T R AT U F F YX T R AA J A XA C M EX T R AT U F F YX T R AC H A M PA J A XA C M EX T R AT U F F YC H A M PA J A XA C M E??????????????????????????????????????? 所以， estimate 語句的使用格式與 contrast 語句非常類同。 輸出的主要結(jié)果見表 。 表 有計劃的均值比較和參數(shù)估計 表 ：顯示了美國品牌均值與非美國品牌均值比較的平方和為 ， F值為 13=，這個 F（ 1， 15）分布 F 值大于 13的概率為 ，因此原假設(shè)是顯著的，拒絕接受，即美國品牌均值與非美國品牌 均值是不同的。效應(yīng)組合的參數(shù)估計 為 =3（ ++）－ 2（ +），對于原假設(shè)參數(shù)是否為 0的 t 檢驗， t統(tǒng)計量為 ，概率為 ，拒絕接受。注意到 t 檢驗的 p 值為 ，與對比分析的 F 檢驗的 p 值相同，這是因為兩種檢驗是相同的， F 值等于 t 的平方。 4． 隨機單位組試驗設(shè)計的方差分析 例 某食品公司對一種食品設(shè)計了四種包裝。為了考察哪種包裝最受歡迎，選了十個有近似相同銷售量的商店作試驗，其中兩種包裝各指定兩個商店，另兩種包裝各 指定三個商店銷售。在試驗期中各商店的貨架排放位置、空間都盡量一致，營業(yè)員的促銷方法也基本相同。觀察在一定時期的銷售量，數(shù)據(jù)見表 。試比較四種包裝的銷售量是否一致。 表 四種包裝在 10 個商店中的銷售量 包裝類型 商店（ block） 商店數(shù) The SAS System Contrast DF Contrast SS Mean Square F Value Pr F US vs . 1 T for H0: Pr |T| Std Error of Parameter Estimate Parameter=0 Estimate US vs . 上海財經(jīng)大學(xué)經(jīng)濟信息管理系 IS/SHUFE Page 22 of 39 （ treat） 1 2 3 n A1 12 18 2 A2 14 12 13 3 A3 19 17 21 3 A4 24 30 2 程序如下： data pack。 input treat $ n 。 do block=1 to n。 input y @@。 output。 end。 cards。 A1 2 12 18 A2 3 14 12 13 A3 3 19 17 21 A4 2 24 30 。 proc glm data=pack。 class block treat。 model y=block treat。 means block treat/snk。 means block treat/dunt(39。139。)。 means block treat。 run。 程序說明：由于包裝 類型 A1 和 A4 在商店 3里沒有進行試驗，所以產(chǎn)生的是一種不平衡的試驗數(shù)據(jù)，請注意隨機單位區(qū)組中有不平衡數(shù)據(jù)集時，程序應(yīng)該如何設(shè)計