參照單位分析適用于等級資料,如按療效分為治愈、好轉(zhuǎn)、無效、惡化,按反應分為-、+、++、+++,和按麻醉效果分為Ⅰ、Ⅱ、Ⅲ、Ⅳ級等的資料。
參照單位分析又稱Ridit分析,常是將等級資料中例數(shù)較多的一組的分布作為一個特定的分布來計算各等級的參照單位值(R值),再參照這些R值計算各組的加權(quán)平均R值并進行假設檢驗。Ridit的前三個字母是Relative to an indentified distribution的縮寫,-it是unit的字尾,故其含義是“相對于某一特定分布的單位”。由于有時將idantified distribution稱為reference distribution,故譯為參照單位 。
參照單位分析一般計算步驟如下:
1.選標準組
標準組的選擇可根據(jù)各組例數(shù)多少以及所研究的問題而定。一般選例數(shù)多的組為標準組。如果各組例數(shù)相近或都較少時,可用合計數(shù)為標準組,(但若各組例數(shù)過少,則不宜用此法,可用秩和檢驗);若研究的是新、舊藥物的療效,則可以舊藥為標準組;若研究的是患者與正常人相比較,則可選正常人為標準組。另外,標準組中的數(shù)字要求分布于各個等級,如果有的等級為0或過于少,將對計算結(jié)果產(chǎn)生影響。本例選辯中西組為標準組。
2.計算標準組的參照單位值,(R值)
計算前最好將各等級按由弱到強的次序排列。計算的步驟與方法是:①計算標準組各等級的1/2值②求標準組累計例數(shù)并下移一行③將①、②求得的值按各等級相加④以標準組總例數(shù)除之,即得標準組各等級的R值。
3.參照標準組的R值,計算各組的平均R值(r )
將每組中各不同療效者例數(shù)與標準組對應的R值相乘,將乘積加總,再除以該組總例數(shù),即為R,公式為
(10.12)
4.計算各組的95%可信區(qū)間進行顯著性檢驗,R值的標準差為,故R值的標準誤為
。10.13)
R的95%可信限為
(10.14)
按此公式算得各組的95%可信區(qū)間后,作兩兩間的比較,凡無重疊者即有顯著差別;有重疊者則無顯著差別。
例10.9 四種療法治療潰瘍病的療效如表10.11,試比較其療效。
表10.11 幾種療法治療潰瘍病的療效
療效 |
例 數(shù) |
百分比(%) | ||||||
西組 | 中組 | 辯中組 | 辯中西組 | 西組 | 中組 | 辯中組 | 辯中西組 | |
治愈 | 15 | 10 | 16 | 214 | 41.7 | 18.5 | 40.0 | 61.1 |
顯效 | 8 | 11 | 13 | 93 | 22.2 | 20.4 | 32.5 | 26.6 |
好轉(zhuǎn) | 9 | 20 | 9 | 39 | 25.0 | 37.0 | 22.5 | 11.2 |
無效 | 4 | 13 | 2 | 4 | 11.1 | 24.1 | 5.0 | 1.1 |
合計 | 36 | 54 | 40 | 350 | 100.0 | 100.0 | 100.0 | 100.0 |
注:西組為單純西藥組,中組為單純中藥組,辯中組為辯證分型中草藥組,辯中西組為辯證分型中西醫(yī)結(jié)合組。
選辯中西組為標準組,并計算其R值如表10.12
表10.12 標準組R值計算表
療效 | 例數(shù)(1) | (1)/2 (2) |
累計例數(shù)并下移一行(3) | (2)+(3) (4) |
R值(4)/n (5) |
無效 | 4 | 2.0 | 0 | 2.0 | 0.006 |
好轉(zhuǎn) | 39 | 19.5 | 4 | 23.5 | 0.067 |
顯效 | 93 | 46.5 | 43 | 89.5 | 0.256 |
治愈 | 214 | 107.0 | 136 | 243.0 | 0.694 |
N=350
R值的計算,也可用百分比進行,結(jié)果相同,見表10.13
表10.13 用百分比求R值的計算表
療效 |
百分比(1) |
(1)/2 |
累計例數(shù)并下移一行(3) |
(2)+(3) |
R值(4)/n |
無效 | 1.1 | 0.55 | 0 | 0.55 | 0.006 |
好轉(zhuǎn) | 11.2 | 5.60 | 1.1 | 6.70 | 0.067 |
顯效 | 26.6 | 13.30 | 12.3 | 25.60 | 0.256 |
治愈 | 61.6 | 30.55 | 38.9 | 69.45 | 0.695 |
N=100.0
再計算各組的R值(加權(quán)平均參照單位值),用式(10.12),式中各符號結(jié)合本例含義是
式中∑為總和,f為各不同療效者例數(shù),R為與f同一療效之標準組的R值。
標準組的R值應為0.5,可用以核對R值是否算對。
本題各組R計算如下:
R辯中西=(4×0.006+39×0.067+93×0.256+2.4×0.694)/350=0.500
R西=(4×0.006+9×0.067+8×0.256+15×0.694)/36=0.363
R中=(13×0.006+20×0.067+11×0.256+10×0.694)/54=0.207
R辯中 =(2×0.006+9×0.067+13×0.256+16×0.694)/40=0.376
將各組R值與例數(shù)代入式10.14,得各組95%可信限如下:
辯中西
西組
中組
辯中組
將各組的95%可信區(qū)間圖示如下:
圖10.1 四種療法治療潰瘍病療效的比較(Ridit分析)
根據(jù)以上四個區(qū)間及其示意圖可以看出,辯中西組的區(qū)間與其它三組都不重疊,其它三組中西組和辯中組兩區(qū)間有較大重疊,而與中組都是稍有重疊,在這種情況下可作較精確的分析,因為上述R值的標準差為,或方差為1/12,只是一個近似值,據(jù)數(shù)理統(tǒng)計研究,R值的方差是一個逐漸接近1/12并以1/12為最大值。它與等級數(shù)有關(guān),隨著等級的增多,最大方差愈來愈接近1/12。表10.14是資料的等級數(shù)與最大方差。
表10.14 各種等級數(shù)時的最大方差
等級數(shù) | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
最大方差 | 1/16 | 2/27 | 5/64 | 2/25 | 30/432 | 4/49 | 21/256 | 20/243 |
實際資料的方差都小于1/12,故用以推算得的可信區(qū)間偏大,假設檢驗結(jié)論比較保守。本例中組與辯中組的區(qū)間重疊部分甚少,用表10.14內(nèi)數(shù)值代替1/12,則可得到更為適當?shù)慕Y(jié)果。該資料等級數(shù)為4,查表10.14得最大方差為5/64故標準誤為,以此計算西組、中組、辯中組的95%可信區(qū)間如下:
西組
中組
辯中組
這樣就可以認為中組與辯中組療效相差也是顯著的,雖然這里用的仍是最大方差,但是等級為4時的最大方差,假設檢驗結(jié)論還是比較保守的,但比用1/12要精確些,當?shù)燃壭∮?時,更以查此表數(shù)值為宜。