統計方法

（一）個人因素之分布與職場危害暴露盛行率和疾病盛行率的計算

以納入統計之「勞動環境安全衛生狀況認知調查」分析對象的性別做分層，

使用百分率呈現受僱者年齡及行業別的分布，再以分析對象之性別與年齡、性別 與行業做分層，使用百分率呈現受僱者個人因素、職場危害暴露之分布和疾病盛 行率的狀況，同時，亦使用卡方檢定年齡和行業別在性別上的差異，以及個人因

素、職場危害暴露盛行率和疾病盛行率在各年齡和各行業分層中性別上的差異，

最後，以性別、年齡、行業別做分層，計算各分層的各項職場危害暴露與脖子、

肩膀、下背或腰部的肌肉骨骼疾病盛行率。以上統計使用 SAS 9.4 版本計算。

（二）職場危害暴露與肌肉骨骼疾病相關性之估計

在雙變項分析的部分，納入主要行業（製造業、營造業、服務業）的分析對 象，以其性別及年齡分層，先行使用卡方兩兩檢定各危險因子對脖子、肩膀、下 背或腰部肌肉骨骼疾病症狀的相關性，若雙變項分析的統計結果為具有統計上顯 著意義的相關性（p 值<0.05），則再投入多變項羅吉斯迴歸中了解危險因子相互 校正後相關性的變化，若某危險因子在多變項羅吉斯迴歸時於大部分的分層統計 結果為不顯著，則排除投入以對數為連結函數的二項式模型（Log-linked binomial model）作參數；本研究將對脖子、肩膀、下背或腰部分別建立以對數為連結函 數的二項式模型，此三個推估模型亦以性別及年齡分層，投入一樣的參數計算各 危險因子對肌肉骨骼疾病的 PR。以上統計使用 SAS 9.4 版本計算。

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

計算 PAF 的基本公式如（公式三）。本研究以前段計算出的 PR 代替 RR 進行 PAF 的估計，若估算的危險因子暴露分為高、低二類，則 i 僅一層，直接使用高 暴露組的暴露盛行率及 PR 進行計算；若危險因子的暴露不只有高、低二類，如 工作時間分為低工時、正常工時、高工時三類，則 i 有二層（以正常工時作為沒 有暴露），計算 PAF 時就需要對應相同暴露程度的 PR 及 P 進行計算。

（公式三）

另外，參考 Miettinen (1974) 在計算 PAF 時所設的條件（RR≥1），當 PR 在統 計上不具顯著意義時，將以「1」代入公式中進行計算。

由於本研究的 PR 是經性別分層、控制年齡而來，因此針對每一個性別、年 齡的分層，使用（公式四）進行校正，例如：計算 20-25 歲男性工作心理負荷對 脖子肌肉骨骼疾病的 PAF，除了使用（公式三）帶入 20-25 歲男性工作心理負荷 對脖子肌肉骨骼疾病的 PR 及暴露於「高」工作心理負荷的 P 之外，亦代入 20-25 歲男性脖子部位有肌肉骨骼疾病症狀的個案占整體男性個案的百分率（ω20-25） 進行調整，而整體男性受僱者在工作心理負荷對脖子肌肉骨骼疾病的 PAFoverall， 即為校正各年齡分層肌肉骨骼疾病個案占整體個案的百分率之後 PAF 的總和。

（公式四）

最後，在使用斯皮爾曼等級相關係數（Spearman’s rank correlation coefficient，

以下以 r_s標示）確認各危險因子的獨立性之後，使用（公式五）計算出聯合所有 職場危害因素的 PAFjoint。

（公式五）

由於行業別與受僱者的職場危害暴露風險和罹患肌肉骨骼疾病的風險有關 (Eurofound, 2017; Marcum & Adams, 2017; Park et al., 2018; X. Wang et al., 2017)，

因此除了整體受僱者的 PAF 估算之外，本研究將再以行業分層，分別計算主要行 業別（製造業、營造業、服務業）各自的 PAF。以上統計使用 Microsoft Office Excel 2007 版本計算。

（四）疾病負荷之計算

疾病負荷的程度以 YLD 呈現。其中，計算 YLD 時所需的疾病嚴重度分布占率 及各嚴重度之失能權重參考 GBD 的參數(GBD 2017 Risk Factor Collaborators, 2018a)，頸部痛的疾病嚴重度分布與各嚴重度之失能權重整理於附表 1，將各疾 病嚴重度的占率及失能權重相乘後加總，即為本研究代入之頸部痛的 YLD 數值；

由於 GBD 並無獨立「肩膀痛」的推估，因此參考「其他肌肉骨骼疾病」的疾病 嚴重度分布及失能權重（整理於附表 2），轉換肩膀痛的疾病嚴重度分率（轉換 後無上下界），再使用與頸部痛相同的計算方式計算代入本研究之肩膀痛的 YLD 數值；下背痛的部分是以年齡分層呈現（附表 3），每個年齡層的各疾病嚴重度 之 YLD 加總後成為本研究各年齡層的下背痛代入之 YLD 數值。各部位肌肉骨骼疾 病的 YLD 分別乘上由全民健保資料庫計算出的相對應個案人數後，再乘上前段估 算之整體的 PAF，作為職場危害暴露對與職業相關的肌肉骨骼疾病之疾病負荷程 度。最後，將男、女性因職場危害造成的肌肉骨骼疾病之 YLD 估計值分別除以當 年度（2016 年）20-65 歲的年中人口數，以方便利用每十萬人口的 YLD 做疾病負 荷的比較。以上統計使用衛生福利部統計處衛生福利資料科學中心提供之 SAS 軟體，以及 Microsoft Office Excel 2007 版本計算。

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

表 3-1、計算與職業相關的肌肉骨骼疾病之 PAF 及 YLD 的參數資料來源對照表

參數 資料來源

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

職場危害暴露與肌肉骨 骼疾病的相關性(RR)

勞動部「勞動環境安全衛生狀況認知調查」，

以計算暴露對疾病的 PR 代替 RR

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

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