利用統(tǒng)計(jì)量作隨機(jī)區(qū)組設(shè)計(jì)資料處理的簡(jiǎn)便方法

摘要:目的:探索一種利用描述統(tǒng)計(jì)量直接計(jì)算隨機(jī)區(qū)組設(shè)計(jì)資料方差分析結(jié)果的簡(jiǎn)便方法。方法:利用數(shù)理統(tǒng)計(jì)的方法推導(dǎo)出描述統(tǒng)計(jì)量與隨機(jī)區(qū)組設(shè)計(jì)的方差分析計(jì)算公式間的關(guān)系，從而找到一種利用描述統(tǒng)計(jì)量直接計(jì)算隨機(jī)區(qū)組設(shè)計(jì)的方差分析檢驗(yàn)統(tǒng)計(jì)量的方法。結(jié)論:通過(guò)實(shí)例運(yùn)算與用原始數(shù)據(jù)算得的結(jié)果一樣，本文提供的在只有描述統(tǒng)計(jì)量的情況下仍然可以作隨機(jī)區(qū)組設(shè)計(jì)的方差分析的方法的是一種簡(jiǎn)便有效的方法。

關(guān)鍵詞:描述統(tǒng)計(jì)量;隨機(jī)區(qū)組設(shè)計(jì)的方差分析;離差平方和

隨機(jī)區(qū)組設(shè)計(jì)的方差分析由于它能有效控制實(shí)驗(yàn)過(guò)程中的混雜因素，而且只用較少的樣本就可得到較理想的結(jié)果，因而在各種科學(xué)實(shí)驗(yàn)，特別是醫(yī)藥衛(wèi)生科學(xué)研究中被廣泛采用，在醫(yī)藥衛(wèi)生統(tǒng)計(jì)中又被稱為配伍組設(shè)計(jì)，盡管該方法效率很高，樣本數(shù)量也不大，實(shí)驗(yàn)數(shù)據(jù)不是很多，但它的計(jì)算仍很麻煩，目前采用計(jì)算機(jī)軟件作數(shù)據(jù)處理的做法比較普遍，各種統(tǒng)計(jì)軟件以其權(quán)威性、規(guī)范性、準(zhǔn)確性已成為當(dāng)前數(shù)據(jù)處理的主流做法，但由于大多數(shù)統(tǒng)計(jì)軟件都是利用原始數(shù)據(jù)作各種統(tǒng)計(jì)分析，很少直接利用已算出的樣本統(tǒng)計(jì)量去作統(tǒng)計(jì)分析。而在實(shí)際工作中，有許多情況是研究者只保留了作描述的各種統(tǒng)計(jì)量，如均數(shù)、標(biāo)準(zhǔn)差等，當(dāng)未保存原始數(shù)據(jù)而又要作進(jìn)一步的統(tǒng)計(jì)分析時(shí)不僅各種統(tǒng)計(jì)軟件不能直接處理，就是手工計(jì)算或利用計(jì)算器處理也找不到適當(dāng)?shù)墓剑@方面已有雜志介紹過(guò)利用計(jì)量資料的描述統(tǒng)計(jì)量作顯著性檢驗(yàn)的方法，但僅對(duì)簡(jiǎn)單t檢驗(yàn)，單因素方差分析可以解決，而在實(shí)際數(shù)據(jù)處理中，特別是醫(yī)藥衛(wèi)生科學(xué)研究中有不少隨機(jī)區(qū)組設(shè)計(jì)的資料，對(duì)它們分析要用到的方差分析，但找不到相應(yīng)的利用描述統(tǒng)計(jì)量作計(jì)算的方法。我們通過(guò)建立描述統(tǒng)計(jì)量與隨機(jī)區(qū)組設(shè)計(jì)的方差分析計(jì)算公式之間的關(guān)系，從而解決在只有描述統(tǒng)計(jì)量時(shí)計(jì)算這類資料方差分析的一種簡(jiǎn)便方法。

1. 方法

1. 1確定描述統(tǒng)計(jì)量的基本數(shù)據(jù)。

通常隨機(jī)區(qū)組設(shè)計(jì)資料的描述統(tǒng)計(jì)量有如表1所示:1. 2利用數(shù)理統(tǒng)計(jì)的定義推導(dǎo)出由描述統(tǒng)計(jì)量計(jì)算隨機(jī)區(qū)組設(shè)計(jì)資料的方差分析檢驗(yàn)統(tǒng)計(jì)量的公式。

由表1提供的統(tǒng)計(jì)量利用數(shù)理統(tǒng)計(jì)的定義，可導(dǎo)出各因素的離差平方和，列于表2。

1. 3列出隨機(jī)區(qū)組設(shè)計(jì)資料的方差分析表。見表3。

由自由度(a-1，N-a-b+2)和(b-1，N-a-b+2)以及事先定好的檢驗(yàn)水準(zhǔn)?琢可查得兩個(gè)F界值:F?琢(a-1，N-a-b+2)與F?琢(b-1，N-a-b+2)，分別與兩個(gè)F統(tǒng)計(jì)量比較便可確定P值了。

2. 實(shí)例計(jì)算

為探索丹參對(duì)肢體缺血再灌注損傷的影響，將30只純種新西蘭實(shí)驗(yàn)用大白兔，按窩別相同、體重相近分為10個(gè)區(qū)組。每個(gè)區(qū)組3只大白兔隨機(jī)采用A、B、C三種處理方案，即在松止血帶前分別給予丹參2ml/kg、丹參1ml/kg、生理鹽水2ml/kg，在松止血帶前及松后1小時(shí)分別測(cè)定血中白蛋白含量(g/L)算出白蛋白減少量，問(wèn)A、B兩方案分別與C方案的處理效果是否不同?

再算出實(shí)驗(yàn)組A與區(qū)組B各平均數(shù)的標(biāo)準(zhǔn)差，然后用公式SSA=b(a-1)S2■A，SSB=a(b-1)S2■B具體計(jì)算見表5。

最后列出方差分析表，見表6。

由表6的結(jié)果作出推斷結(jié)論: 按α=0.05的檢驗(yàn)水準(zhǔn)，實(shí)驗(yàn)組間差異有統(tǒng)計(jì)學(xué)意義，而區(qū)組差異無(wú)統(tǒng)計(jì)學(xué)意義。這與用原始數(shù)據(jù)計(jì)算的結(jié)果完全一樣。

3. 結(jié)論

在無(wú)原始數(shù)據(jù)的情況下，我們利用數(shù)理統(tǒng)計(jì)的基本定義推導(dǎo)出了只用描述統(tǒng)計(jì)量求解隨機(jī)區(qū)組設(shè)計(jì)資料作方差分析的計(jì)算公式，并通過(guò)實(shí)例作了驗(yàn)證。這種簡(jiǎn)便方法在醫(yī)藥工作中有很大的實(shí)用性，雖然現(xiàn)在利用計(jì)算機(jī)采用統(tǒng)計(jì)軟件作數(shù)據(jù)處理已是一種常規(guī)方法，但對(duì)非統(tǒng)計(jì)專業(yè)的實(shí)際工作者來(lái)說(shuō)并非都能作到，他們期望有一種較簡(jiǎn)便的方法用描述統(tǒng)計(jì)量可直接作隨機(jī)區(qū)組設(shè)計(jì)資料的方差分析。另一方面，對(duì)確實(shí)只保留了數(shù)據(jù)的描述統(tǒng)計(jì)量，而無(wú)原始數(shù)據(jù)的實(shí)際工作者來(lái)說(shuō)，上述提供的方法對(duì)作隨機(jī)區(qū)組設(shè)計(jì)資料的方差分析，無(wú)疑是一種很好的彌補(bǔ)辦法。上述方法我們?cè)诮虒W(xué)實(shí)踐中應(yīng)用也獲得了較好的效果，本方法對(duì)無(wú)重復(fù)數(shù)據(jù)的雙向(因素)方差分析(two-way ANOVA)都是可行的。但對(duì)有交互作用或考慮交互作用的雙因素方差分析來(lái)說(shuō)，本方法就有一定的局限性，此時(shí)只有利用統(tǒng)計(jì)軟件來(lái)處理還是適宜的。

參考文獻(xiàn):

[1]高祖新主編.中國(guó)醫(yī)藥科技出版社，2008.6.

[2]祝國(guó)強(qiáng).醫(yī)藥數(shù)理統(tǒng)計(jì)方法[M].高等教育出版社，2004.

[3]方積乾主編.衛(wèi)生統(tǒng)計(jì)學(xué)[M].北京:人民衛(wèi)生出版社，2003.

責(zé)任編輯 羅峰