與職業相關的心腦血管疾病之族群可歸因分率與疾病負荷推估

一、 研究工具與分析對象

（一）計算職場危害暴露盛行率部分

由勞動部勞研所於 2016 年進行的「勞動環境安全衛生狀況認知調查」獲得 所需資料，該調查之抽樣方式、受試者特質等詳細描述，參見前節第一項「研究 工具與分析對象」之說明。本研究僅納入 20-65 歲受僱於政府或私人企業的工作 者進行分析。

（二）推估疾病負荷部分

串連 2000-2016 年全民健保資料庫的處方及治療明細檔（含門急診及住院）

找出目標疾病之新發個案。本研究的心腦血管疾病分為四類，視資料登錄的狀況 以 ICD9 或 ICD10 定義：中風 ICD9 為 434、ICD10 為 I63，心肌梗塞 ICD9 為 410 及 411（除 411.81）、ICD10 為 I21 及 I22，心臟衰竭 ICD9 及 ICD10 分別為 411.81、

412、414 及 I24.0、I25.0、I25.10、I25.2-6，心絞痛 ICD9 及 ICD10 分別為 413 及 I20、I23.7、I25.11、I25.7-9，其中，中風及心肌梗塞的新發個案定義同一個被保 險人自 2000 年起觀察到的第一筆住院紀錄，之後每年追蹤一次，若在年中（7 月 1 日）前 28 天有住院紀錄，則定義為當年度急性期的盛行個案，若無則定義 為慢性期盛行個案，在心肌梗塞的部分，急性期再分為急性初期（第 1、2 天）

及急性後期（第 3-28 天）；心臟衰竭及心絞痛的新發個案則是同一個被保險人自 2000 年起一年內門診檔任一診斷碼出現 2 次以上、或住院檔主診斷碼出現 1 次 以上的第一筆紀錄，之後追蹤至 2016 年，若存活則為本研究之盛行個案。上述

新發個案的年齡及身份別將連結到發病當年的全民健保資料庫承保檔做確認，本 研究僅納入 20-65 歲、第一類一至三目非眷屬身份之被保險人為分析對象。

二、 職場危害暴露與心腦血管疾病相關性之估計

受限於我國並無同時測量職場危害暴露與心腦血管疾病狀況的全國性資料，

無法直接建立模型計算臺灣本土的職場危害暴露與心腦血管疾病之相關性，因此 有別於直接計算職場危害暴露與肌肉骨骼疾病相關性的研究方法（前節），在心 腦血管疾病的部分，是以文獻回顧的方式獲得各職場危害暴露對心腦血管疾病之 RR。首先利用 PubMed 資料庫輸入關鍵字 1) occupational 或 work-related 或 job 或 worker*或 employee*、2) 限定 MeSH Term 為 cardiovascular diseases 或 cerebrovascular disorders、3) meta-analysis，取三者交集蒐集 2005 年之後的文獻，

再依據摘要內容排除非統合分析、目標疾病非心腦血管疾病、非討論暴露與疾病 之相關性、暴露非職場社會心理危害之研究後，若有二篇以上針對相同職場社會 心理危害項目進行分析的文獻，則挑選對職場社會心理危害項目有明確定義，且 使用世代研究資料、有校正共變數再進行統合分析的結果者作為後續計算 PAF 的參數；若僅有一篇文獻則直接將其分析結果納入進行後續計算。

三、 統計方法

（一）職場危害暴露盛行率之計算

依據前節職場社會心理危害之測量方法（頁 44-46），獲得職場社會心理危害 暴露之資料；職場噪音暴露的部分，也透過 2016 年「勞動環境安全衛生狀況認

知調查」的問卷內容獲得：將勾選工作環境「經常有」聲音很大同時「沒有」在 工作中「經常使用」耳塞、耳罩的受僱者定義為有職場噪音的暴露，其他受僱者 定義為沒有職場噪音暴露。接著，依據前段文獻回顧的結果，將調查資料中的職 場危害暴露調整為與文獻相同的定義，最後再利用百分率，將納入統計之「勞動 環境安全衛生狀況認知調查」分析對象的行業別做分層，呈現臺灣受僱者的職場 危害暴露情形，並使用卡方檢定職場危害暴露盛行率在主要行業別（製造業、營 造業、服務業）的差異。以上統計使用 SAS 9.4 版本計算。

（二）族群可歸因分率之計算

計算 PAF 的基本公式如（公式三），其中，RR 由本節第二段挑選出來的統合 分析研究結果代入，P 則以前段計算之暴露盛行率代入。

（公式三）

在使用斯皮爾曼等級相關係數確認各危險因子的獨立性之後，使用（公式五）

計算出聯合所有職場危害因素的 PAF_joint。

（公式五）

同時，本研究除了估算整體受僱者職場社會心理危害對心腦血管疾病的 PAF 之外，亦估算各主要行業別（製造業、營造業、服務業）的 PAF。以上統計使用 Microsoft Office Excel 2007 版本計算。

（三）疾病負荷之計算

疾病負荷的程度以 YLD 呈現，而計算 YLD 時所需的疾病嚴重度分布占率及各

嚴重度之失能權重參考 GBD 的參數(GBD 2017 Risk Factor Collaborators, 2018a)，

其中，因為健保資料庫可根據 GBD 對疾病嚴重度的定義區分出急性初期及急性 後期的心肌梗塞個案，所以心肌梗塞的 YLD 將直接使用由健保資料庫計得的疾病 個案數乘以 GBD 提供的失能權重計算，其他心腦血管疾病因無法根據 GBD 的定 義從健保資料庫中區分不同疾病嚴重度的個案，故將 GBD 提供的各疾病嚴重度 的占率及失能權重相乘後加總成為該疾病分類的 YLD（附表 4）；另外，有別於肌 肉骨骼疾病的疾病嚴重度分布及失能權重，心腦血管疾病的部分 GBD 並未提供 上、下限的估計。最後，各心腦血管疾病的 YLD 分別乘上由全民健保資料庫計算 出的相對應個案人數，再乘以前段估算的整體 PAF，即為與職業相關之心腦血管 疾病的 YLD；另外，也此 YLD 估計值除以當年度（2016 年）20-65 歲的年中人口 數，以方便利用每十萬人口的 YLD 做疾病負荷的比較。以上統計使用衛生福利部 統計處衛生福利資料科學中心提供之 SAS 軟體，以及 Microsoft Office Excel 2007 版本計算。

以上有關計算（公式三）PAF 及疾病負荷 YLD 的參數來源，整理於表 3-2。

表 3-2、計算與職業相關的心腦血管疾病之 PAF 及 YLD 的參數資料來源對照表

參數 資料來源

PAF 職場危害暴露盛行率(P) 勞動部「勞動環境安全衛生狀況認知調查」

職場危害暴露與心腦血 管疾病的相關性(RR)

利用 PubMed 資料庫蒐集 2005 年後的文獻資 料，以有校正共變數之統合分析的 RR 為主

YLD 盛行個案人數 全民健保資料庫「處方及治療明細檔」

疾病嚴重度分布 GBD 2017 Risk Factor Collaborators (2018a) 疾病失能權重 GBD 2017 Risk Factor Collaborators (2018a)

統整以上，本研究推估與職業相關的肌肉骨骼及心腦血管疾病之分析流程，

繪製如圖 3-2：

資料來源 survey¹ survey

or^reference

(MSDs²) (CCVD³) NHIRD⁴ GBD⁵

參 數 P⁶ RR⁷ case⁸ %⁹ DW¹⁰

與職業相關的

疾病負荷推估 PAF YLD

圖 3-2、職業相關肌肉骨骼及心腦血管疾病負荷推估流程圖

【註解】1. survey: 勞動部「勞動環境安全衛生狀況認知調查」；2. MSDs: 肌肉骨骼疾病；

3. CCVD: 心腦血管疾病；4. NHIRD: 全民健保資料庫；5. GBD: 全球疾病負荷；

6. P: 職場危害暴露盛行率；7. RR: 職場危害暴露與心腦血管疾病的相關性；8.

case: 盛行個案人數；9. %: 疾病嚴重度分布；10. DW: 疾病失能權重