1. 相关概念

       临床试验中缺失数据无法避免,常用的填补方法包括单一填补(Single imputation)和多重填补(Multiple imputation)。单一填补方法包括LOCF(Last observation carry forward)、BOCF(Baseline observation carry forward)等。LOCF、BOCF这类方法没有考虑数据的变异性,会降低数据的标准误,增加假阳性的可能性。多重填补针对同一个缺失点产生多条填补记录,之后再合并进行分析。MI可以作为一种填补方法,也可以是实现其他一些填补策略的手段,比如在J2R、CR、CIR等填补方法里,用MI来产生多条填补数据后再进行处理以符合既定的填补要求。

数据的缺失模式:

       1. MCAR(Missing at Complete Random):缺失数据与已观测到的数据及未观测到的数据均无关。
              比较理想的缺失情况,不太可能出现,因此通常不考虑MCAR。
              如受试者因车祸等偶发事件失访。
       2. MAR(Missing at Random):缺失数据与未观测到的数据无关。
              如受试者因疗效不佳退出。
       3. MNAR(Missing not at Random):缺失数据与已观测到的数据及未观测到的数据均有关。
              该缺失模式难以判断,因为与未收集到的数据有关,所以无法判断到底是否MNAR。

       缺失数据未观测到,所以无法假定它们实际应该是什么样的以及到底服从哪种缺失模式。在实际分析中,主要分析通常基于MAR,并且需要增加基于MNAR的敏感性分析。

缺失类型:

       1. 单调缺失(Monotone):在某一访视点及之后的访视中数据均缺失。
       2. 非单调缺失(Non-Monotone):在某一或几个访视点数据缺失,但后续又收集到了其他访视点数据。
       *:临床试验数据中出现非单调缺失可能是因为一些非试验相关的 因素导致,如受试者外出、封控等,导致某几个访视点的数据未收集到,与试验无关,同时该类缺失的比例一般也较小(即使有non-monotone,影响也较小)。因此在后续的缺失数据处理方法中,如果出现非单调缺失数据,通常假定其为MAR,先将其处理为单调缺失,之后再根据相应填补策略(如MNAR下的填补方法)进行处理。
       填补为单调缺失可以通过Proc MI的MCMC命令实现。
       MCMC impute=Monotone和impute=Full 分别对应填补为单调缺失及全填补。

  Proc MI data=XX Nimpute=XX Seed=XXX out=XXX;
  Var  XXX;
  MCMC impute=Monotone;  
  /* MCMC中的其他命令: Chain=Signgle/Mutiple,蒙特卡洛单链或双链;
  /* NITER=100(默认)或 NBITER=200(默认); 设定迭代次数,保证记录间的独立性 */
  run;
敏感性分析:

       主要分析通常基于MAR,敏感性分析进行基于不同缺失假设下进行。

缺失比例:

       不同组别的缺失比例可能造成偏倚,如试验组缺失比例过高,采用Reference-based的方法与阳性对照组相比低估疗效。

贝叶斯后验分布(Bayesian Posterior Distribution):

       多重填补时的参数从后验分布中获得,后验分布依赖于先验分布(Prior distribution)和似然函数(Likelihood function)。常用的先验分布是Jeffrey' prior(无信息先验),也是SAS的MCMC中默认的先验分布,即先验分布中未包含任何有用的信息,后验分布主要依赖于观测到的数据。-Joint Modeling

MCMC(Markov Chain Monte Carlo):

       后验分布无法直接算出时,采用MCMC的方法进行模拟以获得分布相关的信息,并从中获得填补值。MCMC包括两步:
              1. Impute-step:通过已观测到的数据的信息(均值、协方差矩阵)产生一系列填补值。
              ( the I-step draws values for Y_{i(mis)} from a conditional distribution for Y_{i(mis)} given Y_{i(obs)});
                 SAS中的NBITER=X(“Burning-in”)可以定义该步的迭代次数,即迭代X之后的参数采用于填补。
              2. Posterior-step:1步缺失数据填补完整后,根据填补后的数据产生新的均值协方差矩阵迭代步骤1。
       重复步骤1和2直到模型收敛于一个稳定的分布(stationary distribution)。步骤1中的填补值并不用于真正的填补,只有迭代收敛后的值才用于最终填补。模型并收敛到某一个特定的点估计,而是在特定分布周围随机波动,每次填补时都产生不同的参数,并生成填补值。

MNAR下填补的三大框架:
  1. Selection Models
    p(Y_{obs}, Y_{mis},R|X) = p(Y_{obs}, Y_{mis}|X) p(R|X, Y_{obs}, Y_{mis})

  2. Pattern-mixture models
    p(Y_{obs}, Y_{mis},R|X) = p(R|X) p(Y_{obs}, Y_{mis}|R,X)
    PMMs中受试者被分成了不同的序列,同一个序列的受试者具有相同的缺失模式,如受试者的不同特征:缺失的原因、缺失所在的访视。
    p(R|X)为不同缺失模式下的概率分布,p(Y_{obs}, Y_{mis}|R,X)不同缺失模式的数据分析模型。

  3. Shared Parameter models

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 200,961评论 5 473
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 84,444评论 2 377
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 148,009评论 0 333
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,082评论 1 272
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,101评论 5 363
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,271评论 1 278
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,738评论 3 393
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,395评论 0 255
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,539评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,434评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,481评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,160评论 3 317
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,749评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,816评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,038评论 1 256
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,548评论 2 346
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,140评论 2 341

推荐阅读更多精彩内容