很多人問過我敏感性分析的問題，說明大家對(duì)敏感性分析普遍是有疑問的，我之前的推文中也沒有詳細(xì)的講過，因此決定寫這一篇補(bǔ)充一下敏感性分析的工作。
什么是敏感性分析
敏感性分析定義：敏感度分析是研究數(shù)學(xué)模型或系統(tǒng)(數(shù)值或其他)輸出中的不確定性如何在其輸入中被分配到不同的不確定性來源 ——摘自維基百科。
我知道這段拗口的定義讓人看了，感覺跟沒看一樣。
所有的文字都認(rèn)識(shí)，但是連起來就讀不懂是什么意思了。
所以我直白的翻譯一下，敏感性分析就是通過一些手段鑒定結(jié)果是否靠譜的一項(xiàng)分析。
為什么要進(jìn)行敏感性分析
前面回答了敏感性分析的定義，顯而易見，它的作用包括以下三項(xiàng)：
1）評(píng)估結(jié)果是否穩(wěn)健，結(jié)論是否靠譜；
2）評(píng)估結(jié)果是否有潛在的偏倚（比如基因多效性，數(shù)據(jù)異質(zhì)性）；
3）評(píng)估是否存在某一個(gè)工具變量嚴(yán)重影響結(jié)局變量；
敏感性分析包括什么
廣義的敏感性分析包括以下三項(xiàng)：
1）基因多效性檢驗(yàn)；
2）異質(zhì)性檢驗(yàn)；
3）“l(fā)eave-one-out”法；
俠義的敏感性分析指的是“l(fā)eave-one-out”法。
現(xiàn)在大多數(shù)的孟德爾隨機(jī)化研究提到敏感性分析的時(shí)候指的是“l(fā)eave-one-out”法。
基因多效性和異質(zhì)性檢驗(yàn)會(huì)單獨(dú)拎出來講。
敏感性分析的結(jié)果怎么看
1）基因多效性檢驗(yàn)；
關(guān)于基因多效性的文章我之前有講過，見推文使用MR-PRESSO檢驗(yàn)水平基因多效性（孟德爾隨機(jī)化分析）。
我這里就不再贅述了。
2）異質(zhì)性檢驗(yàn)；
來自不同的分析平臺(tái)、實(shí)驗(yàn)、人群等的工具變量（IV，一般指SNP）可能存在異質(zhì)性，最終影響結(jié)果。
異質(zhì)性檢驗(yàn)我之前也寫過的，見推文有相關(guān)性就有因果關(guān)系嗎，教你玩轉(zhuǎn)孟德爾隨機(jī)化分析
命令很簡(jiǎn)單，就一句mr_heterogeneity(dat)
異質(zhì)性檢驗(yàn)的結(jié)果看P值，P值大于0.05說明納入的 IVs 不存在異質(zhì)性，研究結(jié)果不需要考慮異質(zhì)性造成的影響。
3）“l(fā)eave-one-out”法；
“l(fā)eave-one-out”法指的是逐步剔除每個(gè)SNP，計(jì)算剩余SNP的meta效應(yīng)，觀察剔除每個(gè)SNP后結(jié)果是否發(fā)生變化，如果剔除了某一個(gè)SNP后，結(jié)果改變很大，說明存在某一個(gè)SNP對(duì)結(jié)果影響很大，這是我們不愿意看到的情況。
理想的情況應(yīng)該是逐步剔除某個(gè)SNP后結(jié)果變化不大。
如下圖所示，剔除每個(gè)SNP后，總體的誤差線變化不大（所有的誤差線均在0的右側(cè)或者所有的誤差線在0的左側(cè)），說明結(jié)果是靠譜的。
這期就講到這啦，如果你有什么疑問，歡迎進(jìn)群交流（見公眾號(hào)菜單欄）或者向我提問，我會(huì)結(jié)合大家的需求考慮是否發(fā)推文統(tǒng)一解答。
