论文部分内容阅读
在数十年的全基因组关联研究热潮中,大量与人类复杂性状相关联的单核苷酸多态性位点被成功挖掘。然而,这些关联位点仅能解释性状的很少一部分遗传度,我们称之为“遗传度缺失”(missing heritability)现象。该现象的一个合理解释是通用的单位点检验方法较难发现弱效应位点。多位点模型通常具有较高的统计学效能,而大规模高维组学数据下构建多位点模型则面临挑战。第一部分,在适应性条件双指数先验的框架下,构建基于变分推断算法的贝叶斯多位点模型(BAL-VI)。通过适应性条件双指数先验达到对位点效应适应性惩罚的目的;通过变分推断算法达到高维数据下快速、准确进行贝叶斯后验计算的目的。模拟试验中评价了模型在变量选择、参数估计和结局预测方面的表现,并就几个重要问题作探讨:模型区间估计宽度、高相关结构下变分算法与MCMC算法的估计性质以及超参数敏感性分析。实例分析将模型应用于肺癌GWAS数据。模拟试验结果表明,(1)BAL-VI的95%贝叶斯可信区间狭窄,基于该指标的变量选择产生大量的假阳性,借助于额外的遗传度阈值可以有效控制假阳性;(2)综合变量选择、参数估计和预测结果,BAL-VI表现最佳;(3)在高相关结构下,BAL-VI可以避免效应被“稀释”的风险,更有利于阳性位点及其所在区域的识别;(4)模型对超参数取值敏感。实例分析结果表明,(5)模型能够应用于全基因组尺度的数据,并识别多个与肺癌发生风险相关联的位点;(6)模型的速度优势明显,BAL-VI完成分析大约需要半天时间,而对应的MCMC算法则需要5天。第二部分讨论如何将线性EMVS(EM algorithm for Bayesian variable selection)模型推广到Weibull参数生存模型(SurvEMVS)。通过对效应参数施加连续性spike-and-slab先验,促进贝叶斯模型的变量选择;采用EM算法进行快速准确的后验空间探索和参数估计;针对高维下无闭型解(closed-form solution)的问题,算法中嵌套一个改进的循环坐标下降方法来快速更新效应估计值;应用EBIC(extended Bayesian information criterion)指标来选择高维模型下的最佳超参数。通过模拟试验评价模型的变量选择、效应估计和预测准确度,以及不同τ参数下EBIC指标对模型拟合的影响,此外还就不满足Weibull分布假设情形做了分析和讨论。实例分析将模型应用于带有生存结局的肺癌GWAS数据和胃癌基因表达数据。模拟试验结果表明,(1)综合变量选择、参数估计和生存预测的结果,EBIC(τ=0.5)指标下的SurvEMVS优于其它模型;(2)τ越大,越容易产生稀释模型;(3)SurvEMVS在适度偏离分布假设的情况下依然表现稳健。实例分析结果表明,(4)SurvEMVS发现了多个影响癌症预后的潜在位点,其中数个位点通过外部数据验证,表明该模型可以应用于基因组和转录组学数据;(5)模型的EM算法收敛速度快,仅需要几十次到数百次迭代。第三部分对现有研究进行总结,并对未来的研究方向作展望。