2. 各階段樣本量相等的情況(Pocock&OBF)

?????? 各階段等樣本量的情況下瓷耙,前面提到的各階段統(tǒng)計量及協(xié)方差

?????? Z_{k^\sim}為各個階段的統(tǒng)計量圆到,n_{k^\sim}為對應(yīng)的樣本量鲤桥, Z_k^*為累計至階段k時的統(tǒng)計量。
?????? 可以簡化為
?????? ?????? ?????? ?????? ?????? ?????? ?????? Z_k^*=\frac{1}{\sqrt{k}}\sum_{k^\sim=1}^kZ_{k^\sim}拢操,Cov(Z_k^*,Z_{k'}^*)=\sqrt{\frac{k}{k'}} for k≤k'

?????? 常見的方法包括Pocock锦亦、OBF以及Wang and Tsiatis。Wang and Tsiatis通過修改特定的參數(shù)可以產(chǎn)生不同的界值令境,當(dāng)界值為0.5時即為Pocock杠园,當(dāng)界值為0時即為OBF。

  • Pocock

Pocock方法使各個不同階段分析的界值相同舔庶,即針對Z_k^* (k=1,...,K)抛蚁,u_1=...=u_K=u'
使P_{H0}(|Z_1^*|≥u',...,|Z_K^*|≥u')=α
u'與給定的α以及總的分析次數(shù)K有關(guān),因此標(biāo)記為c_p(K,α)惕橙。
即Pocock的接受域為(-u';u'),u'=c_p(K,α), k=1,...,K

  • O'Brien and Fleming

Pocock是各次期中分析的界值均相同瞧甩,而OBF方法隨著期中分析越往后界值越小。
P_{H0}(|Z_1^*|≥u_1,...,|Z_K^*|≥u_K)=α
u'與給定的α以及第幾次分析k有關(guān)弥鹦,因此標(biāo)記為c_{OBF}(K,α)/\sqrt{k}肚逸。
即OBF的接受域為(-u_k;u_k),u'=c_p(K,α)/\sqrt{k}, k=1,...,K

  • Wang and Tsiatis

是在Pocock和OBF方法上的擴展,通過更改位移參數(shù)delta可以獲得不同的界值彬坏。
接受域為(-u_k;u_k),u_k=c_{WT}(K,α,△)k^{△-0.5}, k=1,...,K
當(dāng)△=0.5時為Pocock朦促,此時不同期中分析的界值相同。當(dāng)△=0時栓始,此時為OBF务冕,越接近最終分析界值越小,在樣本量較少的早期期中分析時越不容易拒絕原假設(shè)幻赚。
通過u_k=c_{WT}(K,α,△)k^{△-0.5}計算各個階段的界值以及一類錯誤禀忆。
梳理給出了不同K、α坯屿、△時c_{WT}(K,α,△)的取值以及膨脹因子油湖。

Group Sequential and Confirmatory Adaptive Designs in Clinical Trialsg
三種方法的對比舉例

  • 雙側(cè)檢驗對稱界值

前面的Pocock和OBF僅考慮了有效終止,界值為(-u';u')(-u_k';u_k')繼續(xù)試驗领跛,反之有效終止乏德。這里針對雙側(cè)檢驗進行有效或無效終止,采用對稱的界值(Pampallona and Tsiatis)吠昭,即

  • Z_k^*∈(u_k^0;u_k^1) or Z_k^*∈(u_k^1;u_k^0)時繼續(xù)試驗
  • Z_k^*∈(-u_k^0;u_k^0)時無效中止
  • Z_k^*∈(+∞;u_k^0) or Z_k^*∈(-u_k^0;-∞) 時有效中止
    在最終分析時喊括,u_K^0=u_K^1

---計算方法---
基于前面的Wang and Tsiatis方法,需要考慮無效中止矢棚,因此額外加入β

計算公式:
?????? u_k^0=v_k-c^0(K,α,β,△)k^{△-0.5}, k=1,...,K郑什,
?????? u_k^1=c^1(K,α,β,△)k^{△-0.5}, k=1,...,K
?????? v_k=\sqrt{k} (c^0+c^1)K^{△-1}
c^0,c^1只跟α蒲肋,β蘑拯,△钝满,K有關(guān),給定了想這些參數(shù)后申窘,可以計算出c^0c^1以及v_k弯蚜,進而計算出各個階段的界值。

書中table2.7及2.8給出了不同△剃法、α碎捺、β、K對應(yīng)的c^0贷洲、c^1以及樣本量膨脹因子收厨。

k* denotes the first stage where H0 can be accepted. In parentheses: expected reduction in sample size under H0, the value midway between H0 and H1, and H1, respectively

計算舉例
取△=0,α=0.05优构,1-β=0.8诵叁,K=4。根據(jù)表格2.7俩块,c^0=1.9892, c^1=3.9055黎休。
v_k=\sqrt{k}(c^0+c^1)K^{△-0.5}=\sqrt{k}(1.9892+3.9055)4^{-1}=\sqrt{k}1.4738
????u_1^0=1.4738-1.9892<0
????u_2^0=\sqrt{2}*1.4738-1.9892/\sqrt{2}=0.678玉凯,
????u_3^0=\sqrt{3}*1.4738-1.9892/\sqrt{3}=1.404势腮,

????u_4^0=\sqrt{4}*1.4738-1.9892/\sqrt{4}=1.953
????u_1^1=3.9055漫仆,
????u_2^1=3.9055/\sqrt{2}=2.762捎拯,
????u_3^1=3.9055/\sqrt{3}=2.255
????u_4^1=3.9055/\sqrt{4}=1.953盲厌,
????因此(u_1^0,u_2^0,u_3^0,u_4^0)=(-,0.678,1.404,1.953)署照,(u_1^1,u_2^1,u_3^1,u_4^1)=(-,2.762,2.255,1.953)

  • 單側(cè)檢驗

單側(cè)檢驗:H_0:μ≤μ_0H_1:μ>μ_0
對應(yīng)的
????* Z_k^*∈(-∞;u_k) 時無效中止吗浩,即接受H0
????* Z_k^*∈(u_k;+∞) 時有效中止建芙,即拒絕H0

控制一類錯誤需要滿足:P_{H_0}(Z_1^*≥u_1 or ,...,Z_K^*≥u_K)=α
可以采用前面提到的Wang and Tisatis來控制一類錯誤,產(chǎn)生有效終止的界值懂扼。
無效終止書中提到的方法DeMets and Ware禁荸,此時無效終止的界值為一個常量u^L
????* Z_k^*∈(-∞;u^L) 時無效中止阀湿,即接受H0
????* Z_k^*∈(u^L;u_k) 時繼續(xù)試驗赶熟,
????* Z_k^*∈(u_k;+∞) 時有效中止,即拒絕H0
當(dāng)u^L=-∞時即僅考慮有效中止陷嘴。


在固定樣本量試驗中映砖,單側(cè)和雙側(cè)檢驗通常可以采用同一個界值灾挨。但是在成組序貫單側(cè)檢驗中邑退,兩者不一定相等竹宋。
參考:Group Sequential and Confirmatory Adaptive Designs in Clinical Trials(by Gernot Wassmer &Werner Brannath)

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市地技,隨后出現(xiàn)的幾起案子逝撬,更是在濱河造成了極大的恐慌,老刑警劉巖乓土,帶你破解...
    沈念sama閱讀 206,839評論 6 482
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異溯警,居然都是意外死亡趣苏,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,543評論 2 382
  • 文/潘曉璐 我一進店門梯轻,熙熙樓的掌柜王于貴愁眉苦臉地迎上來食磕,“玉大人,你說我怎么就攤上這事喳挑”蚵祝” “怎么了?”我有些...
    開封第一講書人閱讀 153,116評論 0 344
  • 文/不壞的土叔 我叫張陵伊诵,是天一觀的道長单绑。 經(jīng)常有香客問我,道長曹宴,這世上最難降的妖魔是什么搂橙? 我笑而不...
    開封第一講書人閱讀 55,371評論 1 279
  • 正文 為了忘掉前任,我火速辦了婚禮笛坦,結(jié)果婚禮上区转,老公的妹妹穿的比我還像新娘。我一直安慰自己版扩,他們只是感情好废离,可當(dāng)我...
    茶點故事閱讀 64,384評論 5 374
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著礁芦,像睡著了一般蜻韭。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上宴偿,一...
    開封第一講書人閱讀 49,111評論 1 285
  • 那天湘捎,我揣著相機與錄音,去河邊找鬼窄刘。 笑死蕉汪,一個胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的枕扫。 我是一名探鬼主播,決...
    沈念sama閱讀 38,416評論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼烹骨,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了材泄?” 一聲冷哼從身側(cè)響起沮焕,我...
    開封第一講書人閱讀 37,053評論 0 259
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎拉宗,沒想到半個月后峦树,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,558評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡旦事,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 36,007評論 2 325
  • 正文 我和宋清朗相戀三年魁巩,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片姐浮。...
    茶點故事閱讀 38,117評論 1 334
  • 序言:一個原本活蹦亂跳的男人離奇死亡谷遂,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出卖鲤,到底是詐尸還是另有隱情肾扰,我是刑警寧澤,帶...
    沈念sama閱讀 33,756評論 4 324
  • 正文 年R本政府宣布蛋逾,位于F島的核電站集晚,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏换怖。R本人自食惡果不足惜甩恼,卻給世界環(huán)境...
    茶點故事閱讀 39,324評論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望沉颂。 院中可真熱鬧条摸,春花似錦、人聲如沸铸屉。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,315評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽彻坛。三九已至顷啼,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間昌屉,已是汗流浹背钙蒙。 一陣腳步聲響...
    開封第一講書人閱讀 31,539評論 1 262
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留间驮,地道東北人躬厌。 一個月前我還...
    沈念sama閱讀 45,578評論 2 355
  • 正文 我出身青樓,卻偏偏與公主長得像竞帽,于是被迫代替她去往敵國和親扛施。 傳聞我的和親對象是個殘疾皇子鸿捧,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 42,877評論 2 345