基于马尔可夫链蒙特卡罗算法的我国人群血铅和血镉水平估计研究.pdfVIP

基于马尔可夫链蒙特卡罗算法的我国人群血铅和血镉水平的估计研究

摘要

研究背景：

铅、镉等重金属污染已成为我国重要的公共卫生问题。为了深入探究一般人群

因重金属污染导致的健康风险，国内外相继开展了重金属人体监测项目（Human

biomonitoring，HBM），旨在分析一般人群血液、尿液等生物样本中的重金属水平

的空间分布和时间变化趋势。然而，这些监测项目往往无法追溯项目建立前的历史

数据，并且数据更新也常常存在滞后性。因此，采用重金属文献合并的方法可以弥

补这些不足。

既往研究中，通过公式法（FormulaMethod，FM）和蒙特卡罗模拟方法

（MonteCarloSimulation，MCS）进行文献的模拟及合并，但在模拟偏态分布数据

时，这些方法仍存在偏差。因此，有必要研究新的文献合并方法，以解决既往方法

存在的问题，展现更加接近真实的中国一般人群血铅、镉水平。

研究目的：

1.在蒙特卡罗模拟方法基础上，提出应用于人体生物材料中重金属水平估计的

马尔可夫链蒙特卡罗算法（MarkovChainMonteCarlo，MCMC），以一般人群血铅、

血镉水平为例，对偏态分布数据准确模拟。

2.应用MCMC算法，模拟一般人群的血铅水平、血镉水平原始数据，合并后展

现更加接近真实的一般人群血铅、血镉的时间和空间分布特征。

3.基于MCMC算法模拟的人群血铅、血镉数据，评估无铅汽油政策及《重金属

污染综合防治“十二五”规划》对降低人群中铅、镉重金属负荷的效果。

研究方法：

1.方法构建：基于Metropolis-Hastings算法原理，使用R语言构建适用于人群

重金属水平估计的MCMC算法。模拟数据通过Excel合并，并计算几何均值

（GeometricMean,GM）与几何标准差（GeometricStandardDeviation，GSD）。

2.方法验证：首先，使用马尔可夫链轨迹图、Geweke检验以及将变异系数

（CoefficientofVariation，CV）小于5%作为标准，分别评估MCMC算法的收敛性

和可靠性，并确定最佳模拟次数。其次，利用中国八省份重金属实测数据库检验

MCMC算法的准确度。

3.方法应用：在PubMed

般人群血铅、血镉的相关数据。以无铅汽油政策和《重金属污染综合防治“十二五”

规划》的颁布年份为时间节点，分析比较一般人群血铅、血镉的时空变化趋势。

研究结果：

1.通过马尔可夫链轨迹图和Geweke检验的双重验证，证明MCMC算法具有良

好的收敛性。当血铅数据模拟达到20次，血镉数据模拟达到50次时，其CV值均

低于5%，表明MCMC算法可靠，可作为血铅与血镉数据的参考模拟次数。

2.MCMC算法在模拟对数正态分布数据时相对误差小于5%；模拟样本量较小

的偏态分布数据时相对误差仍能保持在5%以内；即使偏态分布数据的标准差较大，

MCMC算法的相对误差也小于11%，远低于蒙特卡罗模拟将近40%的相对误差。

对于合并的小样本数据，MCMC算法的相对误差在6.25%以内，证明MCMC算法

同样适用于小样本偏态分布数据，由此证明MCMC算法同样适用于偏态分布数据。

3.在时间趋势上，儿童血铅和成人血铅水平均在无铅汽油政策推行后出现快速

下降。儿童血铅水平维持了较快的下降速度，从1998—2002年的70—80μg/L降至

2012年后的40μg/L以下。而成人血铅水平在“十二五”规划后虽有下降，但速度相

对缓慢。同时，一般人群血镉水平的变化趋势则与“十二五”规划前后土壤镉浓度先

升后降保持一致。

4.在性别差异方面，男性与女性的血铅水平均随着无铅汽油政策和“十二五”规

划的实施呈现下降趋势。值得注意的是，女性的血铅水平下降率始终高于男性，最

高达到63.52%。此外，在不同的重金属污染防治政策时期中，男性的血铅、血镉

水平均高于女性，且女性血镉水平的下降率高于男性。但男性和女性血镉水平受

“十二五”规划的影响均不明显。

5.在年龄方面，各年龄段一般人群的血铅水平同样随着政策的实施出现下降。

学龄前儿童的血铅水平下降率最高，到达64.10%；大于60岁的年龄组血铅下降率

最低，为32.14%。不同年龄组血镉水平变化显示，成年组在2012