CHT統(tǒng)計(jì)軟件及其在藥效綜合分析中應(yīng)用效果的評(píng)估

1、    CHT統(tǒng)計(jì)軟件及其在藥效 綜合分析中應(yīng)用效果的評(píng)估        臨床上比較兩種療法或兩種藥物的療效，常從不同角度對(duì)它們進(jìn)行考察，此時(shí)對(duì)比的指標(biāo)往往不止一個(gè)。以往由于種種原因，都采用單因素分析法，不但計(jì)算量大，而且當(dāng)各指標(biāo)的表現(xiàn)趨勢(shì)不一致時(shí)，就難以作出明確的判斷，因此許多學(xué)者采用了Hotelling T2檢驗(yàn)這一多元統(tǒng)計(jì)推斷方法，收到了事半功倍之效。近年來，Hotelling T2檢驗(yàn)的應(yīng)用性研究較多，但作者根據(jù)公式推導(dǎo)和實(shí)例分析發(fā)現(xiàn)，Hotelling T

2、2檢驗(yàn)在某些情況下可能出現(xiàn)假陽性的統(tǒng)計(jì)推斷結(jié)果1，因而在隨后的研究中提出了相應(yīng)的校正方法2,3，并以校正的Hotelling T2檢驗(yàn)作為數(shù)理統(tǒng)計(jì)基礎(chǔ)設(shè)計(jì)微機(jī)程序(另文發(fā)表)，并與該程序相關(guān)的一些文件一起組成CHT統(tǒng)計(jì)軟件。本文采用CHT軟件對(duì)藥效研究中的多變量資料進(jìn)行綜合分析，并與原Hotelling T2檢驗(yàn)進(jìn)行比較，在單變量分析的基礎(chǔ)上分析藥物的綜合藥效，客觀地評(píng)價(jià)CHT軟件的應(yīng)用效果。 1CHT軟件簡(jiǎn)價(jià)1.1數(shù)據(jù)輸入組數(shù)G、變量數(shù)P和每組數(shù)本數(shù)n(i)為鍵盤輸入，其它原始數(shù)據(jù)則根據(jù)使用者的使用目的確定輸入方式，即：隨機(jī)模擬輸入鍵入零，READDATA形式讀入鍵入1，文本文件讀入鍵入3。

3、1.2結(jié)果輸出首先屏幕提示，顯示或打印輸出鍵入1為打印輸出，否則為屏幕顯示，在原始數(shù)據(jù)列表輸出后，接著屏幕詢問：資料來源?操作者可根據(jù)需要，或鍵入資料來源，或直接按回車鍵。接著列表輸出有關(guān)信息，如(1)均差向量：M表示組號(hào)較小的均數(shù)向量減去組號(hào)較大的均數(shù)向量，-M則相反;(2)變量;(3)某變量的均差小于該變量的極小值點(diǎn)均差;(4)極小值點(diǎn)T2值、未校正T2值和校正后的T2值。最后在表下方輸出F值、自由度dF1、dF2、P值。1.3數(shù)據(jù)庫(kù)文件和文本文件的建立采用foxbase 2.0+，以“序號(hào)”及各變量為字段名，生成數(shù)據(jù)庫(kù)結(jié)構(gòu)后，每一個(gè)體給以一定序號(hào)后作為一條記錄輸入，組別以序號(hào)區(qū)別。數(shù)據(jù)輸

4、入完畢后，鍵入quit退出foxbase 2.0+系統(tǒng)。文本文件的生成方法見文獻(xiàn)4。1.4CHT軟件的使用方法在ucdos下鍵入QBASIC，鼠標(biāo)選擇：打開文件裝載CHT.BAS運(yùn)行。此時(shí)屏幕依次出現(xiàn)：組數(shù)，變量數(shù)?鍵入G，P值；每組樣本數(shù)?鍵入n(i)i=1,鍵入n(1)；i=2,鍵入n(1),n(2);余類推。選擇數(shù)據(jù)輸入方式隨機(jī)鍵入0，READDATA鍵入1，文件讀入鍵入2，若為鍵入1方式，則必須事先設(shè)置DATA中的數(shù)據(jù)；若為鍵入2方式，則先建立數(shù)據(jù)庫(kù)并生成相應(yīng)的文本文件。并按提示輸入文本文件名和輸出運(yùn)行結(jié)果。2CHT軟件在藥效綜合分析中應(yīng)用效果的評(píng)估本文采用隨機(jī)模擬和實(shí)例分析相結(jié)合的方

5、法評(píng)價(jià)CHT軟件在藥效分析中的應(yīng)用效果。2.1CHT軟件的隨機(jī)模擬運(yùn)行結(jié)果分析根據(jù)CHT軟件的模擬運(yùn)行方式，共進(jìn)行1170次的統(tǒng)計(jì)推斷，結(jié)果統(tǒng)計(jì)量T2需校正者614次，占52.4%，其中29次的校正值按C=2T2(min)計(jì)算，占4.6%。2.2CHT軟件的實(shí)例運(yùn)行結(jié)果分析本文共收集23例分散于有關(guān)刊物上的原始資料，采用CHT軟件處理，在23例實(shí)例分析中，共進(jìn)行35次統(tǒng)計(jì)推斷(其中6例是三組間兩兩比較)，結(jié)果統(tǒng)計(jì)量T2應(yīng)作校正者23次，占65.7%，其中3次的校正值按C=2T2(min)計(jì)算，占13%，均高于模擬結(jié)果。這可能與應(yīng)用例數(shù)少有關(guān)，也可能與實(shí)例的數(shù)據(jù)相對(duì)復(fù)雜些有關(guān)。2.3典型實(shí)例分析

6、例15用家兔8只(組1)，測(cè)定新藥利明降壓片的降壓效果，測(cè)定項(xiàng)目有降壓開始時(shí)間(變量1)、達(dá)降壓高峰時(shí)間(變量2)、最大降壓值(變量3)、最大降壓持續(xù)時(shí)間(變量4)，并以利血平作對(duì)照組(組2)，用家兔9只，試在單因素分析的基礎(chǔ)上用CHT軟件綜合分析利明降壓片的降壓效果。表1原始數(shù)據(jù)組1變量17.8 5.1 6.2 7.3 8.5 6.9 9.1 6.8組1變量288.3 121.5 89.7 98.5 109.8 113.5 97.9 127.1組1變量337.3 35.6 41.5 39.7 36.6 41.9 42.1 38.1組1變量44.3 5.5 5.7 7.1 6.7 4.8 4.

7、9 5.1組2變量178.7 53.4 85.8 91.7 88.9 62.3 77.4 66.3 58.7組2變量2198.3 201.5 198.2 231.3 211.7 171.9 191.2 221.5 219.4組2變量316.1 17.1 15.9 16.3 16.9 15.8 16 16.1 16.4組2變量46.8 6.4 7.1 7.5 8.1 7.7 7.2 6.9 6.7注：資料來源出自文獻(xiàn)5表2結(jié)果輸出均數(shù)向量變量均差<極值點(diǎn)均差T2(min)M1-66.47639-10.29292763.44722-99.2125-59.62257884.8593MIN T2

8、(min)=763.4472-M3-22.81111-1.819206282.263841.6430562.679425902.549MIN T2(min)=282.2638未校正T2值=911.4947校正后T2值615.3997F(1,2)=123.0799dF1=4dF2=12P(1 2)<0.0001 初步分析原始資料(表1)發(fā)現(xiàn)，對(duì)于變量1、變量2，組2優(yōu)于組1(P<0.0001)；變量3則組1優(yōu)于組2(P<0.05)；變量4則組間差異無統(tǒng)計(jì)學(xué)意義。CHT軟件統(tǒng)計(jì)結(jié)果見表2，綜合降壓效果，組2優(yōu)于組1(P<0.0001)，較之原Hotelling T2檢驗(yàn)，統(tǒng)

9、計(jì)量T2減少296.095，相當(dāng)于未校正T2值的30.8%。例2以膽固醇(變量1)、甘油三酯(變量2)和脂蛋白(變量3)為血脂高低的指標(biāo)，給血脂平以后12例患者的血脂變化值見表3，試用CHT軟件綜合分析降脂效果。表3原始數(shù)據(jù)組1變量117 103 50 70 -25 34 81 41 77 111 9 38組1變量2-52 20 156 -10 -10 -56 56 27 48 -27 -71 -55組1變量3-100 -25 200 -50 -175 -75 -100 0 25 85 -200 -75注：資料來源出自文獻(xiàn)6初步分析原始資料發(fā)現(xiàn)，變量1變化值有極顯著意義(P<0.001)

10、，變量2、變量3則無統(tǒng)計(jì)學(xué)意義(P>0.05)。CHT軟件綜合分析結(jié)果見表4，12例患者給血脂平后綜合降脂作用無統(tǒng)計(jì)學(xué)意義(P=0.3969)；經(jīng)原Hotelling T2檢驗(yàn)P<0.01。統(tǒng)計(jì)量T2減少42.0594，占未校正T2值的91.22%，使統(tǒng)計(jì)結(jié)論相反。表4結(jié)果輸出均數(shù)向量變量均差<極值點(diǎn)均差T2(min)M3-40.8333355.564221.02969MIN T2(min)=21.02969-M1-50.510.453553.7585022-2.16666729.0277440.36709MIN T2(min)=3.758502未校正T2值=46.10709

11、校正后T2值4.047714F=1.103922dF1=3dF2=9P=0.3969 例315例免疫功能低下患者，用胸腺素后，lgA(變量1)、lgG(變量2)和lgM(變量3)的變化值見表5，試在單變量分析的基礎(chǔ)上用CHT軟件綜合分析用藥效果。初步分析原始資料發(fā)現(xiàn)，給藥前后變量1的變化值有極顯著的統(tǒng)計(jì)學(xué)意義(P<0.001)，變量3 P<0.05,而變量2 P>0.05。經(jīng)CHT軟件分析，15例患者用胸腺素后，免疫功能改善(P=0.0004)，與原Hotelling T2檢驗(yàn)比較，統(tǒng)計(jì)量T2僅減少0.144%(見表6)。表5原始數(shù)據(jù)組1變量115617663128 -714

12、21041954202362141397115649組1變量2502081270-1562-74-11546088-11-9315組1變量3491421-918-19244-16382222-11420注：資料來源出自文獻(xiàn)6 表6結(jié)果輸出均數(shù)向量變量均差<極值點(diǎn)均差T2(min)M234.6666737.4907346.90632MIN T2(min)=46.90632-M1-149.4667-9.00052810.677583-12.066675.92451833.86543MIN T2(min)=10.67758未校正T2值=46.94014校正后T2值46.8725F=13.392

13、14dF1=3dF2=12P=0.00042.4CHT軟件應(yīng)用效果的評(píng)價(jià)CHT軟件具有全中文提示，操作簡(jiǎn)單，運(yùn)算快速，中文列表輸出等特點(diǎn)。通過模擬運(yùn)行和實(shí)例應(yīng)用表明，CHT軟件的應(yīng)用可使三分之二實(shí)例的統(tǒng)計(jì)量T2得到校正，可在單變量分析的基礎(chǔ)上使藥效的綜合評(píng)價(jià)不致于出現(xiàn)假陽性的統(tǒng)計(jì)推斷結(jié)果。這一事實(shí)充分表明，對(duì)Hotelling T2檢驗(yàn)進(jìn)行校正是必要的。文獻(xiàn)23根據(jù)原Hotelling T2檢驗(yàn)的偏差程度作相應(yīng)的校正是合理的，也是切實(shí)可行的。作者單位：蘇銀法(溫州市醫(yī)藥科學(xué)研究所溫州325003)參考文獻(xiàn)1蘇銀法，杜樂燕.極值判斷法在Hotelling T2檢驗(yàn)處理多變量資料中的應(yīng)用.數(shù)理醫(yī)藥學(xué)雜志，1995，8(2)：160.2蘇銀法，杜樂燕.Hotelling T2檢驗(yàn)分段函數(shù)的建立及應(yīng)用.數(shù)理醫(yī)藥學(xué)雜志，1996，9(3)：235.3蘇銀法.Hotelling T2檢驗(yàn)校正值計(jì)算的兩條基本原則.數(shù)理醫(yī)藥學(xué)雜志(待發(fā)表).4