3.1 MNAR下的Delta Adj&Tipping P方法(SRI)

  • Delta-Adjustment

      假设受试者在退出试验后,疗效将会降低δ(高优/低优指标分别+/-某一特定的δ值),按照该假设对缺失数据处理后进行分析。在单臂试验中尤为有用,因为此时无法采用Control-based填补方法。同样在等效或者非劣效试验中,对照组一般为阳性标准对照,可以认为它的疗效是大于或等于试验组,特别当两组脱落率接近时,采用control-based方法实际上并不保守。
      Delta值可以基于疾病特征、治疗效应等制定。可以在两组均进行校正、也可以仅在试验组进行校正。
      校正的方式比较灵活,如 a.仅在某一访视(如第一次缺失所在访视)进行校正后续访视仍采用Regression model;或 b.将所有缺失数据都填补完成后,在所有访视均进行校正;c. 在b的基础上不同的访视采用不同的δ。这里采用简单的Sequential Regression Modeling(SRI)实现。

例程序1(上述方法a):

  1. 通过MCMC将Non-monotone缺失处理为Monotone缺失
                MCMC impute=Monotone; 

  2. 假设缺失所在访视为2。
           PROC MI DATA=Unimputed_2 OUT=Imputed_2 NIMPUTE=1 SEED=XXXX; 
                    BY _imputation_; 
                    CLASS Group CatBase; 
                    VAR Group CatBase Conbase AVAL1;
                    MONOTONE REG(AVAL2=Group CatBase Conbase AVAL1); 
            RUN; 

   3. 假设仅针对试验组第2次访视进行delta adjustment,
       则当在第2次访视提前退出时(LastVis=2),在AVAL2上减去delta。
           DATA Unimputed_5; 
             SET Imputed_4; 
             If Group=1 and LastVis=4 then AVAL4=AVAL4-delta;
           run;
   4.  后续的其他访视重复步骤2(访视2的delta校正对后续访视的影响通过回归方程体现)
       
       补充:结局为分类变量时,方法类似。此时采用monotone logistic(AVAL2 = XX XX);进行填补,
            delta为将填补完成的访视结果校正为无应答的概率。
            在第3步中通过一个随机函数对AVAL2进行重新填补, 
            If Group=1 and LastVis=4 then do;
                  If ranuni(seed)<delta then AVAL2 = 0; *0=无应答;
            end;
  • Tipping-Point

     Tipping-Point方法是通过设置一系列的δ值以观察缺失值对最终结论的影响,并找到一个使结论反转(结论成立至不成立)的δ,该值称为Tipping point。通过对tipping point的医学评估,可以判断缺失值对结论的影响,以及不同缺失假设下结论的变化。如Tipping point很大,则说明结论较为稳健。
     Tipping-Point是delta adjustment的一种特殊应用,或者Delta-Adjustment是实现tipping-point的一种方法。可以将delta-adjustment包装在宏里,并设置不同的δ来实现tipping-point方法。

例程序2(仅对某一特定访视进行处理):

  1. 通过MCMC将Non-monotone处理为Monotone

  2. 对Monotone缺失数据进行回归填补,假设仅对第4次访视的结果进行校正。
     (可以仅对某一组进行校正。也可以对不同的组别采用不同的δ,如试验组减一系列的值同时对照组增加一系列的值。通过设置不同的δ找到临界值)
      Proc MI data=Imp_1 out=Imp_2 nimpute=1 seed=XX;
          Var Group CatBase ConBase AVAL1 AVAL2 AVAL3 AVAL4;
          Class Group CatBase;
          MNAR adjust(AVAL4/shift=delta1 adjustobs=(TRTPN=1))
               adjust(AVAL4/shift=delta2 adjustobs=(TRTPN=2));
          Monotone regression;
          By _Imputation_;
      Run;
     校正后的结果为y*=cy+δ,这里c对应SAS中的scale参数,δ对应shift参数。

     补充:a. MNAR adjust(V170/shift=0 adjustobs=(TRTPN=1))就相当于没有进行校正。等同于直接采用Regression进行填补。
              Proc MI 变量同上;
                  Var 同上;
                  Class 同上; 
                  Monotone regression;
              Run;
          b. 采用Proc MI的MANR(2)进行校正,等同于在直接进行Regression填补(a)完成后在填补后AVAL4上减去δ。结果是一样的
           Shift=0或Monotone Regression(a)    |    Shift=-5(b)    |    (a)-5
                    2.9                                2.9          该记录不缺失,不进行处理
                   6.848                              1.848             1.848                  

例程序3:
以MAR填补完成后的组间差值作为初始的delta,并按一定比例递增或递减直至结论反转。

1. 通过FCS(或者回归填补,MCMC)对缺失值进行填补。
    proc MI data=TransV2H seed=XX nimpute=X out=ImpuOut;
      by TRTPN; /*Imputed by Group*/
      class TRTPN country;
      Var TRTPN baseline country AVAL1-AVAL2;  
      FCS REG(AVAL1= country baseline);
      FCS REG(AVAL2= country baseline AVAL1);
      ... ...;
   run;

2. 根据项目要求采用ANOVA或者MMRM进行分析,计算差值

3. 以第2步计算出的差值为基础进行校正,缺失的访视数据按一定比例增加(或减少)K%进行校正,直到结论反转。
    proc sort data=TransH2V; by usubjid; run;
    data TransH2V;
      merge TransH2V /* TransH2V为填补完成后横结构转置为竖结构后的数据集 */
            LastVIS; /* LastVIS用于标记缺失开始的访视 */
      by usubjid;
      if VISIT>LastVIS and TRTPN='R' then  /*可以根据实际要求筛选,如组别、退出原因等*/;
      Adjusted=AVAL-Delta*K%;   *Delta: 第2步中计算的两组差值;
                                *K按10%递增直到P>0.05;
      else AdjustCh=change;
   run;

Ratitch and OKelly Pattern Mixture PharmaSUG-2011-SP04

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

  • 微生物群落研究的主要目标是比较不同群落的组成(β多样性)。在第6章介绍了β多样性,并举例说明了如何计算β多样性指数...
    ZMQ要加油呀阅读 4,786评论 0 5
  • Learning Human-Object Interactions by Graph Parsing Neura...
    Valar_Morghulis阅读 2,010评论 0 0
  • 个人为了看的方便没有广告,转载自泡泡机器人:https://www.sohu.com/a/161346283_71...
    Maxsium阅读 39,562评论 0 8
  • 1.你知道我们公司是做什么的吗?为什么想来我们公司? 网宿科技是做CDN,content delivery net...
    光_武阅读 3,822评论 0 4
  • 1 前言 本章我们将会探讨三维视觉,首先会研究三维到二维的投影变换以及其逆变换,然后会介绍多相机立体景深感知。在这...
    RichardJieChen阅读 5,125评论 0 3

友情链接更多精彩内容