坦桑尼亚儿童疫苗接种覆盖率及其相关决定因素的贝叶斯空间分析

《BMJ Public Health》：Bayesian spatial analysis of child vaccination coverage and associated determinants in Tanzania

　　摘要
引言
儿童免疫接种是最有效的公共卫生干预措施之一，显著降低了可预防疾病的发病率和死亡率。尽管全球付出了努力，但包括坦桑尼亚在内的许多低收入和中等收入国家在实现高疫苗接种覆盖率方面仍面临挑战。本研究考察了坦桑尼亚12至23个月大儿童全面疫苗接种覆盖率的空间分布及其决定因素。

方法
本研究采用了2022年坦桑尼亚人口与健康调查-疟疾指标调查的数据，共纳入2143名12至23个月大的儿童，采用横断面设计方法。通过Moran’s I指数和局部空间关联指标分析空间依赖性以检测聚集现象。利用Besag-York-Mollié模型并通过集成嵌套拉普拉斯近似法进行拟合的贝叶斯空间分层逻辑回归模型，评估了个体和社区层面的疫苗接种覆盖率决定因素。

结果
研究发现，仅有54.6%的儿童完成了全部疫苗接种，低于国家目标及世界卫生组织（WHO）推荐的90%的阈值。母亲的教育水平、产前护理（ANC）参与情况以及家庭财富是影响全面疫苗接种的重要因素。母亲受教育程度较高的儿童完成疫苗接种的概率更高（OR=3.78；95% CI 1.38至12.2），而母亲接受四次或以上产前护理的儿童完成疫苗接种的概率也更高（OR=1.24；95% CI 0.97至1.59）。来自中等收入家庭的儿童完成疫苗接种的概率也较高（OR=1.39；95% CI 1.04至1.85）。母亲年龄和分娩地点在统计上没有显著影响。空间分析显示疫苗接种存在显著地理聚集现象（全球Moran’s I=0.229，p=0.028），其中Kusini Unguja、Kigoma、Mara和Njombe地区的覆盖率较高，而Songwe、Rukwa和Shinyanga地区的覆盖率较低。

结论
本研究强调了提高母亲教育水平、加强产前护理咨询、减少经济和获取障碍以及实施针对不同地区的干预措施的重要性，以改善坦桑尼亚的儿童疫苗接种覆盖率并减少不平等现象。

已有研究
疫苗接种每年挽救了数百万人的生命，但在撒哈拉以南非洲地区，疫苗接种覆盖率仍不均衡。2022年坦桑尼亚仅有53%的儿童完成了全部疫苗接种，这种差异与贫困、地理因素和母亲相关因素有关。以往的研究很少考虑空间变异。

本研究的新发现
通过使用Besag-York-Mollié模型的贝叶斯空间分层逻辑回归，本研究确定了母亲教育水平、产前护理和家庭财富是影响儿童疫苗接种的关键因素，并揭示了各地区之间的明显地理聚集现象。

研究对研究、实践或政策的影响
研究结果为政策制定者提供了证据，以设计有针对性的地区性干预措施并加强母亲教育，从而提高疫苗接种率。

引言
疫苗接种被广泛认为是最有效的公共卫生干预措施之一，能有效降低可预防疾病的发病率和死亡率。疫苗通过刺激免疫系统产生保护性抗体，预防麻疹、脊髓灰质炎、肺炎和破伤风等严重感染。仅儿童免疫接种每年就能防止超过四百万例死亡，并预计在2021至2030年间挽救约520万人的生命。尽管如此，疫苗接种覆盖率仍存在差异，使边缘化和脆弱人群面临更高的疾病风险。全球每年有近2000万婴儿未能接种必要的疫苗，凸显了免疫工作的持续缺口。为应对这些挑战，国际社会发起了多项全球性倡议，其中世界卫生组织于1974年推出的扩大免疫计划为全球常规疫苗接种计划奠定了基础。此外，全球疫苗和免疫联盟等组织以及WHO-UNICEF等合作伙伴关系在资助和支持低收入和中等收入国家的疫苗分发方面发挥了重要作用。最近，《免疫议程2030》设定了雄心勃勃的目标，旨在确保每个人都能在一生中受益于全面的疫苗接种。这些努力显著扩大了疫苗的可及性，降低了许多国家的儿童死亡率。

然而，近年来这些全球性努力面临新的挑战。特别是COVID-19大流行扰乱了全球医疗服务，加剧了免疫差距。截至2023年，从未接种过任何疫苗的儿童人数增加了40%，达到1450万。目前全球三剂白喉、破伤风和百日咳（DTP3）疫苗的接种率为84%，第一剂麻疹疫苗的接种率为83%，均低于大流行前的水平。虽然像人乳头瘤病毒疫苗这样的新疫苗在某些地区的接种率有所提高，但贫困国家由于医疗基础设施薄弱和物流障碍，疫苗分发仍存在持续差异。这些全球趋势在撒哈拉以南非洲尤为明显，该地区5岁以下儿童死亡率居高不下，尽管仅占全球出生人口的30%。坦桑尼亚的疫苗接种覆盖率仍低于全球平均水平，DTP3疫苗的接种率为76%，第一剂麻疹疫苗的接种率为67%。由于医疗系统薄弱、地理隔离以及COVID-19大流行的影响，该地区仍有约800万儿童未接种疫苗。尼日利亚和埃塞俄比亚两国未接种疫苗的儿童数量超过300万，凸显了提供基本免疫服务的系统性挑战。

在坦桑尼亚，该国在实现全民疫苗接种方面取得了进展，但也存在持续挑战。该国为2岁以下儿童制定了全面的常规免疫计划，包括出生时或首次就医时接种卡介苗（BCG）；出生后6周、10周和14周接种口服脊髓灰质炎疫苗（OPV）；14周接种灭活脊髓灰质炎疫苗；以及6周、10周和14周接种轮状病毒疫苗、肺炎球菌结合疫苗（PCV13）和五联疫苗。此外，9个月和18个月还接种麻疹-风疹（MR）疫苗。尽管坦桑尼亚全国范围内都提供儿童免疫服务，但最新数据显示疫苗接种覆盖率有所下降。12至23个月大的儿童中，仅有53%完成了全部基础疫苗接种：91%接种了BCG疫苗，90%接种了三剂DPT疫苗，59%接种了三剂OPV疫苗，87%接种了MR疫苗。更令人担忧的是，只有49%的儿童在1岁前完成了所有基础疫苗接种。虽然1991-1992年和2015-2016年间接种率相对稳定（分别为71%和75%），但2022年骤降至53%。同一时期，未接种疫苗的儿童比例在2%至5%之间波动。鉴于坦桑尼亚政府为扩大接种范围所做的努力（如免除疫苗费用和将免疫接种纳入其他医疗服务），这一下降趋势尤为明显。然而，低收入家庭和农村地区的儿童仍面临医疗设施有限、交通不便和经济困难等障碍。游牧和偏远地区的儿童由于医疗配送的物流问题，疫苗接种率也较低。许多研究都指出了全球和地区层面的类似不平等现象。影响儿童疫苗接种率的因素已被广泛记录，包括社会经济地位、地理可及性和母亲特征。然而，大多数分析依赖于传统的统计方法（如逻辑回归或泊松回归），这些方法无法充分捕捉空间变异。此外，坦桑尼亚的大部分现有研究集中在特定地区、年龄组或种族群体上。最近的空间分析技术进步为检测地理聚集现象、绘制疫苗差距图谱及其决定因素提供了新机会。尽管有少数研究探讨了儿童疫苗接种的空间变异，但没有研究应用灵活的结构化加法回归模型同时估计地区空间效应和控制个体及家庭层面的协变量。本研究利用2022年坦桑尼亚人口与健康调查数据，通过贝叶斯空间分层逻辑回归模型（BYM）探讨了12至23个月大儿童全面疫苗接种的决定因素，考虑了社会人口统计、经济和健康相关因素。研究结果有望为提高疫苗接种覆盖率、减少不平等以及实现可持续发展目标3（确保所有年龄段的健康生活）提供指导。

方法
本研究使用了2022年坦桑尼亚人口与健康调查（TDHS）和疟疾指标调查（TDHS-MIS）的横断面数据。TDHS-MIS自20世纪90年代以来定期开展，提供了全国范围内的健康指标数据。2022年的调查收集了关于母亲和儿童健康、营养以及儿童疫苗接种覆盖率的详细信息。

样本设计和研究对象
本研究聚焦12至23个月大的儿童，因为这一年龄段是评估全面疫苗接种率的关键时期。抽样方法遵循2022年TDHS-MIS的设计，采用分层两阶段抽样。第一阶段选取了629个样本群，其中城市地区211个，农村地区418个；第二阶段从每个样本群中系统选取26户家庭，共计16,354户家庭。经过纳入标准和数据清洗后，最终分析样本包括2143名12至23个月的儿童。

研究变量
主要结果变量是全面疫苗接种覆盖率，这是一个二元变量，表示儿童在12至23个月大时是否完成了所有基础疫苗接种。根据WHO建议，全面基础疫苗接种包括一剂BCG疫苗、三剂DPT-HepB-Hib疫苗、三剂脊髓灰质炎疫苗和一剂麻疹疫苗。

解释变量
解释变量来自先前文献，分为三类：
1. 社会人口统计变量：儿童性别、居住地、出生顺序、母亲年龄、母亲教育水平和家庭户主性别。
2. 经济变量：家庭财富指数、父母就业状况和接触大众媒体的情况。
3. 健康相关变量：产前护理次数、分娩地点和产后护理次数。这些变量的测量结果见在线补充表1。

统计分析
数据使用R软件（V.4.5.0）进行处理和分析。分析样本包括2143名调查时仍在世的12至23个月大的儿童。所有儿童均具有完整的BCG、五联疫苗、脊髓灰质炎疫苗、轮状病毒疫苗和麻疹疫苗接种记录，以及所有分析中考虑的协变量（母亲年龄、母亲教育水平、家庭财富、生育次数、居住地和地区）的完整信息。因此，没有儿童因数据缺失而被排除，从而减少了选择偏差的风险。2022年TDHS-MIS数据集中没有提供具体的接种月份或季度信息，因此无法评估年度内的时间变化。地理位置根据儿童调查时的居住地区确定，而非接种服务的地点。

为了考虑复杂的调查设计，引入了抽样权重、样本群和分层。生成了加权描述性统计数据和可视化结果，随后进行单变量逻辑回归以评估初步关联。单变量分析中p≤0.20的变量或先前文献支持的变量被纳入多变量模型。重要的是，所有具有理论或实证支持的协变量（包括母亲教育水平、家庭财富和产前护理次数）都被纳入多变量贝叶斯空间分层模型，无论其单变量p值如何。这种策略考虑了相关预测因子，减少了模型估计的潜在偏差或不稳定。参考类别代表了最弱势或基线群体（例如，无母亲教育、最贫困的财富五分位数、无产前护理），为解释OR值提供了有意义的基准。这种方法优先考虑政策相关性和可解释性，允许在不同社会人口统计和健康因素之间进行清晰比较。

全球空间自相关性通过Moran’s I指数进行评估，计算公式为：
\[Moran’s I = \frac{\sum_{i=1}^{N} \sum_{j=1}^{N} \frac{d_{ij}}{N(n-1)}\]
其中 \(d_{ij}\) 表示地区i的疫苗接种率， \(\bar{d}\) 表示全国平均疫苗接种率， \(w_{ij}\) 表示地区i和j之间的空间权重，N表示地区总数。正的且统计显著的Moran’s I指数表明疫苗接种率存在空间聚集现象。随后采用了局部空间关联指标（LISA）来识别具有显著局部聚集性的区域，将区域分类为高-高（疫苗接种覆盖率高的区域被高覆盖率邻居包围）、低-低（覆盖率低的区域被低覆盖率邻居包围）、高-低（高异常值区域被低覆盖率区域包围）或低-高（低异常值区域被高覆盖率区域包围），并通过999次蒙特卡洛排列来评估统计显著性。为了实现全面疫苗接种覆盖率，使用了基于BYM规范的贝叶斯空间层次逻辑回归模型来模拟区域特定效应。该模型表示为：其中表示儿童i的全面疫苗接种概率，是全局截距，表示儿童i的协变量值，是具有高斯先验的回归系数，表示模糊的先验信息以确保对后验估计的影响最小。是遵循条件自回归分布的空间结构化随机效应，而则是捕捉区域特定异质性的非结构化随机效应。为这两个组成部分指定了惩罚复杂性（PC）先验以稳定估计。后验分布使用积分嵌套拉普拉斯近似（INLA）进行近似，并报告了后验均值和95%的可信区间。后验空间效应被映射为，以可视化在调整协变量后疫苗接种覆盖率高于或低于预期的区域。

研究参与者的描述性特征和疫苗接种覆盖率
表1展示了2143名12-23个月大儿童的特征，其中54.6%的儿童已完全接种疫苗。城市地区的疫苗接种覆盖率略高于农村地区（分别为56.4%和53.9%），并且在不同地区之间存在差异，从桑给巴尔的近70%到南部地区的42%不等。按出生顺序来看，57.3%的第一个孩子已完全接种疫苗，而后来出生的孩子则较低。覆盖率还受到母亲教育水平的影晌，从没有受过教育的母亲的孩子的41.8%到母亲受教育程度较高的孩子的75%不等，以及母亲年龄的影响，从15-19岁母亲的50%到35-49岁母亲的56.4%不等。来自较富裕家庭的儿童（高达60%）以及母亲有媒体接触的儿童（57.6%）的疫苗接种覆盖率较高。在健康相关因素中，56.3%在医疗机构出生的儿童和57.7%的母亲进行了四次或更多次产前检查的儿童已完全接种疫苗。男性的疫苗接种覆盖率略高于女性（分别为55.5%和53.8%），男性主导的家庭的覆盖率也略高于女性主导的家庭（分别为54.9%和53.6%）。

表1 • 按关键因素划分的儿童疫苗接种覆盖率分布

坦桑尼亚各地区儿童疫苗接种的空间分布和聚集
计算了全球Moran’s I统计量，以评估12-23个月大儿童全面疫苗接种覆盖率的空间自相关性。分析得出的Moran’s I值为0.229，零假设下的期望值为-0.033，方差为0.019，标准差为1.906，p值为0.028，表明在5%的水平上具有统计显著性（见在线补充表2）。这些结果表明存在中等程度的正空间自相关性，意味着疫苗接种覆盖率相似的区域倾向于在地理上聚集而不是随机分布，反映了区域对疫苗接种的影响。为了进一步说明这种空间模式，Moran散点图提供了这种聚集的视觉表示。图1中的Moran散点图显示了区域全面疫苗接种覆盖率（x轴）与空间滞后值（y轴）之间的关系，后者代表邻近区域的平均覆盖率。向上倾斜的回归线证实了这种正空间关系，与Moran’s I的结果一致。右上象限和左下象限的区域表明存在正空间自相关性，即高覆盖率区域被高覆盖率邻居包围，低覆盖率区域与低覆盖率区域相邻。一个显著的异常值Mjini Magharibi位于右上象限，其疫苗接种覆盖率相对于邻近区域异常高。尽管有这种影响，但Moran’s I（I=0.229，p=0.028）所表明的整体正空间自相关性在大多数区域仍然稳健，显示出全国范围内相似疫苗接种覆盖率的统计显著聚集。

基于这些全局模式，基于局部Moran’s I统计量的LISA聚类图（图2）揭示了全面疫苗接种覆盖率的显著局部空间自相关性。在Kusini Unguja和Mjini Magharibi Pemba检测到了具有统计意义的高-高聚类，表明这些区域及其邻居的覆盖率一直很高。相反，在Songwe检测到了低-低聚类，显示出低疫苗接种覆盖率的集中。Kigoma作为一个高-低异常值出现，其特征是相对较高的覆盖率被覆盖率较低的区域所包围。未被归类到这些类别的区域被认为没有显著的空间自相关性，反映了没有检测到局部空间自相关性。这些发现突出了某些特定地理区域，其疫苗接种覆盖率明显高于或低于周围区域，表明有针对性、基于空间信息的公共卫生干预措施的机会。

与全面疫苗接种覆盖率相关的因素
首先进行了单变量逻辑回归分析，以检查每个解释变量与儿童全面疫苗接种覆盖率之间的粗略关联。在单变量分析中显示出显著关联的变量，以及先前文献支持的变量，随后被纳入贝叶斯层次模型中。表2中的结果表明，母亲的教育水平与疫苗接种有很强的关联：母亲具有小学、中学或更高教育水平的儿童完全接种疫苗的概率分别比母亲没有受过教育的儿童高1.54倍（95% CI 1.22至1.96）、1.87倍（95% CI 1.38至2.54）和3.78倍（95% CI 1.38至12.2）。高等教育类别仅包括20名儿童，这在贝叶斯层次模型估计中产生了较宽的可信区间，反映了这一组效应大小的较大不确定性。母亲年龄、分娩地点和大众媒体接触没有显示出统计显著的关联。进行四次或更多次产前检查的儿童完全接种疫苗的概率更高（OR=1.24，95% CI 0.97至1.59），尽管在5%的水平上这种关联并不显著。财富指数显示，来自较富裕和中等收入家庭的儿童完全接种疫苗的概率高于最贫困家庭的儿童，分别为1.33倍（95% CI 1.01至1.75）和1.39倍（95% CI 1.04至1.85）。居住区域也影响了覆盖率：南部和西南高地的儿童完全接种疫苗的概率较低（分别为OR=0.53，95% CI 0.31至0.91；OR=0.66，95% CI 0.45至0.96），而桑给巴尔的概率较高（OR=1.49，95% CI 1.00至2.22）。

表2 • 与儿童疫苗接种覆盖率相关的因素的贝叶斯层次模型结果

贝叶斯层次空间模型的结果显示，坦桑尼亚12-23个月大儿童全面疫苗接种覆盖率存在显著的地理差异。如图3所示，后验空间效应突出了在调整社会人口统计、经济和卫生服务相关协变量后仍然存在的区域差异。

图3 • 坦桑尼亚12-23个月大儿童全面疫苗接种的后验空间效应。正值（红色）表示在调整模型协变量后覆盖率高于平均水平的区域，而负值（蓝色）表示覆盖率低于预期的区域。例如Kilimanjaro、Mara、Njombe、Singida和Kigoma等区域显示出正空间效应。这意味着生活在这些区域的儿童比仅根据其个人和家庭特征预期的接种覆盖率更高。这些正面效应反映了有利的当地条件，如更强的免疫计划实施、更好的卫生服务运作或更高的社区意识和疫苗接种支持。相比之下，Rukwa和Shinyanga等区域显示出负空间效应，表明在调整观察到的因素后接种覆盖率低于预期。这表明这些地区存在额外的挑战，如难以获得卫生设施、交通障碍、卫生系统限制或阻碍疫苗接种的社会文化因素。

Kigoma的空间效应与早期的LISA聚类分析结果一致，该地区表现为高-低异常值，意味着其疫苗接种覆盖率相对于邻近区域较高。贝叶斯空间模型通过显示Kigoma的正空间效应证实了这一模式，表明该地区受益于支持疫苗接种的有利环境或计划相关因素，超出了个人和家庭特征的影响。这些发现与贝叶斯逻辑回归模型的结果一致，后者显示来自较富裕家庭的儿童和母亲受教育程度较高的儿童更有可能完全接种疫苗。然而，在考虑了这些因素后，区域差异的持续存在强调了地理位置本身在塑造疫苗接种结果中的重要作用。

总体而言，结果表明，提高坦桑尼亚的疫苗接种覆盖率需要采取双管齐下的方法：解决个人层面和家庭层面的因素，同时实施针对当地条件的区域特定策略。因此，考虑区域背景的针对性干预措施对于减少全国范围内的持续疫苗接种不平等至关重要。

讨论
本研究使用了2022年TDHS-MIS的数据，通过应用基于BYM规范的贝叶斯空间层次逻辑回归并使用INLA方法，研究了坦桑尼亚12-23个月大儿童全面疫苗接种覆盖率的空间分布和决定因素。通过整合空间效应和后验估计，结合全球Moran’s I和LISA等空间自相关技术，研究识别出显著的区域差异以及与儿童免疫相关的关键社会人口统计、经济和卫生服务相关因素。描述性发现显示，只有54.6%的12-23个月大儿童已完全接种疫苗，远低于坦桑尼亚的国家免疫目标和世界卫生组织推荐的90%的覆盖率阈值。这一表现不佳的情况凸显了加强国家免疫计划和改进外展策略的迫切需要。低覆盖率反映了持续的障碍，如薄弱的卫生基础设施、有限的访问机会和低水平的母亲意识，这些都可能阻碍及时和全面的疫苗接种。其他东非国家（如埃塞俄比亚和肯尼亚）的类似趋势表明，区域合作和知识共享可以为有效的干预措施提供信息。

母亲的教育水平是儿童疫苗接种的强预测因子。母亲具有中学或更高教育水平的儿童完全接种疫苗的概率显著增加。这意味着提高女性教育和健康素养可能对儿童免疫率产生长期影响。母亲的健康素养增强了寻求健康的行为，增强了对免疫信息的理解，并提高了对疫苗接种计划的遵守度，这在尼日利亚和更广泛的撒哈拉以南非洲地区都有观察到。因此，健康促进计划应整合针对母亲的教育干预措施，特别是在覆盖率低的地区。产前检查与全面疫苗接种呈正相关，母亲进行四次或更多次产前检查的更有可能为孩子接种疫苗。这一发现意味着加强产前服务并确保在怀孕期间提供全面的儿童免疫咨询可以改善疫苗接种情况。卫生系统应强调产前检查作为促进儿童健康干预的平台，这与马拉维和埃塞俄比亚的发现一致。家庭财富对疫苗接种覆盖率有适度但积极的影响。来自中等或较富裕家庭的儿童更有可能完全接种疫苗，这表明财务资源可以减少物流障碍，包括交通费用和服务费用。针对较贫困家庭的计划、外展活动或移动疫苗接种服务可以帮助缓解不平等，正如在布基纳法索和尼日利亚所注意到的。然而，最富裕五分位的停滞表明，私人医疗保健的使用（未在国家记录中捕捉到）可能会影响观察到的覆盖率模式。

有趣的是，在调整其他变量后，母亲年龄和分娩地点并未成为统计显著的预测因子。这表明，在坦桑尼亚的背景下，当考虑到更广泛的服务质量和社区层面的决定因素时，这些因素的影响较小。因此，政策应侧重于系统性改进，而不仅仅是针对母亲年龄或分娩地点。空间分析显示了疫苗接种覆盖率的显著区域聚集。全球Moran’s I值为0.229（p=0.028），证实了空间自相关性，表明高和低覆盖率区域在地理上趋于集中。这表明干预措施应具有地理针对性，优先考虑低覆盖率区域以加强外展和资源分配。在莫桑比克和埃塞俄比亚也观察到了类似的模式，那里不均匀的卫生基础设施和项目覆盖范围导致了差异。

LISA分析在Kusini Unguja和Mjini Magharibi Pemba识别出了高-高聚类，在Songwe识别出了低-低聚类，分别指出了项目成功的区域和迫切需要的区域。研究结果非常明确：疫苗接种覆盖率高的地区可以作为最佳实践的典范，而覆盖率低的地区则需要制定针对性的策略，例如加强当地卫生系统的能力、提高社区参与度以及开展疫苗接种活动。贝叶斯空间模型进一步指出，基戈马（Kigoma）、马拉（Mara）和恩乔姆贝（Njombe）地区的空间效应高于平均水平，这表明由于成功的疫苗接种活动和非政府组织（NGO）的支持，这些地区存在积极的、未被充分测量的环境因素影响。相反，鲁夸（Rukwa）和希尼亚加（Shinyanga）地区的空间效应显著为负，这表明这些地区存在持续存在的结构性障碍，如地形恶劣、民众对疫苗接种的抵触情绪以及服务提供有限等问题。这些结果强调了基于地理信息制定政策决策和进行有针对性的资源分配的必要性，以解决不平等问题，这与尼日利亚的研究发现一致——在该国，边缘化地区也面临着类似的健康障碍。值得注意的是，虽然LISA分析根据原始疫苗接种覆盖率将基戈马列为高-低异常值，但贝叶斯空间模型在调整了个体和家庭层面的因素后，仍然显示该地区存在积极的空间效应。这种一致性表明，基戈马相对较高的疫苗接种覆盖率并非仅仅由有利的社会人口统计特征所解释，还得到了积极的区域环境因素的支持，例如有效的疫苗接种服务提供、支持性的当地卫生系统结构或强烈的社区参与度。描述性LISA分析与调整后的贝叶斯估计结果之间的一致性进一步证明了地理环境在该地区疫苗接种结果中起着重要作用。

**研究的优势与局限性**
当前研究具有几个优势。它采用了先进的贝叶斯空间层次建模方法，并考虑了空间依赖性以揭示地理聚集现象；同时使用了具有全国代表性的2022年TDHS-MIS数据，确保了研究的普遍适用性。然而，也存在一些局限性：研究未能涵盖文化因素，并且分析的是整体疫苗接种情况而非特定疫苗的接种情况；由于缺乏关于疫苗接种或调查实施月份或季度的详细信息，无法评估2022年内的时间异质性，这可能会掩盖与COVID-19服务中断和季节性趋势相关的疫苗接种覆盖率变化；地理位置是根据调查时儿童的居住地区来定义的，而非实际接种疫苗的地点，如果家庭在不同地区之间迁移，可能会导致分类错误；受教育程度较高的母亲所生的儿童数量较少，可能导致这一类别的估计结果不够精确且置信区间较宽；未测量因素（如文化习俗或当地服务提供差异）可能产生的残余混杂效应可能会影响观察到的关联结果，从而限制了因果关系的解释。此外，使用二手TDHS-MIS数据限制了通过原始数据收集来测量其他可能相关因素的能力，例如更详细的行为、文化或服务提供决定因素。

**结论**
本研究利用2022年TDHS-MIS数据，探讨了坦桑尼亚12至23个月大儿童的整体疫苗接种覆盖率的空间分布及其决定因素。研究发现，母亲的教育水平、产前护理（ANC）参与情况以及家庭财富是影响疫苗接种覆盖率的重要因素，而母亲年龄和分娩地点在调整后并未显示出统计学上的显著影响。空间分析显示，库西尼乌恩古贾（Kusini Unguja）、基戈马（Kigoma）、马拉（Mara）和恩乔姆贝（Njombe）地区的疫苗接种覆盖率较高，而松韦（Songwe）、鲁夸（Rukwa）和希尼亚加（Shinyanga）地区的覆盖率较低。研究结果表明，受教育程度较高的母亲所生的儿童、来自较富裕家庭的儿童以及母亲参加了四次或更多次产前护理的儿童，其完全接种疫苗的概率更高。区域和地方因素也起到了重要作用，这一点通过Global Moran’s I指数和LISA分析的聚类模式得到了证实。基于这些结果，政策建议包括提高母亲的教育水平和健康素养、加强产前护理服务并提供疫苗接种咨询、解决贫困家庭的财务和获取障碍问题，以及针对低覆盖率地区实施特定干预措施。持续监测和制定针对具体情况的策略对于提高坦桑尼亚儿童的疫苗接种覆盖率并减少不平等现象至关重要。