非线性量效模型逐点置信区间估计及其相关问题研究

来源 :空军军医大学 | 被引量 : 0次 | 上传用户:a504468075
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
对于一个新药的研发,量效关系的研究是一个贯穿整个药物研发阶段的复杂过程。其中着重于剂量探索(Dose-finding)的II期临床研究是药物最适剂量选择的关键性步骤。在采用模型法的剂量探索研究中,量效反应模型(Dose-response model)的参数估计以及用于刻画模型拟合不确定性的整个曲线逐点药物效应的置信区间一起为确定研究药物的最佳剂量提供至关重要的信息。因此准确且充分的曲线及参数的估计信息对于量效关系的研究来说非常重要。用于描述药物量效关系的模型大多是非线性模型,与线性模型相比,非线性量效关系模型在参数及置信区间的估计上均有较大的难度,尤其对于参数置信区间的估计,在目前的文献中采用不正确或者不准确的估计方法并不少见,而这一点并没有引起研究者足够的重视。对于非线性模型参数置信区间的估计,轮廓似然法(Profile likelihood approach)是被推荐使用的一种方法。相比于常用的Wald法,该方法可以提供较好的置信区间覆盖率,也可以避免结果的偏差和错误,而且在逐点置信区间估计时也考虑到了模型的单调性。在某些情况下轮廓似然法的表现还优于Bootstrap方法。但目前该方法在非线性模型量效关系的研究中使用并不常规,而且在对曲线逐点置信区间估计上的应用更少。其主要原因在于轮廓似然法的计算通常来说会比较复杂,大多需要采用优化算法,模型无法收敛的情况比较多见。且对于复杂的非线性模型以及需要计算曲线上多点效应置信区间时,计算过程中不收敛的情况发生频率更高,计算也更为挑战,所以导致轮廓似然法在该领域虽然被推荐但目前仍不被经常使用。以往的研究主要是集中在轮廓似然法对于非线性模型单个参数置信区间的估计以及相对于其它方法的表现的评价,而对于如何基于轮廓似然法计算非线性模型整条曲线逐点效应的置信区间,目前还没有直接相关的研究报道。考虑到轮廓似然法在非线性模型参数置信区间的估计上有很显著的优点,且曲线逐点置信区间在量效关系的研究中会提供非常关键的决策信息,我们认为有必要研究采用该方法在非线性量效关系模型上对整个曲线逐点效应置信区间估计的实现,以促进该方法在这一领域更为常规的应用。另外,由于非线性模型的复杂性,模型拟合过程发生不收敛或收敛至一个不能接受的异常值的情况时有发生,而且较为常见。在对非线性模型参数估计时可以考虑给参数设定界值以改善模型的收敛情况,但这时对于达到设定界值的参数则无法采用常规的方法进行估计。根据以往的研究报道,发生模型不收敛或者参数达到设定界值而无法进行参数估计时,贝叶斯方法是一个可以考虑的方法。但该方法在不同试验设计和不同样本量下,尤其是对灵活性比较高的非线性量效关系模型的参数估计及评价并无相关研究报道。针对以上问题,本研究从以下两个方面展开研究,主要研究内容包括:
  1.基于轮廓似然法非线性量效关系曲线逐点置信区间的计算。
  本研究给出了一种简单易行且结果稳健的基于轮廓似然法估计非线性模型逐点效应置信区间的算法,并开发了相应的SAS宏程序%pwcl。该算法是基于一项二分数据在灵活性较高的4-参数logistic模型上开发的。主要是采用合适的计算顺序应用二分法并当期望反应接近边界时联合搜索法进行计算的方法。该算法考虑到非线性模型逐点置信区间上下限具有单调性的特点,对于各个剂量点在采用二分法设定初始搜索区间上下限时加以应用,使得该方法的计算速度得以提高。并且考虑到发生参数变换的非线性模型在拟合时需要给出各个参数更为合适的迭代初始值以促进模型收敛的情况,在对于某些剂量点采用二分法难于估计的情况下,结合搜索法以方便利用相邻拟合中得到的参数的极大似然估计值,进而为当前迭代提供合适的参数初始值。这些考虑使得在基于轮廓似然法计算非线性模型逐点置信区间的过程中模型收敛率大大提高。我们采用计算机模拟研究对该算法在不同的样本量下进行了测试,该算法在大多数的情况下都相对容易实现并且表现良好,只有在极个别小样本量的模拟数据上计算没有成功。另外也与文献报道中采用的轮廓似然法的算法对曲线逐点置信区间或者单个参数置信区间的计算进行了比对,结果均基本完全一致。与其它常用的方法如Wald法和Bootstrap法相比,对于模拟研究中设定的剂量点,基于轮廓似然法构建的置信区间的覆盖率都更为接近名义水平,尤其是在较低剂量效应接近0%时。这一好的表现在不同的样本量下,包括在样本量较小时都是一致的。本研究中在正态分布的数据上采用Emax模型测试了该算法,对我们编写的SAS宏程序只需进行简单的修改就可以使用,非常方便。所以从应用的角度考虑,本研究中提出的基于轮廓似然法计算非线性量效关系曲线逐点置信区间的算法相对比较简单,结果稳健,具有通用性。该算法以期可以促进轮廓似然法在非线性量效关系模型参数置信区间估计时得到更多的应用。
  2.非线性量效关系模型参数无法估计时贝叶斯方法的研究与应用。
  本研究主要采用计算机模拟的方法,在6种试验设计、7种不同样本量以及3种不同先验信息设定下对贝叶斯方法在灵活性较高的4-参数logistic模型上对单个参数的估计进行了研究。首先在模拟研究中对非线性模型不收敛或者参数达到设定界值时的情况进行了研究和评价。对于本研究采用的模拟数据以及非线性量效关系模型,采用常规频率学派的方法对模型参数进行估计时,对于某些参数会由于模型不收敛或者拟合时达到设定的界值的比例非常高而无法被估计。根据本研究结果,在对这些无法采用常规方法进行估计的参数,贝叶斯方法仍然可以给出一个较好或者尚可的估计,而且通过改善试验设计,或者增加有信息先验分布的设定,可以改善贝叶斯方法构建的可信区间的覆盖率。对于可以采用常规方法进行估计的参数,贝叶斯方法也可以给出一个和常规方法结果相当的估计,而且这种情况下有的参数即使在无先验信息的条件下也可以给出较好的估计。整个研究中,贝叶斯方法在对参数估计的过程中没有明显的数值问题,即使是对试验设计不佳以及先验分布信息不充分的情况,估计过程也没有发现明显的问题,比较稳健。另外对于基于轮廓似然法计算非线性模型逐点置信区间失败的情况下,贝叶斯方法也可以用于计算曲线逐点可信区间。和Bootstrap方法相比,低剂量附近的可信区间覆盖率不如Bootstrap法,不过在其它剂量处表现尚可,而且得到曲线逐点可信区间相较于Bootstrap方法更为平滑。总之,本研究的结果表明在探索非线性量效关系时,如果常规的方法在参数估计时不适用,贝叶斯方法是一个较好的可以替代常规方法用于非线性模型参数估计的方法。
  本研究的主要创新点有:1)提出了一种基于轮廓似然法计算非线性模型曲线逐点置信区间的算法,并开发了相应的SAS宏程序。该算法简单易实现,且结果稳健,在其它类型的数据及非线性量效关系模型上具有通用性;2)对非线性模型拟合不收敛或者参数达到设定界值的情况进行了研究和评价;3)探讨和研究了非线性模型参数采用常规方法无法估计时,贝叶斯方法的应用及相应参数估计时的相关问题,以及对于轮廓似然法在计算曲线逐点置信区间失败时贝叶斯方法的应用,并对贝叶斯方法的应用条件提供推荐。
  综上所述,本研究结合目前非线性模型参数及置信区间估计的方法以及现有的问题展开研究,提出了一种基于轮廓似然法计算非线性模型逐点效应置信区间的算法,以及在常规频率学派的方法对非线性模型参数无法估计时,对贝叶斯方法在非线性模型参数估计中的应用评价及其相关问题进行了探讨。研究结果可以为非线性量效关系模型参数及曲线置信区间的估计在新药临床剂量探索研究中的实际应用提供参考。
其他文献
煤岩声发射信号中包含了大量的煤矿动力灾害前兆信息,通过对其进行检测与定位研究,可实现煤矿井下动力灾害预警,对保障煤矿安全生产和人身安全具有重要意义。但由于煤矿井下强机械噪声的干扰,导致煤岩声发射信号发生畸变甚至淹没在噪声中,严重影响了信号检测与声源定位。  以该问题为背景,本文对强机械噪声干扰下煤岩声发射信号检测与声源定位算法开展了系统的研究。主要研究内容如下:  (1)研究了强机械噪声环境下煤岩
学位
报警系统对现代工业的安全生产与高效运行具有十分重要的作用。然而,现代工业中的报警系统普遍存在着干扰报警过多的问题,严重影响了报警系统的性能。现代工业报警系统经常采用报警死区减少干扰报警的数量。但是,报警死区在实际应用中面临着两个尚未解决的基本问题:如何判断某个信号是否适合使用报警死区?如何设计最优报警死区宽度?本文的研究目标就是初步解决这两个基本问题,从而针对合适的信号,使用相应的报警死区,减少干
学位
热固性树脂基复合材料具有许多优异的性质,如耐高温、硬度大、耐化学试剂、尺寸稳定性好等,广泛应用于防护涂料、胶粘剂、浇铸材料、增强塑料等领域。热固性树脂基复合材料的性能与其分子结构密切相关,而分子结构在很大程度上取决于树脂体系以及与固化过程有关的固化温度、固化时间、固化程度等因素。所以,为了开发高性能热固性树脂基复合材料,必须对热固性树脂的固化过程进行有效监测,研究树脂体系和固化工艺对树脂固化过程和
该论文的主要目的是采用固相合成反应合成新型的[AMM′Q]型多元金属硒化物.经EDS表征,合成了十九个新的[AMM′Q]型多元金属硒化物.测定了KGeSe和KCuInSe的晶体结构.通过测定KGeSe,KCuInSe的反射光谱,探讨了这两种晶体的光吸收特性.探讨了固相合成反应的规律.
学位
该文使用截流法和流动法着重考察了O和CS、COS在光照和无光照条件下的反应.无光照条件下CS和O的反应对于组分CS和O分别表观为一次,-d[O]/dt=k[CS][O].无光照条件下COS和O的反应也有-d[O]/dt=k[COS][O],但是相对于CS和O的反应而言,COS和O的反应非常缓慢.在红外光谱中读出CS和O的反应,在CS和O不同的浓度比例时,产物是不同的.当[O]]时,O和CS的反应产
学位
该论文首次报道了碳化钨负载的SO/ZrO(PSZ)固体超强酸催化剂(WC/PSZ),并研究了该催化剂上的正戊烷反应情况及影响催化剂活性的各种因素,与载Pt的SO/ZrO催化剂相比显示出更佳的催化性能,应用BET、XRD、XPS和FTIR等分析手段对催化剂进行了表征,得到以下结论:1.常温下PSZ固体超强酸催化剂对正戊烷具有异构化和裂解的双重作用,主要产物分别为异戊烷和异丁烷.2.高温条件下,PSZ
学位
大气环境是一个巨大、复杂、多因素相互作用的系统。PM10指空气动力学直径DP≤10μm的气体或固体颗粒物。几乎所有的大气污染物都会成为PM10污染的组成部分。由于大气污染中可吸入颗粒物PM10的特殊组成结构、理化性质以及对人体产生的越来越明显的伤害,PM10的监测与治理在整个大气环境整治中具有重要的地位与作用。但是,由于PM10在时间、空间与数量上的随机性和不确定性,增加了政府与企业如何选择有效的
学位
第一部分RIG-I样受体基因多态性与HCV感染慢性化的关联研究  丙型肝炎病毒(hepatitis C virus,HCV)感染严重危害全人类的健康。根据感染持续的时间,HCV感染可被分为急性和慢性感染。急性感染指感染者不经任何治疗即可在感染6个月之内自行清除病毒,通常无明显症状。仅15%-45%的感染者可以自限清除HCV,其余60%-80%的感染者会发展为慢性感染。固有免疫系统是识别病毒的第一道
铁是人体所必需的微量元素,参与血红蛋白、肌红蛋白、氧化呼吸链等重要生命物质的合成。铁稳态代谢调控对维持机体健康十分重要。铁缺乏会出现缺铁性贫血;而铁的蓄积会导致细胞死亡、组织损伤,甚至引发疾病,如遗传性血色病、地中海贫血等。铁调素是由肝脏分泌,具有调控机体铁代谢作用的核心激素,可以通过诱导小肠上皮细胞和巨噬细胞表面的铁外排转运蛋白降解来抑制机体铁吸收和铁循环,维持机体的铁稳态平衡。  铁稳态失衡,
学位
药物的全球同步研发,已经成为一种共享全球研发资源的开发模式,可以控制不必要的资源浪费,有助于解决在关键市场上市时间的滞后问题,使患者有更早的机会获得更新更好的治疗方案。国际多区域临床试验(Multi-regional Clinical Trials, MRCT)的目的在于多区域同步研发、同步注册上市,提高药物研发的效率,以利于更多地区的患者接受新的更有效的治疗方案。然而,由于不同区域之间存在种族、