《spss方差分析》課件

流行病學(xué)與衛(wèi)生統(tǒng)計(jì)學(xué)教研室1方差分析2方差分析單因素方差分析廣義線性模型AnalyzeCompareMeansOne-WayANOVAAnalyzeGeneralLinearModelUnivariate3單因素方差分析理論復(fù)習(xí)多個(gè)樣本均數(shù)比較時(shí)不能用兩兩t檢驗(yàn)方差分析的基本思想是將變異分解總變異=隨機(jī)變異+處理因素導(dǎo)致的變異要求：正態(tài)性、方差齊性組內(nèi)變異反映隨機(jī)變異組間變異則是處理因素影響造成4單因素方差分析應(yīng)用案例某職業(yè)病防治院對31名石棉礦工中的石棉肺患者、可疑患者及非患者進(jìn)行了用力肺活量（L）測定，問三組石棉礦工的用力肺活量有無差別？（lung.sav）H0及H1對適用條件的檢查：正態(tài)性、方差齊性方差分析表的閱讀單因素方差分析Analyze——CompareMeans——One-WayANOVA單因素方差分析單因素方差分析8單因素方差分析練習(xí)某醫(yī)生為了研究一種四類降糖新藥的療效，以統(tǒng)一的納入和排除標(biāo)準(zhǔn)選擇了60名2型糖尿病患者，按完全隨機(jī)設(shè)計(jì)方案將患者分為三組進(jìn)行雙盲臨床試驗(yàn)。其中，降糖新藥高劑量組21人、低劑量組19人、對照組20人。對照組服用公認(rèn)的降糖藥物，治療4周后測得其餐后兩小時(shí)血糖的下降值，問治療4周后，餐后2小時(shí)血糖下降值的三組總體平均水平是否不同？（jiangtang.sav）9單因素方差分析當(dāng)多組間存在差異，還需要進(jìn)行兩兩比較兩兩比較方法的選擇：事先計(jì)劃好時(shí)：LSD、Bonferroni

事先沒有計(jì)劃好：試驗(yàn)組與對照組比較Dunnett

樣本含量相同時(shí)Tukey

樣本含量不同時(shí)Scheffe結(jié)果閱讀LSD法、SNK法、均數(shù)圖單因素方差分析PostHoc

對話框單因素方差分析單因素方差分析13單因素方差分析各組間精細(xì)比較（contrast按鈕）多重比較歸結(jié)為對均數(shù)的線性組合中L＝0的假設(shè)檢驗(yàn)。 如：a1,a2,a3分別為1，－1，0，就相當(dāng)于對第一、第二組均數(shù)是否相等的兩兩比較。

a1,a2,a3分別為1，-1/2，-1/2，可以檢驗(yàn)第一組均數(shù)是否為后兩組均數(shù)的平均。單因素方差分析單因素方差分析16單因素方差分析應(yīng)用案例比較4種強(qiáng)心劑的毒性，將稀釋過的藥物注入豚鼠的體內(nèi)，記錄導(dǎo)致豚鼠死亡時(shí)藥物的劑量，每種藥物各進(jìn)行了10次試驗(yàn)，事先已經(jīng)采取了有關(guān)措施以保證每只豚鼠試驗(yàn)時(shí)環(huán)境和測量方法盡量相同，問這四種強(qiáng)心劑的效力是否相同？第一種藥物的效力是否相當(dāng)于后三種藥的效力的平均？（qiangxinji.sav）17廣義線性模型18廣義線性模型模型概念(單因素)

總變異＝處理因素導(dǎo)致的變異＋隨機(jī)變異 其中Yij表示第i組的第j個(gè)觀察值；μ表示總體的平均水平；αi表示影響因素在i水平下對隨機(jī)變量的附加效應(yīng)，假設(shè)所有αi之和為0；εij為服從正態(tài)分布的隨機(jī)變量，代表隨機(jī)誤差。 一般檢驗(yàn)的是檢驗(yàn)各個(gè)αi是否均為0，即原假設(shè)為各組均數(shù)相等。19廣義線性模型分析實(shí)例強(qiáng)心劑類型與死亡時(shí)藥物劑量（qiangxinji.sav）廣義線性模型Analyze－GeneralLinealModel－Univariate廣義線性模型第一行為對整個(gè)模型的檢驗(yàn)，第二行為截距的檢驗(yàn)，第三行為對模型中各因素的檢驗(yàn)23廣義線性模型模型概念（兩因素）

αi、βj分別表示A因素i水平和B因素j水平的附加效應(yīng)，αiβj則為二者的交互效應(yīng)。若要說明因素A有無影響，即檢驗(yàn)如下假設(shè)：

H0：αi＝0，H1：至少有一個(gè)αi≠0若要說明因素B有無影響，即檢驗(yàn)如下假設(shè)：

H0：βj＝0，H1:至少有一個(gè)βj≠0模型無統(tǒng)計(jì)學(xué)意義則是指上面兩個(gè)H0同時(shí)成立24廣義線性模型分析實(shí)例（隨機(jī)區(qū)組設(shè)計(jì)）為探索丹參對肢體卻血再灌注損傷的影響，將30只純種新西蘭實(shí)驗(yàn)用大白兔，按窩別相同、體重相近劃分為10個(gè)區(qū)組。每個(gè)區(qū)組3只大白兔隨機(jī)采用A、B、C三種處理方案，即在松止血帶前分別給與丹參2ml/kg、丹參1ml/kg、生理鹽水2ml/kg，在松止血帶前及松后1小時(shí)分別測定血中白蛋白含量（g/L），問各區(qū)組間及各處理間無有差異。（danshen.sav）廣義線性模型Analyze－GeneralLinealModel－Univariate廣義線性模型廣義線性模型考慮了交互作用的結(jié)果如下：廣義線性模型不考慮交互作用：Model對話框選擇maineffects廣義線性模型不考慮交互作用的結(jié)果如下：廣義線性模型輪廓圖（Plot對話框）廣義線性模型廣義線性模型33廣義線性模型練習(xí)為了研究克拉霉素的抑菌效果，某實(shí)驗(yàn)室對28個(gè)短小芽孢桿菌平板依據(jù)菌株的來源不同分成了7個(gè)區(qū)組，每組4個(gè)平板用隨機(jī)的方式分配給標(biāo)準(zhǔn)藥物高劑量組（SH）、標(biāo)準(zhǔn)藥物低劑量組（SL），以及克拉霉素高劑量組（TH）、克拉霉素低劑量組（TL）。給予不同的處理后，觀察抑菌圈的直徑。試對該資料進(jìn)行分析。（kelameisu.sav）析因設(shè)計(jì)的方差分析析因?qū)嶒?yàn)設(shè)計(jì)是將每個(gè)因素的所有水平都互相組合，交叉分組進(jìn)行實(shí)驗(yàn)，可尋找最佳組合。析因設(shè)計(jì)的方差分析2×2析因例題利血平可以使小鼠腦中去甲腎上腺素（NE）等遞質(zhì)下降，現(xiàn)考察某種新藥MWC是否具有對抗利血平使遞質(zhì)下降的作用，將24只小鼠隨機(jī)等分為四組，并給予不同處理后，測定腦中NE的含量。析因設(shè)計(jì)的方差分析小鼠經(jīng)不同處理后腦中NE的含量蒸餾水組利血平組MWC組利血平＋MWC組630181715407760103663397687138638378676141887363892197625438523193648412析因設(shè)計(jì)的方差分析Analyze－GeneralLinealModel－Univariate析因設(shè)計(jì)的方差分析析因設(shè)計(jì)的方差分析重復(fù)測量資料的方差分析重復(fù)測量資料的定義同一受試對象的同一觀察指標(biāo)在不同時(shí)間點(diǎn)上進(jìn)行多次測量所得的資料。41廣義線性模型分析實(shí)例（重復(fù)測量資料）為研究減肥新藥鹽酸西布曲明片和鹽酸西布曲明膠囊的減肥效果是否不同，以及肥胖患者服藥后不同時(shí)間的體重隨時(shí)間的變化情況。將體重指數(shù)BMI≥27的肥胖患者40名隨機(jī)等分成兩組，一組給予鹽酸西布曲明片＋模擬鹽酸西布曲明膠囊，另一組給予鹽酸西布曲明膠囊＋模擬鹽酸西布曲明片。所有患者每天堅(jiān)持服藥，共服藥24周，期間禁用任何影響體重的藥物，且行為、飲食及運(yùn)動(dòng)與服藥前的平衡期保持一致。分別與平衡期（0周）、服藥后8周、16周、24周測定肥胖患者的體重。（jianfei.sav）結(jié)果解釋廣義線性模型Analyze－GeneralLinealModel－Repeatedmeasures預(yù)定義對話框：重復(fù)測量的指標(biāo)為受試者內(nèi)因素（Within-subjectFactor），重復(fù)測量的次數(shù)為受試者內(nèi)因素的水平數(shù)廣義線性模型廣義線性模型下表為受試者內(nèi)因素、受試者內(nèi)因素與自變量的一級交互作用的多元方差分析統(tǒng)計(jì)學(xué)檢驗(yàn)結(jié)果。廣義線性模型球