本文關(guān)鍵詞：單因素方差分析的數(shù)學(xué)模型及其應(yīng)用，由筆耕文化傳播整理發(fā)布。

網(wǎng)友lu6966近日為您收集整理了關(guān)于單因素方差分析的數(shù)學(xué)模型及其應(yīng)用的文檔，希望對您的工作和學(xué)習(xí)有所幫助。以下是文檔介紹：2013年 6月第 27卷第 2期陰山學(xué)刊 YINSHAN ACADEMIC JOURNAL Jun.2013 V01.27 No.2 單因素方差分析的數(shù)學(xué)模型及其應(yīng)用田 兵(包頭師范學(xué)院《陰山學(xué)刊》編輯部,內(nèi)蒙古包頭 014030) 摘要:本文主要介紹了單因素方差分析的數(shù)學(xué)模型。同時結(jié)合實(shí)例演示了利用單因素方差分析通過 R 軟件編程解決實(shí)際問題的過程。關(guān)鍵詞:單因素方差分析;模型;平方和分解;檢驗統(tǒng)計量中圖分類號:0212.1 文獻(xiàn)標(biāo)識碼:A 文章編號:1004—1869(2013)02-0024一o4 在日常的工作和生活中,影響一件事的因素有很多,人們希望根據(jù)各種實(shí)驗來判斷不同的因素對實(shí)驗結(jié)果的影響。例如:不同的生產(chǎn)廠家、不同的原材料、不同的操作規(guī)程及不同的技術(shù)指標(biāo)等對產(chǎn)品的質(zhì)量、性能都會有影響,然而不同因素的影響大小不等。方差分析是研究一種或多種因素的變化對實(shí)驗結(jié)果的觀測值是否有顯著影響,從而找出較優(yōu)的試驗條件或生產(chǎn)條件的一種常用數(shù)理統(tǒng)計的方法。在實(shí)驗中所關(guān)注的數(shù)量指標(biāo)如產(chǎn)量、性能等稱為觀測值。影響觀測值的條件稱為因素。因素的不同狀態(tài)稱為水平,一個因素可以采用多個水平。在一項試驗中,可以得到一系列不同的觀測值。引起觀測值不同的原因是多方面的。有的是因為處理方式不同或條件不同引起的,稱作因素效應(yīng);有的是實(shí)驗過程中偶然性因素的干擾或觀測誤差所導(dǎo)致的,稱作試驗誤差。方差分析的主要工作是將測量數(shù)據(jù)的總變異按照變異原因的不同分解為因素效應(yīng)和試驗誤差,并對其作出數(shù)量分析,比較各種原因在總變異中所占的重要程度,作為統(tǒng)計推斷的依據(jù),由此確定進(jìn)一步的工作方向在研究實(shí)際問題時,我們通常是從最簡單的情形人手。單因素方差分析是方差分析的最簡單情形——影響其因素只有一個。 1 單因素方差分析的數(shù)學(xué)模型通常假設(shè)試驗只有一個因素 A在發(fā)生變化,其余的因素沒有變化。A有 r個水平 A A …,A,,在水平 A 下進(jìn)行凡次獨(dú)立觀測,得到試驗指標(biāo)如表 1所示。表 1:單因素方差分析數(shù)據(jù)其中 表示在因素 A的第 i個水平下的第次試驗的試驗結(jié)果。收稿日期:2013—05—21 作者簡介:田兵(1982一),男,山西五臺人,碩士,研究方向:數(shù)理統(tǒng)計。將水平 A 下的試驗結(jié)果 …, 看作是來自第 i正態(tài)總體置一N(/x , )的樣本觀測值,其中 , 都是未知的。而且對于每個總體置是相互獨(dú)立的,考慮線性統(tǒng)計模型』 +占, = ,2,…,r = ,2,…凡, (1) 【 —N(0, 。)且相互獨(dú)立, 其中 是第個總體的均值, f是相應(yīng)的試驗誤差。比較因素 A的 r個水平的差異歸結(jié)為比較這 r總體的均值。即檢驗假設(shè)風(fēng):/xl= 2=…=/x,,H1: l, 2,…, 不全相等。 (2) 記, r r 1 1f1 rt/]..ti,n ni, /.ti一這里表示總和的均值, 為水平Ai所對應(yīng)指標(biāo)的效應(yīng)。因此有∑trio/ =0. 上述模型又可以等價的寫成 r if= + i+8 ,,=1,2,…,r,.,=1,2,…,n , J ~N(O, )且相互獨(dú)立, (3) l, 【∑ i=0. 稱模型(3)為單因素方差分析的數(shù)學(xué)模型,其是一種線性模型。 2 方差分析假設(shè)(2)式等價于風(fēng): 1= 2=…= ,=0,H1: l, 2,…, ,不全為零。如果風(fēng)被拒絕,那么就說明因素 4的個水平的效應(yīng)之間有顯著的差異;否則,差異不明顯。為了導(dǎo)出 的檢驗統(tǒng)計量,方差分析法建立在平方和分解和自由度分解的基礎(chǔ)上,考察統(tǒng)計量 r ni 1 r ni |s,=∑∑( 一) , =÷∑∑‘ 1 J I ’ ‘ I , l 稱|s,為總離差平方和(或稱為總變差),其實(shí)全部試驗數(shù)據(jù) 與總平均值差的平方和,描述了所有觀測數(shù)據(jù)的離散程度,可以證明如下的平方和分解公式: Sr=SE+SA (4) 其中 r ni . ni r i r .s =∑∑( 一) ,扎=÷∑Xij,S =∑∑( 一) =∑rti( 一) . 這里 Js 是代表隨機(jī)誤差的影響。這是因為對于固定的 i來講,觀測值 ,…, 是來自同一個正態(tài)總體 ni N(I.t ,0.2)的樣本,因此,它們之間的差異是有隨機(jī)誤差多引起的。而∑( 一.) 是這n 個數(shù)據(jù)的變動平方和,正是它們差異大小的度量。將 r組這樣的變動平方和相加,就得到了 Js ,一般稱 Is 為誤差平方或組內(nèi)平方。.s 表示在A 水平下的樣本均值和總平均值之間的差異之和,它反映了r個總體均值之間的差異,因為. 是第個總體的樣本均值,是 的估計,因此 r個總體均值。, ,…, 之間的差異越大,這些樣本均值, :, …, 之間的差異也就越大。平方和∑17,i( i.一) 正是這種差異大小的度量,這里ni反映了第i個總體樣本大小在平方和.s 中的作用,稱 Is 為因素 A的效應(yīng)平方和或組間平方和。式(4)表明,總平方和 Js 可按其來源分解成兩部分,一部分是誤差平方和 5 ,是由隨機(jī)誤差所致;另一部分是因素 A的平方和 5 ,是由因素 的各水平的差異一起的。 25 由模型假設(shè)(2)經(jīng)過統(tǒng)計分析可以得到 E(S )= (n—r)or ,即 Is /( —r)是 or 的一個無偏估計,且~ z(n一). 如果原假設(shè) 成立,則有,即此時 S /(r一1)也是的無偏估計,且一 (r一1), 并且 s 與 S 相互獨(dú)立,因此當(dāng) 成立時有 F = ~F(r-1,n-r . 于是 F可以作為風(fēng)的檢驗統(tǒng)計量,對給定的顯著性水平,用 F (r一1,n—r)表示 F分布的上 ot分為點(diǎn)。若 F >F (r一1,n—r),則拒絕原假設(shè),認(rèn)為因素 A的 r個水平有顯著差異。也可以通過計算 P值的方法來決定是接受還是拒絕原假設(shè) Ho。P值為P=P{F(r一1,n—r)>F},它表示的是服從自由度為(r一1,n— r)的 F分布的隨機(jī)變量取值大于 F的概率。顯然,P值小于 ot等價于 F>F (r一1,n—r),表示在顯著性水平 o/下的小概率事件發(fā)生了。這意味著應(yīng)該拒絕原假設(shè) 。當(dāng)P值大于,則無法拒絕原假設(shè) ,所以應(yīng)接受原假設(shè) 。將上述分析整理成表的形式,就可以得到方差分析表 2。表 2:單因素方差分析表 3 相關(guān)實(shí)例 3個工廠生產(chǎn)同一種零部件�，F(xiàn)在從各廠產(chǎn)品中分別出取 4件產(chǎn)品做檢測,其檢測強(qiáng)度見表 3所示。表 3:產(chǎn)品檢測數(shù)據(jù)判斷 3個廠生產(chǎn)的產(chǎn)品零件強(qiáng)度是否有顯著差異。利用 軟件解決上述問題,相應(yīng)的 軟件計算過程如下。解:設(shè)零部件的強(qiáng)度為為所考察的因素。3個不同工廠生產(chǎn)能力為 3個水平,從各廠取出的產(chǎn)品所做得檢測強(qiáng)度值,視為來自3個正態(tài)分布總體的樣本觀測值。問題歸結(jié)為檢驗: 風(fēng):IZl= 2= 3= ,Hl:/zl, 2, 3, 4不全相等。程序如下: > intensity< 一data.framef X=C(115,l16,98,83,120,1l4,99,101,106,107,103, 108,107,118,116,104,117,121,111,10

12>

  本文關(guān)鍵詞：單因素方差分析的數(shù)學(xué)模型及其應(yīng)用，由筆耕文化傳播整理發(fā)布。