8. GSD-條件把握度的計(jì)算

隨機(jī)縮減(Stochastic Curtailment):

?????? 指在當(dāng)前累積數(shù)據(jù)的基礎(chǔ)上雀哨,預(yù)測(cè)試驗(yàn)的最終結(jié)果有較大的可能性有效或無效市俊,因此可以提前中止試驗(yàn),對(duì)應(yīng)的過程稱為Stochastically Curtailed Procedure(SCP)。在試驗(yàn)進(jìn)行過程中的任意一個(gè)階段,通過計(jì)算SCP的相關(guān)指標(biāo)與某一提前設(shè)定的固定值γ來判斷試驗(yàn)是否可以提前終止谒府。
?????? 相關(guān)指標(biāo)包括條件把握度:Conditional power tests(frequentist approach);預(yù)測(cè)把握度:Predictive power tests(mixed Bayesian-frequentist approach)浮毯;預(yù)測(cè)概率值 Predictive probability tests(Bayesian approach)完疫;非參數(shù)法:Parameter-Free Approach。等等债蓝。常用的是條件把握度壳鹤,在適應(yīng)性設(shè)計(jì)的樣本量重估中也會(huì)用到。
?????? SAS在Proc SEQTEST里可以進(jìn)行計(jì)算饰迹,PASS芳誓、R(gsdesign package)中也可以計(jì)算讯嫂。

  • Conditional Power
    CP指基于當(dāng)前累積的數(shù)據(jù),試驗(yàn)結(jié)束時(shí)拒絕原假設(shè)的概率兆沙。SAS中的命令為Condpower <(options)>
    • options包括:
      Type=Allstages | Finalstage
      ?????? 默認(rèn)為Allstages欧芽。CP的原始定義以及實(shí)際臨床試驗(yàn)中使用的的Conditional power為FinalStage。
      CREF=c1... ...葛圃,默認(rèn)CREF=0 0.5 1.0 1.5千扔;
      ?????? 指CP計(jì)算所基于的效應(yīng)量相對(duì)于備擇假設(shè)的比例。假設(shè)θ為CP計(jì)算所基于的參數(shù)值库正。θ_k為第k次期中分析時(shí)的效應(yīng)值曲楚,θ_1為備擇假設(shè)值,CREF=θ/θ_1褥符。CRef=1時(shí)CP計(jì)算基于備擇假設(shè)龙誊。如θ_k=-0.0444,θ=0.09喷楣,則此時(shí)CRef=-0.0444/0.09=-0.4938趟大。CRef相當(dāng)于設(shè)定最終分析時(shí)的效應(yīng)量與備擇假設(shè)之間的比例,來計(jì)算不同情況下的CP铣焊。
    ???? - Allstages: 計(jì)算未來所有階段分析拒絕H0的總概率逊朽。

???? ???? 相當(dāng)于將各階段拒絕H0的概率相加,z_k為第k次分析統(tǒng)計(jì)量曲伊,θ為備擇假設(shè)叽讳,a_k,b_k為有效或無效終止時(shí)的界值。

SAS Help

???? - FinalStage:計(jì)算最后一次分析拒絕H0的概率坟募。

???????基于第k次期中分析的統(tǒng)計(jì)量Zk計(jì)算最后一次分析時(shí)拒接原假設(shè)的條件概率岛蚤。PASS和SAS中給出的雙側(cè)檢驗(yàn)CP公式,兩個(gè)公式換算之后是一樣的懈糯。

??????? θ:假設(shè)檢驗(yàn)的參數(shù)涤妒, 原假設(shè)下θ=0,備擇假設(shè)下θ=θ1昂利; I_k:第k次分析時(shí)的信息量
??????? I_K:最終分析時(shí)的信息量届腐; Z_k:第k次分析時(shí)的統(tǒng)計(jì)量; Z_{a/2}:a/2對(duì)應(yīng)的分位數(shù)

???????PASS中給出的一個(gè)根據(jù)公式直接計(jì)算單樣本率單側(cè)檢驗(yàn)的例子蜂奸,α=0.025犁苏,Z_{1-α}=1.96p_1=0.65扩所,p_0=0.55围详。備擇假設(shè)θ=p_1-p_0=0.1σ^2=\hat{p}(1-\hat{p})\hat{p}=(p_1+p_0)/2助赞,z_k=2买羞。


???????SAS中給出的例子:
??????? Example 104.2,共4次分析雹食。H0:θ=0 vs. H1:θ=0.10畜普, θ=p-p0。Stage=2時(shí)群叶,Pdiff=-0.064789吃挑。
??????? 這個(gè)例子用上述公式直接算和SAS輸出的結(jié)果有一些差異,CP差了0.03左右街立,不確定是什么原因?qū)е碌牟俺摹AS中率的I_k=n_k/σ^2,σ的\hat{p}似乎是基于備擇假設(shè)p_1赎离,而不是(p_1+p_0)/2逛犹。

proc seqtest Boundary=Test_Prop1
         Data(Testvar=PDiff Infovar=NObs)=Data_Prop2
         infoadj=prop
         boundarykey=both
         boundaryscale=mle
         condpower(cref=1 type=finalstage)
         predpower
         plots=condpower
         ;
   ods output test=Test_Prop2;
run;
* Nstages=,設(shè)定一共有幾次分析梁剔,Nstages=3則輸出基于共3次分析計(jì)算的各階段樣本量虽画、界值等。
  • Predictive Power
    PredPower <(options)>
    在目前已觀測(cè)到的統(tǒng)計(jì)量及給定一個(gè)先驗(yàn)分布的情況下憾朴,最后一次分析時(shí)拒絕H0的后驗(yàn)概率狸捕。P_k(θ)為條件把握度,π為給定當(dāng)前期中分析數(shù)據(jù)后D(K)及θ的先驗(yàn)分布后众雷,θ的后驗(yàn)分布,作為權(quán)重做祝,計(jì)算加權(quán)平均的概率值砾省。

參考: - SAS help、
???????????? - PASS help混槐、
???????????? - Group Sequential and Confirmatory Adaptive Designs in Clinical Trials

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末编兄,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子声登,更是在濱河造成了極大的恐慌狠鸳,老刑警劉巖,帶你破解...
    沈念sama閱讀 218,682評(píng)論 6 507
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件悯嗓,死亡現(xiàn)場(chǎng)離奇詭異件舵,居然都是意外死亡,警方通過查閱死者的電腦和手機(jī)脯厨,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,277評(píng)論 3 395
  • 文/潘曉璐 我一進(jìn)店門铅祸,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人,你說我怎么就攤上這事临梗∥卸螅” “怎么了?”我有些...
    開封第一講書人閱讀 165,083評(píng)論 0 355
  • 文/不壞的土叔 我叫張陵盟庞,是天一觀的道長(zhǎng)吃沪。 經(jīng)常有香客問我,道長(zhǎng)什猖,這世上最難降的妖魔是什么票彪? 我笑而不...
    開封第一講書人閱讀 58,763評(píng)論 1 295
  • 正文 為了忘掉前任,我火速辦了婚禮卸伞,結(jié)果婚禮上抹镊,老公的妹妹穿的比我還像新娘。我一直安慰自己荤傲,他們只是感情好垮耳,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,785評(píng)論 6 392
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著遂黍,像睡著了一般终佛。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上雾家,一...
    開封第一講書人閱讀 51,624評(píng)論 1 305
  • 那天铃彰,我揣著相機(jī)與錄音,去河邊找鬼芯咧。 笑死牙捉,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的敬飒。 我是一名探鬼主播邪铲,決...
    沈念sama閱讀 40,358評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼无拗!你這毒婦竟也來了带到?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 39,261評(píng)論 0 276
  • 序言:老撾萬榮一對(duì)情侶失蹤英染,失蹤者是張志新(化名)和其女友劉穎揽惹,沒想到半個(gè)月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體四康,經(jīng)...
    沈念sama閱讀 45,722評(píng)論 1 315
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡搪搏,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,900評(píng)論 3 336
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了箭养。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片慕嚷。...
    茶點(diǎn)故事閱讀 40,030評(píng)論 1 350
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出喝检,到底是詐尸還是另有隱情嗅辣,我是刑警寧澤,帶...
    沈念sama閱讀 35,737評(píng)論 5 346
  • 正文 年R本政府宣布挠说,位于F島的核電站澡谭,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏损俭。R本人自食惡果不足惜蛙奖,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,360評(píng)論 3 330
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望杆兵。 院中可真熱鬧雁仲,春花似錦、人聲如沸琐脏。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,941評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽日裙。三九已至吹艇,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間昂拂,已是汗流浹背受神。 一陣腳步聲響...
    開封第一講書人閱讀 33,057評(píng)論 1 270
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留格侯,地道東北人鼻听。 一個(gè)月前我還...
    沈念sama閱讀 48,237評(píng)論 3 371
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像联四,于是被迫代替她去往敵國和親精算。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,976評(píng)論 2 355

推薦閱讀更多精彩內(nèi)容