第十一方差分析演示文稿

第十一方差分析演示文稿當(dāng)前1頁(yè)，總共47頁(yè)。（優(yōu)選）第十一方差分析當(dāng)前2頁(yè)，總共47頁(yè)。一、單因素方差分析

（完全隨機(jī)設(shè)計(jì)資料的方差分析）

One-wayANOVA當(dāng)前3頁(yè)，總共47頁(yè)。

例1為了解燙傷后不同時(shí)期切痂對(duì)肝臟三磷酸腺苷（ATP）含量的影響，將30只雄性大隨機(jī)分為3組，每組10只；A組為燙傷組，B組為24小時(shí)切痂組，C組為96小時(shí)切痂組。全部動(dòng)物統(tǒng)一在燙傷后168小時(shí)處死，并測(cè)定其肝臟的ATP含量，結(jié)果如下表，問(wèn)不同時(shí)期切痂對(duì)ATP含量有無(wú)影響？當(dāng)前4頁(yè)，總共47頁(yè)。A組B組C駔Xij7.766.6711.1414.1610.858.817.7111.7311.606.948.588.228.435.7811.4213.017.199.958.476.6113.8514.189.3611.2610.306.9713.5317.729.598.68例數(shù)ni10101030(N)均數(shù)8.0412.769.2518合計(jì)80.43127.5592.49300.47(∑X)平方和676.321696.96868.933242.21(∑X2)合計(jì)當(dāng)前5頁(yè)，總共47頁(yè)。幾個(gè)基本術(shù)語(yǔ)1處理因素：將所研究的對(duì)象分為多個(gè)處理組,施加不同的干預(yù)，施加的干預(yù)稱為處理因素（factor），例如影響農(nóng)作物產(chǎn)量的因素有氣溫、降雨量、日照時(shí)間等。2水平：因素的不同取值等級(jí)稱作水平，例如性別有男、女兩個(gè)水平。處理因素至少有兩個(gè)水平(level)。用這類資料的樣本信息來(lái)推斷各處理組間多個(gè)總體均數(shù)是否存在差別，常采用的統(tǒng)計(jì)分析方法為方差分析(analysisofvariance,ANOVA)。由英國(guó)統(tǒng)計(jì)學(xué)家首創(chuàng)，為紀(jì)念Fisher，以F命名，故方差分析又稱F檢驗(yàn)（Ftest）。當(dāng)前6頁(yè)，總共47頁(yè)。離均差離均差和離均差平方和6方差幾個(gè)基本術(shù)語(yǔ)當(dāng)前7頁(yè)，總共47頁(yè)?？傋儺愔纸饨M內(nèi)誤差組間誤差當(dāng)前8頁(yè)，總共47頁(yè)?？傋儺悾═otalvariation）：全部測(cè)量值Xij與總均數(shù)間的差別

組間變異（betweengroupvariation）各組的均數(shù)與總均數(shù)間的差異組內(nèi)變異（withingroupvariation)每組的10個(gè)原始數(shù)據(jù)與該組均數(shù)的差異

試驗(yàn)數(shù)據(jù)有三個(gè)不同的變異下面先用離均差平方和(sumofsquaresofdeviationsfrommean，SS)表示變異的大小

當(dāng)前9頁(yè)，總共47頁(yè)。1.總變異SS總反映了所有測(cè)量值之間總的變異程度，

SS總=各測(cè)量值Xij與總均數(shù)差值的平方和當(dāng)前10頁(yè)，總共47頁(yè)。SS組間反映了各組均數(shù)間的變異程度組間變異＝①隨機(jī)誤差+②處理因素效應(yīng)

2.組間變異mi

mj當(dāng)前11頁(yè)，總共47頁(yè)。

在同一處理組內(nèi)，雖然每個(gè)受試對(duì)象接受的處理相同，但測(cè)量值仍各不相同，這種變異稱為組內(nèi)變異。SS組內(nèi)僅僅反映了隨機(jī)誤差的影響。也稱SS誤差3.組內(nèi)變異

m

i當(dāng)前12頁(yè)，總共47頁(yè)。2、變異來(lái)源變異原因總

變異

組內(nèi)變異組間變異處理因素隨機(jī)誤差TreatmentErorr?三種變異之間的關(guān)系當(dāng)前13頁(yè)，總共47頁(yè)。均方(meansquare，MS)當(dāng)前14頁(yè)，總共47頁(yè)。均方之比＝Fvalue當(dāng)前15頁(yè)，總共47頁(yè)。F分布F分布概率密度函數(shù)：當(dāng)前16頁(yè)，總共47頁(yè)。F分布曲線當(dāng)前17頁(yè)，總共47頁(yè)。F界值表附表4F界值表（方差分析用，單側(cè)界值）上行：P=0.05下行：P=0.01分母自由度υ2分子的自由度，υ1123456

1161200216225230234

405249995403562557645859

218.5119.0019.1619.2519.3019.33

98.4999.0099.1799.2599.3099.33

254.243.392.992.762.602.49

7.775.574.684.183.853.63

當(dāng)前18頁(yè)，總共47頁(yè)。

(二)、分析步驟

1、建立假設(shè)：2、計(jì)算檢驗(yàn)統(tǒng)計(jì)量F3、求P值、下結(jié)論當(dāng)前19頁(yè)，總共47頁(yè)。1、建立假設(shè)、確定檢驗(yàn)水準(zhǔn)H0：三個(gè)均總體數(shù)全部相等

1=2=3H1：三個(gè)總體均數(shù)全不相等或不全相等

=0.05當(dāng)前20頁(yè)，總共47頁(yè)。2、計(jì)算檢驗(yàn)統(tǒng)計(jì)量計(jì)算各組的∑X、∑X2及總的∑X和∑X2。計(jì)算校正數(shù)C

式中N為各組樣本含量之和。本例C=（300.47）2/30=3009.4074當(dāng)前21頁(yè)，總共47頁(yè)。、計(jì)算總的變異及總的自由度=30-1=29當(dāng)前22頁(yè)，總共47頁(yè)。、計(jì)算組間變異及相應(yīng)的自由度=3-1=2當(dāng)前23頁(yè)，總共47頁(yè)。、計(jì)算組內(nèi)變異及相應(yīng)的自由度

=232.8026-119.8314=112.9712=30-3=27當(dāng)前24頁(yè)，總共47頁(yè)。

、列出方差分析表變異SSMSFP總變異232.80

29組間變異119.83259.9214.32P<0.05誤差變異112.97274.18當(dāng)前25頁(yè)，總共47頁(yè)。3、確定P值、下結(jié)論

從上表得F=14.32，查附表4（方差分析界值表，單側(cè)），自由度相同時(shí)，P值越大,F界值越小。因F0.05，2，27=3.35；故P<0.05，按=0.05水準(zhǔn)拒絕H0，可認(rèn)為三個(gè)不同時(shí)期切痂對(duì)ATP含量有影響。當(dāng)前26頁(yè)，總共47頁(yè)。方差分析與t檢驗(yàn)的比較兩個(gè)樣本均數(shù)的比較可用方差分析，計(jì)算得到的統(tǒng)計(jì)量F值是用t檢驗(yàn)計(jì)算的t統(tǒng)計(jì)量的平方，而F,(1,2)=t2,，因此兩種檢驗(yàn)方法的結(jié)果是完全一致的。多個(gè)樣本均數(shù)的比較不能用t檢驗(yàn)。當(dāng)前27頁(yè)，總共47頁(yè)。完全隨機(jī)設(shè)計(jì)的變異分解總變異組間變異

組內(nèi)變異（誤差變異）當(dāng)前28頁(yè)，總共47頁(yè)。二、配伍組方差分析（隨機(jī)區(qū)組方差分析）當(dāng)前29頁(yè)，總共47頁(yè)。例11-3:為比較不同濃度的血水草總生物堿對(duì)血吸蟲尾蚴的殺滅作用，試驗(yàn)動(dòng)物為48只雌性小鼠。試驗(yàn)方法是將小鼠感染40只血吸蟲尾蚴后，分別接受4種處理，其中甲處理為對(duì)照，其余3種處理分別為不同濃度的血水草總生物堿濃度，用試驗(yàn)后小鼠體內(nèi)的尾蚴的存活率評(píng)價(jià)不同濃度的血水草總生物堿的作用。可能影響試驗(yàn)結(jié)果的因素是小鼠體重，試分析。當(dāng)前30頁(yè)，總共47頁(yè)。不同濃度血水草總生物堿處理后小鼠體內(nèi)尾蚴存活率當(dāng)前31頁(yè)，總共47頁(yè)。配伍組設(shè)計(jì)的變異分解總變異組間變異配伍組變異

處理?+隨機(jī)配伍?+隨機(jī)誤差變異(隨機(jī))當(dāng)前32頁(yè)，總共47頁(yè)。(一)、建立假設(shè)、確定檢驗(yàn)水準(zhǔn)處理間：H0：1=2=3=4（4種處理組小鼠體內(nèi)的尾蚴存活率的總體均數(shù)相等）H1：1、2、3不相等或不全相等

=0.05當(dāng)前33頁(yè)，總共47頁(yè)。（二）、計(jì)算檢驗(yàn)統(tǒng)計(jì)量F1、先計(jì)算校正數(shù)C：2、計(jì)算SS總及總的自由度當(dāng)前34頁(yè)，總共47頁(yè)。3、計(jì)算SS處理和SS配伍當(dāng)前35頁(yè)，總共47頁(yè)。4、列出方差分析表變異SSMSFP總變異2.054947處理組間0.947830.315916.7294P<0.01區(qū)組間0.4840110.0442.33P<0.05誤差變異0.6231330.0189當(dāng)前36頁(yè)，總共47頁(yè)。（三）、確定P值、下結(jié)論處理間差別的推斷：因P<0.01，按=0.05水準(zhǔn)拒絕H0，故可認(rèn)4個(gè)處理組的小體內(nèi)尾蚴存活率的總體均數(shù)不全相等。區(qū)組間P<0.01，按=0.05水準(zhǔn)拒絕H0，故認(rèn)為各區(qū)組間的總體均數(shù)有差別（體重對(duì)小鼠體內(nèi)尾蚴存活率有影響）。此設(shè)計(jì)減少了誤差，較之完全隨機(jī)設(shè)計(jì)，試驗(yàn)效率提高了。當(dāng)前37頁(yè)，總共47頁(yè)。配伍組設(shè)計(jì)的基本思想總變異組間變異誤差變異配伍組變異組內(nèi)變異當(dāng)前38頁(yè)，總共47頁(yè)。多個(gè)均數(shù)兩兩比較如果有K組,則有k(k-1)/2對(duì)均數(shù)需作比較,這些比較并非獨(dú)立的,如用t檢驗(yàn)作兩兩比較,勢(shì)必增大Ⅰ類誤差。宜用多重比較(multiplecomparison）方法。具體方法較多，最常見的是q檢驗(yàn)法/SNK法。當(dāng)前39頁(yè)，總共47頁(yè)。Student-Neman-Keuls（SNK）檢驗(yàn)各處理組的樣本含量ni相等時(shí)：ni不相等：當(dāng)前40頁(yè)，總共47頁(yè)。H0：任兩對(duì)比組的總體均數(shù)相等即A=BH1：A≠B

=0.05當(dāng)前41頁(yè)，總共47頁(yè)。將3個(gè)樣本均數(shù)從小到大順序排列，并編上秩次：均數(shù)

8.049.2512.76

組別A組C組B組秩次123

列出兩兩比較計(jì)算表，共進(jìn)行3次兩兩比較。當(dāng)前42頁(yè)，總共47頁(yè)。多個(gè)均數(shù)兩兩比較q值表當(dāng)前43頁(yè)，總共47頁(yè)。二LSD-t檢驗(yàn)

（Leastsignificantdifference）當(dāng)前44頁(yè)，總共47頁(yè)。24小時(shí)切痂組和對(duì)照組的比較建立假設(shè)并確定檢驗(yàn)水準(zhǔn)。H0：24h=0H1：24h≠0H0：=0.05;求LSD-t值MSE=4.184,v=27,=0.05水準(zhǔn)，拒絕H0當(dāng)前