样本率多重比较方法的模拟研究
中国卫生 统计 20 0年 94 月第 2 卷6第 2
期样
本率 重多 较比方 的模 拟 研 究
法广州南 医方大科公共学卫学 生生院统物 系计 ( 11 ) 钱俊 陈 平雁 5 05 5
【提 】 要的 目用 采Mn ao o t Cr方法模拟研 样本究多率 比较方 重特法 性方。 法e l不 同算 数参条 下件种各 本样率多 重较比方法总的I型错误 率 检、验 效和判对 能。率结果
用 S
S3.A 9 1软 编件 程, 拟 模 不估同 方法具各 特色 ,能较 控
好制
总I 型误错率且检效 验能较高方的法 R 有 a、me 、o t H ̄ e 、o hS p 、 t方法。结B论如 严需格控 总制 ny P r Bo o lH 和 、ctS S pe e I 型错率误 在水 准,内o Hmm、 oehS p、 B效能较t 高; l c 、tH SS ee p如 许总 I允 型错误 率溢略 出 (.5L 11 o内) R以a e 、n和 ,y nP
B no法能较效 。高 t o
【
关 键 】词多 重比较
总I
错误率型
检
验效能 对判
率多
重比较是样 本率 数 据 理 中常处 见的问 题 ,本样
评 指价标 要 主包 括 FR 、效验 能 对、 率及 以 E 检 判W稳
健性 。
率 多 重比较 方 的法有 近2 0 种之 … 多 见 , 方法有 常B ernn 法 、N法 、Kef 等c¨。 在 众方 多 of 法i S or Sfh e 法中, 法方 用运 简便 、哪 些果结稳 健 可 ?信它们 结的 各 果具 什么特点 ?本 文 助计 算 借机 拟 , 模索不 同方 法 探的优 劣性
, 讨论在不 同数 据 类型下 些哪方 法 值 得 更实际
用。 应 方
3 法模 拟数据 . 在多 重 较 中 比,组 数为 k 多 重较 比次的数 k 总为, ( 1 一 k / )次 , 2拟 数据要 考虑 主以下 三 种情 :形 模
() 个 多体率 全 部总相等 , 各 察 方控法 制 F 1考 w— E的能 力 。能否 将F R 的 值 控 制在 检 验 水准 RWE , 这是内 评 多价 比较重方法 的一 重要 项指 标J
。
用 运Mot r nCeal o方 法 以,二 分 项布 随机的数函 数
进 行 重抽样 ,复复 次 数 0 50次 S ,S9 1 重0 0用 A . 软3 件编程 拟 , 计算各 种 方 的 法 I型总 错误率 a (i fm模l y w s—r r aW, 、E ier V tR)检 验 能效判和对 率。 e o r e '1 多重
较比 法 方.
()多 总个率全部体不等, 2 考察 方法的检验各 效 能 : 算其 完 全 效 能 (—pl io e
) 例 效 比能 计a l a psw r 和 r
( e— aoe )p r ppi wr 大小 。的 全效完能 指在 多 重 较比 中
r 现发 了有 所差 异 能力 ,的比 例 效 能 差指异 被发 的现平 均 概率 ¨ 。 J
)(部 总分 率体相 等 ,合考 各 察 法方 F R的、3 综E W 检 效验能 及 判对以率 的高 低。 判对率指 多对重比较 务任 推的断完 全正确 的 率 概】
4 参。 数设置 .
本样率多重 比 的较方 法按算 法 的 不同 , 划可分
4为 类 :¨
) 1 (调界整值表 :法 杜养法(志u 和 罗海文法(u ) o Lo)。(
)2调整 检验 水 准 法:步 整 的调 法方 有B n-e o单f r n r法 ( ) Sd ko 法 ; 调 步 的方整法 有 S eo B in和 ai 逐 p—
dt w f o oi法 (t B ) t noB n enr S e pSe、p— wd ia o Sn k法d
检
水验 d:00 准;率 差80=0 ;数 k 取3 ~.5 .9 组 961 ; 、2、总 体取 率值 0 0. 9整数 倍 的:90.1 …09 00、 .8.
;样 量 本为n 5 、1 0 4 、01 两0 间样 本组量 例 比 、02 、 0 、 68 ;
为1 11 2 1 3 1 4 :、 :、: :。、 结 果
( S t、SeS )et pp—u oh eg法 ( oh S、e H cpr bHc )t pu— p
Hm l( o o m e 法Hm me )t—pu o法 ( o1 Se、pR mR )和 B nm an— crb( —。 jmeiiHo eh法gB )H
() 于检验 的方 : 法KN 法、c efM aaic3基 S S h fe —rsu— l c e (e、 c e e—Ch 法 ( on e) S f h o S 法h )f S hf e Coh n ce、fe
S —a 法 ( nea ) v Lem aS a n、 ey 和法R m yan 法
。1多
个总 体率全部相 等 .在模 拟的 数 组参合 条件下 , F R的 大 值最 按W
E(xF ) 这 R些方法 为 分类3: x W MaE W E将 MaF R 超 过
) 于重(抽样复 方法的 :tro 4 B基 o p和 Pt ru— s aet ma n it抽 样方 。法 文中 模 拟 o tB p 抽样 采 用 中 了o ot sr Femaa e n—rT kyD u l As re Ts 法 ,记 为 eu o eb c et算i n 简B o e; unt运用中 Fs了e'atee 算 ,法o trti Pm a
o ir c t ths x s
记简 为r P e 。m
2 ; Wx E M a c xF R于介 与 2 之a ;间 x E MW aF R 在以 内 。按再 F 与 的R对误差绝大的小排序, W 绝E对
差较 小误 的 种 几方 法 是:e < B o S<< R a P r o tmN K ny
<
L v<ey B— H;较 大的几 种法方为 L o> o De>uu Cnh
>Sf ( 1h。 c f 表e )e
2
价评 标指 .△通讯作
者 :平雁 陈,— i:h yn9 2 o . .maE ec p 9@6 1 cr nl
aM WEx R大F于 2 [4种 方 法 中 :u 的 0 L 法oF W-
RE值的高达 0 7 5, .7 大于 的比8是例 9% , 0 明其不
能控说制F , R 于率适 的多重比较 。 DEW不 u 的 法FR E
W
・
21 3・
C
iee uJ fHe l tt at s A p 0o. 6V , 2o h s onr o la taSi , rin h s c0 92, 21N .
最大 值为 0 18,W E. 0 F5 落 在R 与 2r间 的 情形 o之也 较 多( 0左 右 从 7 ) % 分位数 看 , ema 和 7% 。 5 S an hS f 部在情分 下形 不能 制 F控R。 ce ef W EMa
WE x F介 R 与 2于的 x种 方 法 7的对 绝 误差 c
F 非 常R近 接 ,WE 于 O大比例 很小 , 明 其检 分 布验 r . 说 非常接 近实际分 布 引。
M aW xF E在 R 内以方 的 法 ,中o e的绝对 Ch 法
n
误差较大 均, 和值7 %分 位数最 小 , 该 明 方法 的 且5 说果结很保守 。
值较都小 P,r Bo 和 BH— 这 3 方种 的 M法 a em且、 tox
表
1 个多总体全率部相时 , 等拟结模果 按a MWExF 排序R
各参
数条 对 F件 值 R的影响 是 :随着 组 数k W E①
增的大 ,EW F R的值一 般会递 减。C h n和 Shf eoc e法 e
平 水。③随样本量 n着增 大,的WE RF一般逐渐
增大 (值或减 小) 于 。 如 S aa 趋 n和 ehS mfc ef , n e法随
最 显为著 , 组随 数k增 大而的 剧急下降 。但 a S ae n法 m 在
数达 到组 后 , 出6现 FR值 迅速 上升 出超0 的 会 WE【
的 增大 F 值R逐 渐降 低 趋 于仅,而 控 制 F R 。 EW W从 E④ 不 平衡设 计 对 v L影的 显响 著 ,本量 悬殊样 e时法y在 其F R 下值 降剧 ( 烈
2 。 E W表)
象。②随着总现体接近 率 0 各5方的法F R都值 ., WE 逐 渐 大增, 以 Du法 S a和 最a为 显 , 明远超 出 e n法 m
远表
2多 个 体率全 总相 部时 等,总体 率为 4 05的平.衡设计 F下 的模R结拟 果 W
E
2多个体率总全部不等 .此参在条数下件, 比效例均能值小排序有大u :按 L > Dou> S KNB>H >yn > o tPn — R aB o > e >n Hml oe。m 完全效能较 的高法是方 u: >u >NB — Lo D SK H>> oh H / coHm eR m> t o e ( 3。由 于比 例效能 和 ml o B o>P 表 ) /n 全效能完度相高(关 =.9 ,00 0) 后面列仅出 059P<r. 1 ,0 比例效能 的果。结各参数条对检件效能验的响影:是 总
率体 的置位靠 0近会降5低验检效 能 ;量 n增会 大 样本
.
高检提效能;验不平 设衡会计降低检效验能在。上述上
升
和下降变化过的程中,各 方种法检验效大小的能对 相顺序基本不
变 。且并拟模 结的也 果证 验结了论 :逐步调
整检验水准 S p 、的 tS法方结的优果相于 应单的 t B pSe e
步法 方 nB Sd法 。k 、oi a 部3分总体 率 相等 .
中国 卫生 统 计20 0 9年 4月 第2 6卷第 期2
3 表多个总体率全不等时,部 拟结果按检模效验的降能序排 列
此条件 下,u L o法和 D法 u仍不 能 制控F R, WE 远 远 超出仪 平水 ; Rm法 不 能 既 控制 R且 检F 验效 o EW 很能低 ;3 方法种不 宜用 于本 率样 的重多比 较。
其他方法此判对率按均的降序排值有:列N K> BS
—
H>Ry n Pe >m>Bo >tHo n l>H oh> S ea r o mr ec tpS>
SBe , t 比p例效 能的结 果 排也类 序似, 只 S是 an的 em a法 比 效 能有例所 高提 另。 ,外时 Sa a 此 en Sh和 fm ef e可 c 控以 F e制( 4R。 r E表) 各参数 件条模对 拟结 果影 响显著 的 是 : 组数 ① k 总体:来 说比例 效能 随着组 数 的增 大 而呈 k下降 趋 势 , 但 B—H方 法比 平较 稳, = k 6时其当检 效验能 超 过已
S
K 。 N②样 本量n 大部 方法分着随 n的大增 :,检验 效能
提高, FxR也 增 。大 S和 B— H方 法 当n M aE NKW 很时 大, 能 制控F R , 不W 如 E1=时0 ,NK 和—BH6 S 方法
的 RF最 大 达值到 0 0 1EW 24和. 00 ;5. 99n =8 0 ,时两 的者 RF 最值大 00为 5WE .8 和2 005 。而 R- .67 y
一直a 严能格 控制 Z R( 5。 法nW 表 E)
4 部分表总体不率时 ,等 模拟结果按判率对降的序排
列 讨
论
影响很 大 果, 不 稳 健S。hf结 cf e的 3e种 方 ,法 提能 虽 供信置区 间 ,只有 Ch 法n能严 控 制格 FR, 其 但 o eEW且 结 果太 保守 ;mea aSn和 S hfc e 法在fn 较 ( I4 )大e 方n 0 > 时, 较 易容控 F制 R在。 数组为 3 ,e时a EW aS 法n m 的果 较结 敏感, H过 ml 法 方R。a 在各 情 种 超 om 等ey n法 下 形 验检效 能 较高 能 且较 好 控 制F R ( 高 W 最
00 7E) 52 .。
_l调 界值整 法 表: u和 uD , 然检验 效 L 能法 法 虽 较 o高, 而 大 在多 数情 形 其下 F R的 都值远 超 出 然W E
O, r不能 控 F R, 制实 际应用 。 . 以WE 难
2 .调整检
水准验 法o :法不 适用于 样 率本 多的 m
R重比较
; 逐步 方法(t SeBS优于 相应 的单步方法 S p 、 t ) e
( p 、ioa ) BSnd k逐 。步 方 法 :B 中H 方—法 检 效 验最能
4
于重 复抽 样的 方 :法 o Ptn. B基o 和 en方 法结
高的 ,在但 大较 的时 候,不能严 格控 制F R (WE最 高 05 0 )9 . 9。其 4 种方法他都严格控能制 FR, WE其 验检 效 能显 H示m l oh>SeHS>S eB o。 me > ctp t p3 于
验检 的 方法: K N 法的 果 结检验 效 能最 . 基S
较果接为近, 检 验效能高于 H m , 其 eom l能较好 控制
 ̄VR( 高 004 。 E 最 5. )
上综述,所在 样本率多重的 较比 中, 判对率高较
方的法 :有SN 、K —、 yHn nP 、ot Hm em、 B R 、a e oB、o l
,高但 在, z 80容易时现 F出 值较R高 情的形( > IWE 最高 0 1 0 )r法e 检验 效 较 低 且 受能样 本 量比 例 . 。L 2、,4
H)c 、 tSS p 检验且能效也致大依顺此递减。 序o S h pt、 , B e
其 中 H me loH heSSS e B 严能格 制 控 RFo me、 c 、tp 、t pW E
・134 ・
hC nso amf H el t tt s . r 20 Vo .6,. eeJi
u ol a taS iiA p 0。912 o 2 Nh c
s在
仅 以 内 ; y n和 P r B Rao e、mo t能 较好 控 制 RFW
E量 较本情 小下形, 才 能 基本 控制F R( .o 以 ) 内E 1 3W r 如模, 中拟S NK 需n O 且 , H尼则 需n 8 。≤ 4≤5 B—<  ̄0
(5.1 a以 内 1 S );KN 、—H虽 效 能较 高 ,只 能在 B样但
表
部5分 体总率相时 , 等平设衡下计比例效能和 Ma EWxF 的R模拟果结
te t
cn lg fs mauo a n npS 1 rag .eR lsu heh o yo o iltn s f gm ASi9 3. po r m s i si t
结
论
E
ah me oa erc a a ts c r o . cr d r saes n i v n a c t d h st ih ce t sS mep o e u e l e ts ae d cn h h r i i i
综
上所 述 ,据数处 理 中 各 种 法方的 应 用可根据 对F R的 制控 要求进 行选 。 择E W1 如.严需 格 控 F制R 在 水 准 内 ,E 按敏感 性从W 大 小选 到有 :择 oH mm、 e ho、 t S eS等 B逐步调 l H ceS、 t p 整p检 验 水 准 。其法中 o Hm m l的e 结果 最 敏感 而, Se SS eB方法 大 则较熟家悉, 计易算 。 t、pt p 容
cn tFl R Etnwe hls c a eP ,m ,om e t.mn- o r W o o h oe u ,haR y, r Bn too s Hec lCo c ui n I WER rqsie sr lc o t ,Hoon loH ,htp l os F i f uredt t c rn ily rm e c Se ,s
,tS pcn opi e a eut or f FE Rs al w el hy o e foe B a r v ddq ae p we .I W i l ods i t rv wl g
l
(tn 1 1 R c ,ean dBaotr es v t o wi .5i y t nPn no senie me  ̄ .h ), i a th
【
e o d】 Mi ecu msna miF—w es ore Kyw rs l o lrop ; al t pi y ir ra e r ;
tP o ewr; r T e ouem ld
2 如. 允 F许 R略 出超o 水 准( 拟 中1 1 EW l .5模 内 )
按敏 性从 感到小大 选 有 :y择nP n , Ra 、ne 和 oB 方to 。P法 rB o 方 需 法在S 环S境 下运 实行现 ,y enm、ot A aR 法 可手计 工 。 算
26 l_ O 22.
参
考
文
献
1钱 ,俊. 陈平 雁.多 个本样 的多重率比 较. 中国 卫生计统 0, 28( ) 20 5 2,
:2颜
主虹编 . 医学. 统计 学 .京 : 民 生 卫版出社 , 50 810 1北 人 02, 5 —6 . 3徐 勇 勇主 编. 统 计学学 . :京 等教 育出 版社 ,0 ,18 2. . 医北 高 2 08 - 3
3 能提供.置信间区 的 种3方法结果不敏都感, 如
需 提确供 信 区 间置, ,4 在 l≥ 的0情 形 ,下当组 数 3 为 时,择 S a 法n, 之就 运 用Shf方 法 。 当 < 选ne am反 ce f e 4 0时 , 采用 Cn法h ,可o e但其 结果非 常保 。
守 Sm ua niS ud f i-Ps mp sr n n P o ot n i lt ot yo a rw e Ci oai o is r p r sio
4 钱 , 俊. 陈雁.平 验组实对与组比较照的样本 率重 比多较方法的 模拟研
究 . 国 卫 生统计 ,0 ,8 55 :0- 0 中 20. 2( 5 0) 5 3
.5e afl PM ut l oa sn n utp e ts u i e S - W sy l H. l pe mc r po a sdm l tls ii i es :g n t SA s s
het . mrN : ACS I s it I , c99 a y, SC n t u ,e 1n9 .t
陈平雁6 , 王斌会 . .四表资格料 三种常用统计方法的 比 .较国 卫生 统 中
计 9 6, 1( ) 1— 69 1, 32:4 1.
. 7ema Sa nMA w. e oeme t av s r plec m p sr n: o N d e lvp i np i ,em u t l o ais S me sn i o pio ru p a dtc ebr c du e . y h lg c l B ltnu 9,1, 1 0wel a c ar lp o er Pss coo ia l i 1 f 9n ie 1
Q
au , h nP nyn Deat o is nts St u miJ Cne iga p r.t f Bo ii , t ̄eo n em asc
M ie U l rii (1 5 ,5 u no z d cnv ys5 0 1 ) aGgh u a t
e( )
57 65
:37- 8.
. 8bi iCr b e R .Ai iw l pe aion :A d e c m ps n pr ar P smut l mc rps se oi om l ao i a- o
【
b t 】c0 jc v T r ̄ t arcc op rw eAs a tr etb oe xo e h ae s sf a — i i e eh c trt i s i
i cm r po r sdc r rp ro s uig M otC rto s Me h oi na o e e ui po t n s nea l med . t -s psn o i n oh
soF R , o a dr Ime e o o c r cd r cr pmrd u g di WE P ew nu dl f ah po e eae o aue s e r n
p o c es tp s r dcrs Brt so ra f Mama c l ad ra h rvssu ei wp eou e .i i J hun l oe t het a i
nStt c sPls co o y ,03 5:,7 3. a at y lh g 0 252 9 - ii