济南市复合环境暴露对急性症状的影响:基于LASSO变量筛选与广义加法混合模型评估

Impact of complex environmental exposures on acute symptoms in Jinan: Based on LASSO variable selection and generalized additive mixed models

  • 摘要:
    背景 空气污染与气象因素对人群急性症状有复杂的非线性影响,各项指标间也存在复杂的交互关系,传统统计方法很难同时处理复杂非线性变化和多重共线性干扰问题。
    目的 防止多重共线性对研究结果的干扰,深入探究空气污染物和气象因素在不同浓度水平和环境参数下对3种人群急性症状的动态影响,为相关健康风险因素的防控提供科学依据。
    方法 采用时间序列研究设计,收集2023年6—12月济南市空气污染(每日平均气温、每日降水量、每日平均相对湿度和每日平均风速)、气象因素空气质量指数(AQI)、细颗粒物(PM2.5)、可吸入颗粒物(PM10)、二氧化硫(SO2)、二氧化氮(NO2)、一氧化碳(CO)以及8 h最大臭氧(O3)和发热、咳嗽、咽痛3种急性症状数据。应用最小绝对收缩和选择运算回归(LASSO回归)筛选出关键变量后,引入广义加法混合模型(GAMM)分析空气污染与气象因素复合环境暴露对健康效应的影响,线性关系变量采用线性混合效应函数,非线性关系变量采用薄板回归样条平滑函数,存在交互作用变量采用低秩尺度不变张量积平滑函数。将自变量符合正态分布的波动部分视为抽样误差,作为GAMM随机效应项。
    结果 对于发热,日均气温、日均相对湿度、日均风速和SO2的影响具有统计学意义(均P<0.05),其中日均风速为线性影响因素。日均气温<3 ℃时每升高10 ℃对应RR值为2.64(95%CI:2.50~2.79);≥3 ℃时每升高10 ℃对应RR值为0.86(95%CI:0.83~0.89)。日均相对湿度每升高10%对应RR值为0.93(95%CI:0.89~0.97)。日均风速每增加1 m·s−1对应RR值为1.06(95%CI:1.02~1.10)。SO2在<10 μg·m−3,10~<12.5 μg·m−3,≥12.5 μg·m−3 3段浓度区间内,每升高1 μg·m−3对应RR值分别为1.01(95%CI:0.98~1.05)、1.21(95%CI:1.17~1.24)和0.97(95%CI:0.94~0.99)。对于咳嗽,日均气温、日均相对湿度、PM10和SO2的影响具有统计学意义(均P<0.001),其中PM10为线性影响因素。日均气温<1 ℃时每升高10 ℃对应RR值为1.47(95%CI:1.42~1.52),≥1 ℃时每升高10 ℃对应RR值为0.85(95%CI:0.82~0.87)。日均相对湿度每升高10%对应RR值为0.95(95%CI:0.92~0.98)。PM10每升高50 μg·m−3对应RR值为1.05(95%CI:1.02~1.08)。SO2在<10 μg·m−3、10~<12.5 μg·m−3、≥12.5 μg·m−3 3段浓度区间内,每升高1 μg·m−3对应RR值分别为1.00(95%CI:0.97~1.03)、1.12(95%CI:1.09~1.16)和0.98(95%CI:0.95~1.00)。对于咽痛,日均气温、日均相对湿度、日均风速、PM10和SO2的影响具有统计学意义(均P<0.05),其中日均风速、PM10为线性影响因素且具有交互作用。日均气温<2 ℃每升高10 ℃对应RR值为1.82(95%CI:1.69~1.96);≥2 ℃时每升高10 ℃对应RR值为0.81(95%CI:0.77~0.87)。日均相对湿度每升高10%对应RR值为0.94(95%CI:0.88~1.00)。SO2在<10 μg·m−3、10~<12.5 μg·m−3、≥12.5 μg·m−3 3段浓度区间内,每升高1 μg·m−3对应RR值分别为1.02(95%CI:0.97~1.08)、1.13(95%CI:1.08~1.19)和0.98(95%CI:0.94~1.02)。日均风速和PM10每升高1 m·s−1和50 μg·m−3对应RR值分别为1.06(95%CI:1.00~1.12)和1.04(95%CI:0.98~1.10)。日均风速与PM10对咽痛的影响存在交互效应,日均风速升高时非线性降低PM10影响,PM10对日均风速无影响。
    结论 本研究基于LASSO与GAMM的联合应用,基本消除了指标间的多重共线性影响,揭示了济南市空气污染物和气象因素对不同人群急性症状存在复杂的非线性影响和交互作用。日均气温和SO2浓度与发热、咳嗽和咽痛等症状的发生风险呈非线性关联,而PM10和风速等因素则呈现线性关系或交互效应。这些发现为精准防控相关健康风险因素提供了新的依据。

     

    Abstract:
    Background Air pollution and meteorological factors exert complex nonlinear effects on acute symptoms in the population, with intricate interactions among these factors. Traditional statistical methods struggle to simultaneously address complex nonlinear relationships and multicollinearity issues.
    Objective To delineate the dynamic effects of air pollutants and meteorological parameters on acute symptoms in three distinct populations with the multicollinearity being addressed and to generate reliable scientific evidence for prevention and control of health risk factors.
    Methods A time-series study design was employed to collect data on air pollution (daily mean temperature, daily precipitation, daily mean relative humidity, and daily mean wind speed), meteorological factors Air Quality Index (AQI), fine particulate matter (PM2.5), inhalable particulate matter (PM10), sulfur dioxide (SO2), nitrogen dioxide (NO2), carbon monoxide (CO), and 8-hour maximum ozone (O3), and acute symptoms such as fever, cough, and sore throat in Jinan from June to December 2023. Key variables were selected using least absolute shrinkage and selection operator (LASSO) regression, followed by generalized additive mixed modeling (GAMM) to analyze the health effects of combined environmental exposures to air pollution and meteorological factors. Linear variables were modeled using linear mixed-effects function, nonlinear variables were smoothed using thin-plate regression splines, and variables with interaction effects were smoothed using low-rank scale-invariant tensor product splines. Fluctuations in independent variables following a normal distribution were treated as sampling errors and incorporated as random effects in the GAMM.
    Results For fever, the daily mean temperature, daily mean relative humidity, daily mean wind speed, and ambient SO2 were statistically significant (P<0.05), with daily mean wind speed being a linear influencing factor. When the daily mean temperature was below 3 °C, each 10 °C increase corresponded to a relative risk (RR) of 2.64 (95%CI: 2.50, 2.79). When the daily mean temperature was ≥3 °C, each 10 °C increase corresponded to an RR of 0.86 (95%CI: 0.83, 0.89). Each 10% increase in daily mean relative humidity was associated with an RR of 0.93 (95%CI: 0.89, 0.97). Each 1 m·s−1 increase in daily mean wind speed corresponded to an RR of 1.06 (95%CI: 1.02, 1.10). Within the concentration ranges of <10 μg·m−3, 10–<12.5 μg·m−3, and ≥12.5 μg·m−3, each 1 μg·m−3 increase in ambient SO2 corresponded to RR values of 1.01 (95%CI: 0.98, 1.05), 1.21 (95%CI: 1.17, 1.24), and 0.97 (95%CI: 0.94, 0.99), respectively. For cough, the daily mean temperature, daily mean relative humidity, PM10, and SO2 were statistically significant (P<0.001), with PM10 being a linear influencing factor. When the daily mean temperature was below 1 °C, each 10 °C increase corresponded to an RR of 1.47 (95%CI: 1.42, 1.52). When the daily mean temperature was ≥1 °C, each 10 °C increase corresponded to an RR of 0.85 (95%CI: 0.82, 0.87). Each 10% increase in daily mean relative humidity was associated with an RR of 0.95 (95%CI: 0.92, 0.98). Each 50 μg·m−3 increase in PM10 concentration corresponded to an RR of 1.05 (95%CI: 1.02, 1.08). Within the concentration ranges of <10 μg·m−3, 10–<12.5 μg·m−3, and ≥ 12.5 μg·m−3, each 1 μg·m−3 increase in ambient SO2 corresponded to RR values of 1.00 (95%CI: 0.97, 1.03), 1.12 (95%CI: 1.09, 1.16), and 0.98 (95%CI: 0.95, 1.00), respectively. For sore throat, the daily mean temperature, daily mean relative humidity, daily mean wind speed, PM10, and SO2 were statistically significant (P<0.05), with daily mean wind speed and PM10 being linear influencing factors. When the daily mean temperature was below 2 °C, each 10 °C increase corresponded to an RR of 1.82 (95%CI: 1.69, 1.96). When the daily mean temperature was ≥2 °C, each 10 °C increase corresponded to an RR of 0.81 (95%CI: 0.77, 0.87). Each 10% increase in daily mean relative humidity was associated with an RR of 0.94 (95%CI: 0.88, 1.00). Within the concentration ranges of <10 μg·m−3, 10–<12.5 μg·m−3, and ≥12.5 μg·m−3, each 1 μg·m−3 increase in ambient SO2 corresponded to RR values of 1.02 (95%CI: 0.97, 1.08), 1.13 (95%CI: 1.08, 1.19), and 0.98 (95%CI: 0.94, 1.02), respectively. Each 1 m·s−1 increase in daily mean wind speed and each 50 μg·m−3 increase in PM10 concentration were associated with RR values of 1.06 (95%CI: 1.00, 1.12) and 1.04 (95%CI: 0.98, 1.10), respectively. An interaction effect was observed between daily mean wind speed and PM10: increasing daily mean wind speed non-linearly reduced the impact of PM10, on sore throat whereas PM10 had no significant effect on wind speed.
    Conclusion This study, by combining LASSO and GAMM, largely eliminates the multicollinearity among selected variables. It reveals complex non-linear effects and interactions between air pollutants, meteorological factors, and acute symptoms in different population groups in Jinan. The symptoms like fever, cough, and sore throat are non-linearly associated with daily mean temperature and SO2 concentration, while PM10 and wind speed show a linear relationship or interactive effects. These findings provide a new basis for the precise prevention and control of health risk factors.

     

/

返回文章
返回