讲一讲孟德尔随机化分析中的敏感性分析（sensitivity analyses）

很多人问过我敏感性分析的问题，说明大家对敏感性分析普遍是有疑问的，我之前的推文中也没有详细的讲过，因此决定写这一篇补充一下敏感性分析的工作。

什么是敏感性分析

敏感性分析定义：敏感度分析是研究数学模型或系统(数值或其他)输出中的不确定性如何在其输入中被分配到不同的不确定性来源 ——摘自维基百科。

我知道这段拗口的定义让人看了，感觉跟没看一样。

所有的文字都认识，但是连起来就读不懂是什么意思了。

所以我直白的翻译一下，敏感性分析就是通过一些手段鉴定结果是否靠谱的一项分析。

为什么要进行敏感性分析

前面回答了敏感性分析的定义，显而易见，它的作用包括以下三项：

1）评估结果是否稳健，结论是否靠谱；

2）评估结果是否有潜在的偏倚（比如基因多效性，数据异质性）；

3）评估是否存在某一个工具变量严重影响结局变量；

敏感性分析包括什么

广义的敏感性分析包括以下三项：

1）基因多效性检验；

2）异质性检验；

3）“leave-one-out”法；

侠义的敏感性分析指的是“leave-one-out”法。

现在大多数的孟德尔随机化研究提到敏感性分析的时候指的是“leave-one-out”法。

基因多效性和异质性检验会单独拎出来讲。

敏感性分析的结果怎么看

1）基因多效性检验；

关于基因多效性的文章我之前有讲过，见推文使用MR-PRESSO检验水平基因多效性（孟德尔随机化分析）。

我这里就不再赘述了。

2）异质性检验；

来自不同的分析平台、实验、人群等的工具变量（IV，一般指SNP）可能存在异质性，最终影响结果。

异质性检验我之前也写过的，见推文有相关性就有因果关系吗，教你玩转孟德尔随机化分析

命令很简单，就一句mr_heterogeneity(dat)

异质性检验的结果看P值，P值大于0.05说明纳入的 IVs 不存在异质性，研究结果不需要考虑异质性造成的影响。

3）“leave-one-out”法；

“leave-one-out”法指的是逐步剔除每个SNP，计算剩余SNP的meta效应，观察剔除每个SNP后结果是否发生变化，如果剔除了某一个SNP后，结果改变很大，说明存在某一个SNP对结果影响很大，这是我们不愿意看到的情况。

理想的情况应该是逐步剔除某个SNP后结果变化不大。

如下图所示，剔除每个SNP后，总体的误差线变化不大（所有的误差线均在0的右侧或者所有的误差线在0的左侧），说明结果是靠谱的。

image

这期就讲到这啦，如果你有什么疑问，欢迎进群交流（见公众号菜单栏）或者向我提问，我会结合大家的需求考虑是否发推文统一解答。

image