統(tǒng)計分析思路課件

統(tǒng)計分析思路

及案例分析常用統(tǒng)計方法選擇思路反應變量是單變量，雙變量還是多變量資料類型單因素還是多因素單樣本、兩樣本還是多樣本是否配對、配伍設計、重復測量是否滿足檢驗方法的前提條件單反應變量計量資料的分析例1

某醫(yī)生測量了36名從事鉛作業(yè)男性工人的血紅蛋白含量（原始數(shù)據(jù)略），算得其均數(shù)為130.83g/L，標準差為25.74g/L。問從事鉛作業(yè)工人的血紅蛋白是否不同于正常成年男性平均值140g/L？

思路：單反應變量：血紅蛋白含量資料類型：計量資料單因素：是否從事鉛作業(yè)工人一份樣本該設計為樣本均數(shù)與總體均數(shù)的比較分析方法：

如果服從正態(tài)分布，用單樣本t檢驗（onesamplettest)

如果不服從正態(tài)分布，非參數(shù)檢驗方法（WilcoxonsignedranktestWilcoxononesampletest）例2為比較兩種方法對乳酸飲料中脂肪含量測定結(jié)果是否不同，隨機抽取了10份乳酸飲料制品，分別用脂肪酸水解法和哥特里－羅紫法測定其結(jié)果如下表第(1)~(3)欄。問兩法測定結(jié)果是否不同？思路：單反應變量：脂肪含量資料類型：計量資料單因素：脂肪含量的測定方法兩份樣本配對設計分析方法：差值如果服從正態(tài)分布，用配對t檢驗（pairedttest)

如果不服從正態(tài)分布，非參數(shù)檢驗方法（WilcoxonsignedranktestWilcoxononesampletest）例3為研究國產(chǎn)四類新藥阿卡波糖膠囊的降血糖效果，某醫(yī)院用40名Ⅱ型糖尿病病人進行同期隨機對照試驗。試驗者將這些病人隨機等分到試驗組(用阿卡波糖膠囊)和對照組(用拜糖平膠囊)，分別測得試驗開始前和8周后的空腹血糖，算得空腹血糖下降值見表，能否認為該國產(chǎn)四類新藥阿卡波糖膠囊與拜糖平膠囊對空腹血糖的降糖效果不同？

思路：單反應變量：血糖下降值資料類型：計量資料單因素：不同藥物兩份樣本兩獨立樣本設計

分析方法：如果兩樣本都服從正態(tài)分布，且方差齊性，用兩樣本t檢驗（twosamplettest)；如果方差不齊，但服從正態(tài)分布，用t’檢驗，或兩樣本秩和檢驗（Wilcoxontwosampletest/MannWhitneytest)；皆不滿足，兩樣本秩和檢驗（Wilcoxontwosampletest/MannWhitneytest)例4某醫(yī)生為了研究一種降血脂新藥的臨床療效，按統(tǒng)一納入標準選擇120名高血脂患者，采用完全隨機設計方法將患者等分為4組（安慰劑組和不同劑量組），進行雙盲試驗。6周后測得低密度脂蛋白(mmol/L)，問4組患者的低密度脂蛋白含量有無差別?例58名受試對象在相同實驗條件下分別接受4種不同頻率聲音的刺激，他們的反應率（%）資料見表。問4種頻率聲音刺激的反應率是否有差別？8名受試對象對4種不同頻率聲音刺激的反應率（%）

受試號頻率A頻率B頻率C頻率D反應率秩反應率秩反應率秩反應率秩1

8.41

9.62

9.83

11.74211.6112.7411.8212.0339.429.1110.449.8349.828.719.9312.0458.328.018.63.58.63.568.619.839.6210.6478.919.0210.6311.4487.818.228.5310.84分析方法：若方差齊性，且各樣本均服從正態(tài)分布，用雙向方差分析（two-way-ANOVA)

若不滿足，用Friedman秩和檢驗若有統(tǒng)計學意義，進一步做多重比較例5某研究者為了比較甲、乙、丙、丁、戊、己6種藥物給家兔注射后產(chǎn)生的皮膚皰疹大?。╩m2），采用拉丁方設計，選用6只家兔、并在每只家兔的6個不同部位進行注射，試驗結(jié)果如下。分析方法拉丁方設計：一個處理因素，兩個區(qū)組因素，用three-way-ANOVA，不能分析交互作用家兔編號（行區(qū)組）注射部位編號（列區(qū)組）1234561C（87）B（75）E（81）D（75）A（84）F（66）2B（73）A（81）D（87）C（85）F（64）E（79）3F（73）E（73）B（74）A（78）D（73）C（77）4A（77）F（68）C（69）B（74）E（76）D（73）5D（64）C（64）F（72）E（76）B（70）A（81）6E（75）D（77）A（82）F（61）C（82）B（61）拉丁方設計與試驗結(jié)果（皮膚皰疹大小，mm2）

例6下表是A、B兩種閃爍液測定血漿中3H-cGMP的交叉試驗結(jié)果。第I階段1、3、4、7、9號用A測定，2、5、6、8、10號用B測定；第II階段1、3、4、7、9號用B測定，2、5、6、8、10號用A測定。試對交叉試驗結(jié)果進行方差分析。分析方法二階段交叉設計：處理因素，受試者，試驗階段，不分析交互作用。受試者階段ⅠⅡ1A（760）B（770）2B（860）A（855）3A（568）B（602）4A（780）B（800）5B（960）A（958）6B（940）A（952）7A（635）B（650）8B（440）A（450）9A（528）B（530）10B（800）A（803）

兩種閃爍液測定血漿中3H-cGMP的交叉試驗

例7觀察A,B兩種鎮(zhèn)痛藥物聯(lián)合運用在產(chǎn)婦分娩時的鎮(zhèn)痛效果。A藥取3個劑量：1mg，2.5mg，3mg；B藥也取3個劑量：5，15，30；共9個處理組。將27名產(chǎn)婦隨機等分9組，記錄分娩時的鎮(zhèn)痛時間，見表。試分析A,B兩藥聯(lián)合運用的鎮(zhèn)痛效果。

分析方法析因設計：兩個處理因素，各有3個水平，可分析交互作用。例8將手術要求基本相同的15名患者隨機分3組，在手術過程中分別采用A，B，C三種麻醉誘導方法，在T1（誘導前）、T2、T3、T3、T4，五個時相測量患者的收縮壓，數(shù)據(jù)記錄見表。試進行方差分析。

分析方法重復測量設計：一個處理因素有3個水平，一個重復因素有5個水平，可分析交互作用。（相應設計的協(xié)方差分析更合理）

方法序號麻醉誘導時相T0T1T2T3T4A1120108112120117A2118109115126123A3119112119124118A4121112119126120A5127121127133126B6121120118131137B7122121119129133B8128129126135142B9117115111123131B10118114116123133C11131119118135129C12129128121148132C13123123120143136C14123121116145126C15125124118142130

不同麻醉誘導時相患者的收縮壓（mmHg）

計數(shù)資料分析思路例9根據(jù)大量篩查得知某沿海省份成人的α-地中海貧血基因攜帶率為7.6%。某研究在本省某山區(qū)隨機抽樣調(diào)查成人125例，檢出α-地中海貧血基因攜帶者12例，問該山區(qū)與本省一般成人的α-地中海貧血基因攜帶率有無差異？分析方法：樣本率與總體率的比較：單樣本u檢驗（條件）例10根據(jù)以往經(jīng)驗，新生兒染色體異常率一般為1％，某醫(yī)院觀察了當?shù)?00名新生兒，只有1例染色體異常，問該地新生兒染色體異常是否低于一般。分析方法：樣本率與總體率的比較，發(fā)生數(shù)X較小，且總體率π接近于0或1的資料，用基于二項分布的確切概率法。例11某醫(yī)師研究用蘭芩口服液與銀黃口服液治療慢性咽炎療效有無差別，將病情相似的80名患者隨機分成兩組，分別用兩中藥物治療結(jié)果如下。兩種藥物療效比較藥物有效無效存合計有效率（%）蘭芩口服液41(36.56）4（8.44）4591.11銀黃口服液24(28.44)11(6.56)3568.57合計65158081.25分析方法：兩獨立樣本率比較，2檢驗，F(xiàn)isher確切概率法，兩樣本率u檢驗例12有260份血清樣品，每份樣品一分為二，用兩種不同的免疫學檢測方法檢驗類風濕因子，結(jié)果如下。問兩種檢測結(jié)果的陽性率是否相同。A法B法合計+-+1728180-126880合計18476260例13某抗癌新藥的毒理研究中，將78只大鼠按性別、窩別、體重、年齡等因素配成39對，每對大鼠經(jīng)隨機分配分別接受甲劑量和乙劑量注射，實驗結(jié)果見下表。試分析該新藥兩種不同劑量的毒性有無差異。

甲劑量乙劑量合計死亡(+)生存(-)死亡(+)6（a）12（b）18生存(-)3（c）18（d）21合計93039分析方法：配對設計率的比較：McNemar檢驗例14黑色素瘤患者隨訪術后三年生存情況

治療組三年生存例數(shù)三年死亡例數(shù)合計三年生存率（%）I7710818541.6II8910319246.4III9910420348.8IV909118150.0合計35540676146.6分析方法：多個樣本率（構(gòu)成比）的比較，2檢驗例15三種療法治療慢性喉炎的結(jié)果治療方法例數(shù)治愈顯效好轉(zhuǎn)無效合計綜合治療186722624308電子治療儀3224162092清音丸2214222078分析方法：單向有序行列表，反應變量有序，秩和檢驗

例16黑色素瘤患者隨訪術后三年生存情況

年齡組三年生存例數(shù)三年死亡例數(shù)合計三年生存率（%）30-7710818541.640-8910319246.450-9910420348.860-909118150.0合計35540676146.6分析方法：單向有序行列表，分組變量有序，視為多組樣本率的比較，用2檢驗例17對150名冠心病患者用兩種方法檢查室壁收縮運動的情況，檢查結(jié)果見表，試比較兩方法結(jié)果的概率分布有無差別。甲乙合計正常減弱異常正常603265減弱042951異常891734合計685428150分析方法：雙向有序行列表，配對設計（行列性質(zhì)一樣），用MarginalHomogeneity檢驗有無差別(McNemar法由二分類到多個有序分類的擴展方法）,或Kappa檢驗兩種方法的一致性。例18不同期次矽肺患者肺門密度級別分布，二者相關否？

矽肺期次肺門密度級別

合計＋＋＋＋＋＋Ⅰ4318814245Ⅱ19672169Ⅲ6175578合計50301141492分析方法：雙向有序行列表，單樣本設計（行列指標不同），用秩相關分析。例19不同年齡組英國男性醫(yī)生肺癌死亡率與吸煙與否資料分析方法：分層計數(shù)資料，Mantel-Haenszel分層2檢驗，Logistic

回歸分析。雙變量資料簡單相關分析線性回歸分析曲線回歸分析多變量資料例20糖尿病人的血糖可能受血清總膽固醇、甘油三酯、空腹胰島素、糖化血紅蛋白、的影響。27名患者的有關測量值如下，試建立血糖與其他幾項指標關系的多元線性回歸方程。

27名糖尿病人的血糖及有關變量的測量結(jié)果總膽固醇甘油三脂胰島素糖化血血糖（mmol/L）(mmol/L)(μU/ml)紅蛋白(%)(mmol/L)X1X2X3X4Y15.681.904.538.211.223.791.647.326.98.836.023.566.9510.812