金志寧 張肇宇 馬俊奇 庫恩都斯·庫瓦提江 吳月康 森干

摘要：目的，建立個體腦卒中發(fā)病風(fēng)險評估模型。 方法? 對有關(guān)腦卒中發(fā)病危險因素的文獻(xiàn)進(jìn)行meta分析，計算個體腦卒中發(fā)病風(fēng)險值，隨機(jī)生成10000例模擬數(shù)據(jù)，確定風(fēng)險等級的分界值。結(jié)果，共納入28篇文獻(xiàn)，Hcy、肥胖不具有統(tǒng)計學(xué)意義（P>0.05），吸煙、體育鍛煉、腦卒中家族史等具有統(tǒng)計學(xué)意義。討論，已有的數(shù)據(jù)表明，致使腦卒中的多因素包括：吸煙、體育鍛煉、腦卒中家族史、高血壓、高尿酸、SBP等。

關(guān)鍵詞：腦卒中;危險因素;meta分析;評估模型

中圖分類號：TP311? ? ? ? 文獻(xiàn)標(biāo)識碼：A? ? ? ? 文章編號：1009-3044（2019）01-0240-03

腦卒中是一組以出現(xiàn)意識障礙和局灶性神經(jīng)功能缺失為共同特征的發(fā)病急促的腦血管慢性病[1]。 由于影響人群腦卒中發(fā)病的因素較多，針對不同的人群，同一種發(fā)病因素的強(qiáng)度也存在差異，因此長期的健康管理，可以讓健康人群了解可能引起腦卒中發(fā)生的因素，并通過對高危人群的健康風(fēng)險評估進(jìn)行發(fā)病預(yù)測，從而提供有效的預(yù)防控制措施。本研究運(yùn)用meta分析的方法對腦卒中發(fā)病風(fēng)險的影響因素進(jìn)行評價，并采用Logistic建立模型為腦卒中的預(yù)防提供依據(jù)[2]。

1? 資料與方法

1.1 收集系統(tǒng)評價的資料

1.1.1 文獻(xiàn)檢索

手工檢索萬方數(shù)據(jù)知識服務(wù)平臺、維普期刊資源整合服務(wù)平臺、中國生物醫(yī)學(xué)文獻(xiàn)服務(wù)系統(tǒng)、百度學(xué)術(shù)、PubMed相關(guān)會議論文，納入國內(nèi)外發(fā)表的關(guān)于腦卒中患者發(fā)病風(fēng)險的研究。英文檢索詞有：“risk factor”“ischemic”“stroke”“hemorrhagic”;中文檢索詞有：腦卒中、高血壓、危險因素、發(fā)病因素等。

1.1.2? 納入排除標(biāo)準(zhǔn)

納入標(biāo)準(zhǔn)：

（1） 研究對象：來源于社區(qū)、醫(yī)療單位或其他機(jī)構(gòu);

（2） 性別、年齡以及腦卒中類型不限;

（3） 以隊列研究或病例對照研究為研究類型;

（4） 腦卒中的診斷依據(jù)采用國際上通用的診斷標(biāo)準(zhǔn);

（5） 以發(fā)生腦卒中為觀測終點(diǎn)。

排除標(biāo)準(zhǔn)：

（1） 排除重復(fù)發(fā)表、存在嚴(yán)重缺失數(shù)據(jù)或數(shù)據(jù)報道不全的研究;

（2） 排除腦卒中復(fù)發(fā)的研究。

1.1.3 質(zhì)量評價和資料提取

對納入的文獻(xiàn)進(jìn)行資料提取，閱讀最終納入文獻(xiàn)，提取文獻(xiàn)基本信息，包括：第一作者、發(fā)表年份、研究地點(diǎn)、研究時間、病例來源等。同時提取引起腦卒中發(fā)病的相關(guān)因素如吸煙、體育鍛煉、高血壓、糖尿病、腦卒中家族史等[3]。

1.2? 統(tǒng)計學(xué)方法

1.2.1? meta分析

運(yùn)用Review Manager 5.3統(tǒng)計軟件對納入的28篇文獻(xiàn)提取出的與腦卒中發(fā)病有關(guān)的危險因素的RR值（某一危險因素的相對危險度）、95%CI（95%可信區(qū)間）等數(shù)據(jù)進(jìn)行meta分析。當(dāng)異質(zhì)性檢驗結(jié)果P≤0.05時，采用隨機(jī)效應(yīng)模型進(jìn)行meta分析;反之當(dāng)異質(zhì)性檢驗結(jié)果P>0.05時，采用固定效應(yīng)模型進(jìn)行meta分析，最終得到合并后的RR值和95%CI。比較兩種模型的結(jié)果，進(jìn)行敏感性分析，評價其穩(wěn)定性。最后采用Beeg檢驗，Egger檢驗判斷是否存在發(fā)表偏倚。

1.2.2? Logistic回歸模型

用Logistic回歸模型對疾病的發(fā)病情況進(jìn)行預(yù)測，選用優(yōu)化后的 Logistic回歸模型作為本次研究的基本模型，結(jié)合meta分析的結(jié)果，建立個體腦卒中發(fā)病風(fēng)險評估模型[5]。

公式（1）

其中P是個體發(fā)病風(fēng)險，[x1、x2、x3…xi…xn] 分別為第1、2、3…i…n個發(fā)病因素，[βi]是以[LnRRi]估計，常數(shù)項[α]是通過腦卒中發(fā)病率來估計的，根據(jù)Logistic回歸公式，

其中P為腦卒中發(fā)病率，[x1、x2、x3...xi、xn]，分別為第1、2、3…i…n個發(fā)病因素在人群中的平均暴露水平，但在實際情況下，一般很難得到，所以忽略，優(yōu)化后得到：

所以得到：

[Logit（P）=α+β1x1+β2x2+β3x3+...+βixi+βnxn]? 公式（3）? ?2? 結(jié)果

2.1 文獻(xiàn)納入結(jié)果

初檢出文獻(xiàn)232篇，對獲取的232篇文獻(xiàn)進(jìn)行查重，獲得195篇文獻(xiàn)。根據(jù)納入排除標(biāo)準(zhǔn)進(jìn)行初篩，排除65個研究。閱讀全文進(jìn)行復(fù)篩，排除102篇研究對象。最終納入28篇文獻(xiàn)[6-33]進(jìn)行meta分析，研究時間從2006-2016年。

2.2? meta分析結(jié)果

共納入9項因素，有7項因素的meta分析合并結(jié)果具有統(tǒng)計學(xué)意義（P<0.05），包括吸煙、體育鍛煉、腦卒中家族史、收縮壓（SBP）、高尿酸、規(guī)律服藥。合并RR值為1.957、0.387、1.878、2.575、1.797、2.328、2.702。此外，同型半胱氨酸（Hcy）和肥胖不具有統(tǒng)計學(xué)意義。因此未納入模型研究（見表1）。

2.3? Rothman-Keller模型參數(shù)

Rothman-Keller模型參數(shù)[4]計算方法如下：

（1）基準(zhǔn)發(fā)病比例：

Pi（某一危險因素在全人群中的暴露率），RRi（某一危險因素的相對危險度）。

（2）危險分?jǐn)?shù)

[S=ρ×RRi]

（3）組合危險分?jǐn)?shù)

[?=（P1-1）+...+（Pn-1）+q1×q2×...qn]

P為危險分?jǐn)?shù)≥1的項之值，q為危險分?jǐn)?shù)<1的項之值。

（4）個體發(fā)病風(fēng)險

個體發(fā)病風(fēng)險=人群腦卒中發(fā)病率 [×?]

通過研究獲取某一危險因素在全人群中的暴露率Pi;通過meta分析合并危險因素的風(fēng)險值，得到某一危險因素的相對危險度RRi。腦卒中預(yù)測模型的參數(shù)值見表2。

2.4? 確定等級劃分界值

本研究根據(jù)7個危險因素的二項分布函數(shù)，確定低發(fā)病風(fēng)險、中發(fā)病風(fēng)險、高發(fā)病風(fēng)險的分界值，一共生成10000例患者7個危險因素分布情況。將隨機(jī)生成的這10000例數(shù)據(jù)，代入Rothman-Keller模型，分別計算相應(yīng)的個體發(fā)病風(fēng)險。將發(fā)病風(fēng)險從小到大排序，以個體ID為橫坐標(biāo)，腦卒中風(fēng)險為縱坐標(biāo)繪制聯(lián)合分布圖。

根據(jù)風(fēng)險的變化趨勢，選取0.048為低發(fā)病風(fēng)險和中發(fā)病風(fēng)險的分界值，0.079為中發(fā)病風(fēng)險和高發(fā)病風(fēng)險的分界值。即腦卒中發(fā)病風(fēng)險<0.048，屬于低危人群;腦卒中發(fā)病風(fēng)險在0.048—0.079，屬于中危人群;腦卒中發(fā)病風(fēng)險>0.079，屬于高危人群。（見圖1）

2.5? 建立Logistic模型

有研究指出腦卒中的發(fā)病率為2.57%[3]，查找表2，可知各項危險因素的<E：＼知網(wǎng)文件＼電腦＼電腦01-02＼電腦01＼7xs201901＼Image＼image17.png>值，計算 <E：＼知網(wǎng)文件＼電腦＼電腦01-02＼電腦01＼7xs201901＼Image＼image18.png>值，代入公式（3），得最終的Logistic回歸模型：

3? 討論

本次研究對腦卒中的發(fā)病因素進(jìn)行了發(fā)病風(fēng)險值的合并，建立了我們自己的腦卒中個體發(fā)病風(fēng)險健康評估模型，并運(yùn)用二項分布函數(shù)進(jìn)行了10000例數(shù)據(jù)的隨機(jī)模擬，對腦卒中個體發(fā)病風(fēng)險進(jìn)行了等級劃分。

本研究發(fā)現(xiàn)致使腦卒中發(fā)病的主要影響因素有：吸煙、體育鍛煉、腦卒中家族史、收縮壓、高尿酸、規(guī)律服藥。可以通過加強(qiáng)體育鍛煉、戒煙戒酒、規(guī)律服藥等健康的生活方式或飲食習(xí)慣降低腦卒中的發(fā)病風(fēng)險。根據(jù)趨勢圖最終確定0.048為低發(fā)病風(fēng)險和中發(fā)病風(fēng)險的分界值，0.079為中發(fā)病風(fēng)險和高發(fā)病風(fēng)險的分界值。將腦卒中劃分為低危、中危、高危三個風(fēng)險等級，針對不同的人群，要進(jìn)行不同的健康教育。本研究結(jié)果可以作為腦卒中發(fā)病風(fēng)險預(yù)測的參考依據(jù)。

參考文獻(xiàn)：

[1] 吳智鑫，蔣龍元.B型腦鈉肽對急性缺血性腦卒中患者TOAST病因?qū)W分型的應(yīng)用進(jìn)展[J].嶺南急診醫(yī)學(xué)雜志，2012，17（5）：399-400.

[2] 周亮，程迪祥，張彥琦，等.腦卒中發(fā)病風(fēng)險logistic回歸綜合評估模型研究——基于亞洲人群研究數(shù)據(jù)[J].重慶醫(yī)科大學(xué)學(xué)報，2014，38（4）：446-450.

[3] 黃艷紅，周亮，陳虹汝，等.中國高血壓患者并發(fā)腦卒中風(fēng)險的Rothman-Keller模型研究[J].第三軍醫(yī)大學(xué)學(xué)報，2017，（20）：2042-2050.

[4] Rothman K，Keller A.The effect of joint exposure to alcohol and tobacco on risk of cancer of the mouth and pharynx[J].J Chronic Dis，1972，25（12）：711-716.

[5] 孫穎，楊君慧.Logistic回歸與分類樹模型的比較[J].西安工業(yè)大學(xué)學(xué)報，2014（9）：689-692.

[6] 來春林，邢金平，胡風(fēng)云，等.脂蛋白相關(guān)磷脂酶A_2與原發(fā)性高血壓患者發(fā)生缺血性卒中的相關(guān)性分析[J].中國腦血管病雜志，2016，13（09）：455-459.

[7] 陳銀燕，胡飛娥，陳榮浩，等.2012-2015年某院高血壓性腦出血的臨床特點(diǎn)及危險因素[J].中國老年學(xué)雜志，2016，36（11）：2748-2749.

[8] 閆麗，曹燁，汪秋實，覃玉，等.高血壓患者腦卒中發(fā)病影響因素Cox回歸分析[J].中國公共衛(wèi)生，2016，32（06）：748-752.

[9] 張金武.老年高血壓患者動態(tài)血壓參數(shù)與缺血性腦卒中的相關(guān)性[J].中國實用神經(jīng)疾病雜志，2016，19（06）：72-73.

[10] 龐慧，韓冰，李先池，等.性別差異對于原發(fā)性高血壓患者腦卒中風(fēng)險的影響[J].中華腦科疾病與康復(fù)雜志（電子版），2016，6（01）：14-19.

[11] 李曉乾，孔新興，王寶艷，等.高血壓性腦出血發(fā)病的相關(guān)因素分析[J].廣西醫(yī)學(xué)，2016，38（01）：43-46+49.

[12] 黃少芬，肖軍，胡崢，等.高血壓患者同型半胱氨酸與腦卒中的關(guān)系分析[J].實用預(yù)防醫(yī)學(xué)，2015，22（10）：1184-1188.

[13] 趙潔，張小培，鐘愛群，等.健康行為對原發(fā)性高血壓患者腦卒中發(fā)病影響的病例對照研究[J].中國全科醫(yī)學(xué)，2012，15（34）：3963-3967.

[14] 周潔.腦卒中患者危險因素分析及社區(qū)干預(yù)模式探討[J].中外醫(yī)學(xué)研究，2012，10（14）：28-29.

[15] 魏凱，張念民，吳天華.缺血性腦卒中發(fā)病因素的病例對照分析[J].西南軍醫(yī)，2012，14（03）：471-473.

[16] 宋晴，郭艷華，高彩霞，等.缺血性腦卒中發(fā)病危險因素的病例對照研究[J].職業(yè)與健康，2012，28（01）：4-6.

[17] 邢立影，張慧英，郭來敬，等.高血壓并發(fā)腦卒中的危險因素分析[J].中國醫(yī)藥導(dǎo)報，2011，8（33）：144-145.

[18] 薛蕾，唐義平.高原藏區(qū)腦卒中危險因素的相關(guān)性分析[J].實用醫(yī)院臨床雜志，2011，8（05）：62-64.

[19] 周國芳，王會明，劉志明，等.男性腦梗死患者首次發(fā)作的危險因素分析[J].山東醫(yī)藥，2011，51（29）：61-62.

[20] 常彬賓，肖欣榮，呂柏寧，等.高齡老年高血壓人群血尿酸濃度與腦梗死的相關(guān)性研究[J].中華保健醫(yī)學(xué)雜志，2011，13（03）：183-185.

[21] 王燦，時景璞，付凌雨，等.遼寧省彰武縣農(nóng)村地區(qū)血清脂聯(lián)素水平與腦卒中關(guān)系的病例對照研究[J].中華疾病控制雜志，2011，15（01）：11-14.

[22] 黃中紅，楊洪江，那開憲.高血壓患者并發(fā)腦卒中危險因素分析[J].現(xiàn)代中西醫(yī)結(jié)合雜志，2009，18（31）：3813-3814.

[23] 關(guān)紅軍，崔新剛，郭毓鵬，等.牡丹江市腦卒中危險因素的病例對照研究[J].牡丹江醫(yī)學(xué)院學(xué)報，2009，30（04）：16-18.

[24] 董兆舉，崔玉蘭，高梅蘭，等.農(nóng)村社區(qū)人口腦卒中發(fā)病危險因素的病例對照研究[J].中國農(nóng)村衛(wèi)生事業(yè)管理，2008（07）：514-515.

[25] 任雁京.貴陽市省直機(jī)關(guān)干部高血壓并發(fā)腦卒中的危險因素分析[J].中國臨床保健雜志，2007（04）：394-395.

[26] 黃久儀，王桂清，郭吉平，等.腦血管血流動力學(xué)積分對高血壓患者卒中相對危險度的評估[J].中國腦血管病雜志，2006（12）：537-541+576.

[27] 王娓娓，張紅苗，張瑩，等.老年高血壓患者高尿酸血癥與缺血性腦卒中相關(guān)性分析[J].中華老年醫(yī)學(xué)雜志，2012，（6）：462-464.

[28] 楊艷，肖玉嫻.高血壓腦出血的危險因素分析及治療體會[J].中國實用神經(jīng)疾病雜志，2015，（19）：72-73，74.

[29] 呂鴻杰.高血壓腦出血相關(guān)危險因素影響分析[J].山西醫(yī)藥雜志，2016，45（14）：1661-1663.

[30] 楊亞莉，李艷紅，劉艷瓊，等.高血壓患者頸動脈斑塊和急性缺血性腦卒中的相關(guān)性研究[J].中國實用神經(jīng)疾病雜志，2015，18（18）：11-13.

[31] 常旺，崔建忠，王凱杰，等.高血壓性腦出血臨床危險因素分析[J].河北聯(lián)合大學(xué)學(xué)報（醫(yī)學(xué)版），2015，17（05）：100-103.

[32] 韋彤雁，蒙蘭青，李天資，等.高血壓伴腔隙性腦梗死患者危險因素關(guān)聯(lián)性分析[J].右江醫(yī)學(xué)，2015，43（04）：401-404.

[33] 王純，梁雨露.高血壓患者頸動脈斑塊和急性缺血性腦卒中的相關(guān)性[J].貴陽醫(yī)學(xué)院學(xué)報，2014，39（1）：80-83.