秦秉杰
(山西藥科職業(yè)學(xué)院,山西 太原 030031)
藥物的體內(nèi)藥效學(xué)實(shí)驗(yàn)用以研究藥物對(duì)臨床疾病的功效以及藥物劑量對(duì)藥效的影響,是藥物臨床前評(píng)價(jià)的重要組成部分,也是藥物臨床藥理評(píng)價(jià)和指導(dǎo)臨床合理用藥的基礎(chǔ)[1]。醫(yī)用數(shù)理統(tǒng)計(jì)在實(shí)驗(yàn)數(shù)據(jù)的收集整理、統(tǒng)計(jì)分析與推斷解釋等方面發(fā)揮了重要的作用。SPSS(Statistical Package for the Social Science)統(tǒng)計(jì)分析軟件,可以更簡(jiǎn)單、快捷地完成數(shù)據(jù)的統(tǒng)計(jì)分析,將繁瑣的統(tǒng)計(jì)計(jì)算變得簡(jiǎn)單,同時(shí)使較短時(shí)間內(nèi)準(zhǔn)確完成大量的數(shù)據(jù)統(tǒng)計(jì)分析成為可能。SPSS在處方篩選[2]、藥效學(xué)研究[3-4]、毒性研究[5]等方面均有廣泛的應(yīng)用。下面應(yīng)用SPSS軟件對(duì)某藥物的體內(nèi)藥效學(xué)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。
研究某藥物的體內(nèi)藥效學(xué)。動(dòng)脈粥樣硬化小鼠尾靜脈注射低、中、高劑量的藥物或者生理鹽水,治療兩周。治療結(jié)束后,安樂死動(dòng)脈粥樣硬化小鼠,取小鼠主動(dòng)脈根并制備冰凍切片,采用油紅O對(duì)斑塊脂質(zhì)進(jìn)行染色后,定量統(tǒng)計(jì)動(dòng)脈粥樣硬化斑塊面積與血管腔面積的比值,收集實(shí)驗(yàn)數(shù)據(jù)如表1所示。通過研究各給藥劑量組的斑塊面積與血管腔面積的比值有無顯著性差異,判斷藥物的抗動(dòng)脈粥樣硬化作用及其抗動(dòng)脈粥樣硬化作用與給藥劑量之間的關(guān)系。
表1 各給藥劑量組斑塊面積與血管腔面積的比值
研究各給藥劑量組的斑塊面積與血管腔面積的比值有無顯著性差異,即驗(yàn)證各給藥劑量組的斑塊面積與血管腔面積的比值對(duì)應(yīng)總體的分布是否有顯著性差異。本實(shí)驗(yàn)含有四組相互獨(dú)立的實(shí)驗(yàn)數(shù)據(jù),為多獨(dú)立樣本。同時(shí)由于各組的總體分布未知,樣本量較少,故采用多獨(dú)立樣本的非參數(shù)檢驗(yàn)方法。多獨(dú)立樣本的非參數(shù)檢驗(yàn)方法包括:一是Kruskal-Willis H檢驗(yàn),重點(diǎn)考察各組秩的均值;二是Median檢驗(yàn),檢驗(yàn)各組的中位數(shù);三是Joneckheere-Terpstra檢驗(yàn),可用于檢驗(yàn)各組之間是否有順序效應(yīng)。設(shè)置零假設(shè)H0:各給藥劑量組的斑塊面積與血管腔面積的比值對(duì)應(yīng)總體的分布無顯著性差異。
首先采用SPSS統(tǒng)計(jì)軟件設(shè)置變量“組別”與“斑塊面積與血管腔面積比值”,錄入實(shí)驗(yàn)數(shù)據(jù)如圖1所示。在SPSS軟件數(shù)據(jù)視圖主菜單中,選擇分析→非參數(shù)檢驗(yàn)→舊對(duì)話框→K個(gè)獨(dú)立樣本進(jìn)行多個(gè)獨(dú)立樣本的非參數(shù)檢驗(yàn)如圖2所示。設(shè)定分組變量,單擊定義范圍,設(shè)定分組變量的最小值與最大值如圖3所示。選擇檢驗(yàn)方法進(jìn)行多獨(dú)立樣本的非參數(shù)檢驗(yàn)如圖4所示。若各給藥劑量組的斑塊面積與血管腔面積的比值對(duì)應(yīng)總體的分布存在顯著性差異,可以進(jìn)一步地進(jìn)行多個(gè)獨(dú)立樣本間的兩兩比較的非參數(shù)檢驗(yàn),進(jìn)一步地檢驗(yàn)?zāi)膬山M的總體分布存在顯著性的差異。采用兩個(gè)獨(dú)立樣本的非參數(shù)檢驗(yàn),Mann-Whitney U檢驗(yàn),檢驗(yàn)四組樣本間兩兩之間的顯著性差異。
圖1 各給藥劑量組斑塊面積與血管腔面積的比值數(shù)據(jù)集
圖2 多個(gè)獨(dú)立樣本的非參數(shù)檢驗(yàn)
圖3 多個(gè)獨(dú)立樣本檢驗(yàn)(定義范圍)
圖4 選擇檢驗(yàn)類型
表2 秩統(tǒng)計(jì)表
表3 Kruskal-Wallis H檢驗(yàn)統(tǒng)計(jì)量
由表2可知,各治療組的斑塊面積與血管腔面積比值的平均秩依次為16.00、13.00、7.80、2.50。由表3 Kruskal-Wallis H檢驗(yàn)統(tǒng)計(jì)量可知,卡方為14.434,概率P值為0.002。對(duì)顯著性水平α=0.05,概率P值0.002<0.05,拒絕零假設(shè),表明各治療組的斑塊面積與血管腔面積比值的平均秩存在顯著性差異,各治療組的斑塊面積與血管腔面積比值對(duì)應(yīng)的總體分布存在顯著性差異。
表4 頻數(shù)表
表5 Median檢驗(yàn)統(tǒng)計(jì)量
a. 8個(gè)單元格 (100.0%) 的期望頻數(shù)小于5,最小的期望頻數(shù)為 2.0。
由表4與表5可知,四組共同的中位數(shù)為29.82,卡方為10.267,概率P值為0.016。對(duì)顯著性水平α=0.05,概率P值0.016<0.05,拒絕零假設(shè),表明各治療組的斑塊面積與血管腔面積比值對(duì)應(yīng)的總體分布存在顯著性差異??紤]到卡方檢驗(yàn)對(duì)期望頻數(shù)小于5的單元格數(shù)的限制以及SPSS給出的提示,該檢驗(yàn)只能作為參考。
表6 Joneckheere-Terpstra檢驗(yàn)表
由表6可知,Joneckheere-Terpstra統(tǒng)計(jì)量的概率P值為0.000。對(duì)顯著性水平α=0.05,概率P值0.000<0.05,拒絕零假設(shè),表明各治療組的斑塊面積與血管腔面積比值對(duì)應(yīng)的總體分布存在顯著性差異,且各治療組的斑塊面積與血管腔面積比值存在順序效應(yīng),表明隨著給藥劑量的增加,藥物的抗動(dòng)脈粥樣硬化療效增加。
在得到上述結(jié)果后,對(duì)多個(gè)獨(dú)立樣本間進(jìn)行兩兩比較的非參數(shù)檢驗(yàn),實(shí)驗(yàn)結(jié)果如下:
表7 生理鹽水VS低劑量組
表8 生理鹽水VS中劑量組
表9 生理鹽水VS高劑量組
表10 低劑量組VS中劑量組
表11 低劑量組VS高劑量組
表12 中劑量組VS高劑量組
由表7可知,精確概率P值0.177,對(duì)顯著性水平α=0.05,精確概率P值0.177>0.05,接受零假設(shè),表明生理鹽水組與低劑量組的斑塊面積與血管腔面積比值對(duì)應(yīng)的總體分布沒有顯著性差異。由表8—表12可知,對(duì)顯著性水平α=0.05,精確概率P值均<0.05,拒絕零假設(shè),表明生理鹽水組與中劑量組、高劑量組;低劑量組與中劑量組、高劑量組;中劑量組與高劑量組的斑塊面積與血管腔面積比值對(duì)應(yīng)的總體分布均存在顯著性差異,表明藥物具有良好的抗動(dòng)脈粥樣硬化作用,且其抗動(dòng)脈粥樣硬化作用與給藥劑量呈正相關(guān)。
應(yīng)用SPSS軟件更簡(jiǎn)單、快捷的完成對(duì)某藥物的體內(nèi)藥效學(xué)數(shù)據(jù)的統(tǒng)計(jì)分析。采用多獨(dú)立樣本的非參數(shù)檢驗(yàn)方法,進(jìn)行了Kruskal-Willis H檢驗(yàn)、Median檢驗(yàn)和Joneckheere-Terpstra檢驗(yàn),經(jīng)檢驗(yàn)認(rèn)為各治療組的斑塊面積與血管腔面積比值對(duì)應(yīng)的總體分布存在顯著性差異。同時(shí),對(duì)多個(gè)獨(dú)立樣本間進(jìn)行兩兩比較的非參數(shù)檢驗(yàn)和Mann-Whitney U檢驗(yàn),經(jīng)檢驗(yàn)認(rèn)為生理鹽水組與低劑量組的斑塊面積與血管腔面積比值對(duì)應(yīng)的總體分布沒有顯著性差異,而生理鹽水組與中劑量組、高劑量組;低劑量組與中劑量組、高劑量組;中劑量組與高劑量組的斑塊面積與血管腔面積比值對(duì)應(yīng)的總體分布均存在顯著性差異。統(tǒng)計(jì)分析結(jié)果表明,該藥物具有良好的抗動(dòng)脈粥樣硬化作用,且其抗動(dòng)脈粥樣硬化作用與給藥劑量呈正相關(guān)。