大气污染的健康危害是全球关注的公共卫生问题。大气污染与人群健康效应研究证据是环境空气质量标准与大气污染治理方案制定的重要依据。近年来，因果推断方法逐步引入到环境流行病学观察性研究中，为大气污染与人群健康效应的因果关系研究提供了更多的统计方法选择。观察性研究的复杂混杂偏倚问题是因果推断的重大挑战。本研究主要介绍了目前应用于大气污染与人群健康效应研究的混杂因素识别与控制的各因果推断方法，旨在为我国大气污染与人群健康效应因果推断研究提供方法学上的参考和依据。

目的:运用双样本双向孟德尔随机化法，探讨人免疫细胞与增生性瘢痕（HS）之间的因果关系。方法:该研究基于双样本双向孟德尔随机化（MR）法，分别从全基因组关联分析catalog数据库和FinnGen数据库获取731种人免疫细胞与HS的数据集。设置显著阈值，筛选与免疫细胞及HS显著相关的单核苷酸多态性（SNP）并排除弱工具变量偏倚的影响。采用逆方差加权（IVW）法［同时采用错误发现率（FDR）中的Benjamini-Hochberg法校正其 P值］初步检测免疫细胞与HS的因果关系并筛选与HS有显著因果关系的免疫细胞。进一步地，采用双样本MR分析的5种方法：IVW法、加权中位数法、简单模式法、加权模式法和MR-Egger法绘制散布图检测筛选出的免疫细胞与HS的因果关系。针对符合假设的免疫细胞SNP，进行Cochran Q检验评估异质性，进行MR-Egger回归和MR-PRESSO法排除水平多效性，进行留一法分析来评估显著结果是否由单个SNP决定。采用双样本MR分析中IVW法反向检测HS与免疫细胞之间的因果关系。 结果:731种免疫细胞达到显著阈值的SNP数量从7个到1 786个不等，HS中达到显著阈值的SNP数量有119个，这些SNP的 F值均>10，提示具有弱工具变量偏倚的可能性较小。IVW法分析显示，60种免疫细胞与HS具有潜在因果关系（ P值均<0.05），经Benjamini-Hochberg法校正后，仅CD45RA和CD39双阳性调节性T细胞（Treg）与HS具有潜在因果关系（ PFDR<0.05）。IVW法（比值比为1.16，95%置信区间为1.08~1.24， P<0.05、 PFDR<0.05）、加权中位数法（比值比为1.16，95%置信区间为1.05~1.28， P<0.05）、加权模式法（比值比为1.14，95%置信区间为1.02~1.27， P<0.05）和MR-Egger 法（比值比为1.18，95%置信区间为1.07~1.30， P<0.05）的散布图均显示，CD45RA和CD39双阳性Treg的14个SNP与HS患病风险存在因果关系；仅简单模式法的散布图显示，CD45RA和CD39双阳性Treg的14个SNP与HS患病风险的因果关系不明显（ P>0.05）。Cochran Q检验显示，CD45RA和CD39双阳性Treg与HS的因果关系不存在异质性（ P>0.05）。MR-Egger回归、MR-PRESSO法分析显示，CD45RA和CD39双阳性Treg与HS的显著因果关系不存在水平多效性（ P>0.05）。留一法分析显示，CD45RA和CD39双阳性Treg与HS的显著因果关系在逐个剔除SNP后结果稳定。双样本反向MR分析显示，HS与731种免疫细胞均无潜在因果关系（ P>0.05）。 结论:从遗传学的角度揭示，免疫细胞CD45RA和CD39双阳性Treg可能会增加HS患病风险。

观察和实验及其设计、相关的概念和内涵仍然模糊，影响着研究设计的正确认识及因果推断有效性的判断。本文借用物理学的相变概念，结合因果思维与因果图，首先建立了事物的属性及其状态、事件与现象之间的关系；接着识别了两种现象获取方式——人为观察和人为实验之间相反的因果结构，及其与暴露的干预、分配在因果推断研究中受各自的因果机制影响；最后，识别了干预系因果概念，是自然因果律基础上联系已知的现象和/或需经测量而获得的现象之间的核心纽带；基于此分析了研究设计二分类的策略，以干预法和非干预法的分类更为清晰。观察和实验提供了所有科学知识的基础，应是内涵统一的概念。基于因果律及其测量过程可能是实现研究设计精准分类的最佳选项之一，值得深入研究。

目的:分析老年原发性高血压患者的血清25羟基维生素D［25-（OH）-D］水平，探讨血清25-（OH）-D水平与老年原发性高血压患者靶器官损伤的相关性。方法:选取浙江省荣军医院2019年1月至2021年4月收治的老年原发性高血压患者163例为研究对象，收集患者一般资料及实验室指标。根据靶器官损伤数目、患者25-（OH）-D水平、单一靶器官损伤部位分别进行分组。采用Pearson相关分析进行两个变量间的相关性分析。采用多因素logistic回归分析进行靶器官损伤影响因素分析。结果:不同数目靶器官损伤患者年龄、病程、25-（OH）-D水平、收缩压以及舒张压比较均差异有统计学意义（ F=16.95、14.39、14.95、33.18、20.88，均 P < 0.001）；不同25-（OH）-D水平患者三酰甘油（TG）以及25-（OH）-D水平比较均差异有统计学意义（ F=2.31、178.48，均 P < 0.05），且25-（OH）-D水平与TG水平呈负相关。心脏、肾脏、颈动脉单一靶器官损伤患者组25-（OH）-D水平与无靶器官损伤组比较显著降低（ t=9.24、6.15、6.09，均 P < 0.05），但单一靶器官损伤患者组间比较均差异无统计学意义（均 P > 0.05）。年龄偏大、病程长、收缩压和舒张压高以及25-（OH）-D水平低是原发性高血压患者靶器官损伤的独立危险因素（ P < 0.05）。 结论:血清25-（OH）-D水平与老年原发性高血压患者靶器官损伤具有一定程度的相关性。但维生素D缺乏与原发性高血压靶器官损伤的证据尚且不足，需要进一步深入研究。本研究具备显著创新性和科学性。

近年来，孟德尔随机化的应用在全球范围内越来越广泛。它以基因作为工具变量，可以排除体系中混杂因素和反向因果的干扰，有效地判断暴露因素与结局的因果关联，对于医学领域的病因探索具有独特优势。本文对国内外孟德尔随机化在眼科和其他医学领域的应用进行综述，为相关性研究和因果研究在眼科的发展提供新思路。 （中华眼科杂志，2021，57：791-795）

目的:探究人肠道菌群与增生性瘢痕（HS）之间的因果关系。方法:该研究为基于双样本孟德尔随机化（TSMR）分析的研究。从全基因组关联分析数据库获得人肠道菌群（18 473个样本）和HS（208 248个样本）的数据。提取已知肠道菌群门、纲、目、科和属5个水平的遗传变异基因，即单核苷酸多态性（SNP）作为工具变量，并进行连锁不平衡（LD）分析。使用PhenoScanner V2数据库进行人类基因型-表型关联分析，排除与HS不相关的肠道菌群SNP并分析获得的SNP是否为弱工具变量。采用TSMR分析的4种方法，即逆方差加权（IVW）、MR-Egger回归、加权中位数和加权模式，对肠道菌群SNP与HS之间的因果关系进行分析。绘制前述4种分析方法得出重要结果的散布图，分析肠道菌群SNP与HS的相关性。采用IVW检验和MR-Egger回归检验评估肠道菌群SNP的异质性，采用MR-Egger回归检验和MR-PRESSO离群值检验评估肠道菌群SNP的水平多效性，采用留一敏感性分析确定HS是否由肠道菌群中的单个SNP引起。对HS SNP与肠杆菌属或HS SNP与瘤胃球菌2属分别进行逆向TSMR分析，检测它们之间是否存在反向因果关系。结果:共获得196个已知肠道菌群，属于9门、16纲、20目、32科、119属，从每个菌群中均获得多个SNP作为工具变量。LD分析显示，除rs1000888、rs12566247、rs994794外，其余肠道菌群SNP均符合遗传变异与暴露因素密切相关的假设。人类基因型-表型关联分析显示，经LD分析后获得的SNP均未被排除且均不是弱工具变量。IVW、MR-Egger回归、加权中位数和加权模式的TSMR分析显示，肠杆菌属及瘤胃球菌2属均与HS存在因果关系。其中，IVW和MR-Egger回归的森林图分析还显示，肠杆菌属的16个SNP（该菌属SNP个数下同）和瘤胃球菌2属的15个SNP（该菌属SNP个数下同）均为HS的保护因素。进一步地，IVW分析显示肠杆菌属SNP（比值比为0.62，95%置信区间为0.41~0.93， P<0.05）和瘤胃球菌2属SNP（比值比为0.62，95%置信区间为0.40~0.97， P<0.05）均与HS发生风险呈负相关。散布图显示，肠杆菌属和瘤胃球菌2属的SNP均是HS的保护因素。IVW检验和MR-Egger回归检验均显示，肠杆菌属SNP（ Q值分别为5.73、5.76， P>0.05）和瘤胃球菌2属SNP（ Q值分别为13.67、15.61， P>0.05）均无异质性。MR-Egger回归检验显示，肠杆菌属和瘤胃球菌2属的SNP均无水平多效性（截距分别为0.01、0.06， P>0.05）；MR-PRESSO离群值检验显示，肠杆菌属和瘤胃球菌2属的SNP均无水平多效性（ P>0.05）。留一敏感性分析显示，无单个肠道菌群SNP驱动HS的发生。逆向TSMR分析显示，HS SNP与肠杆菌属之间、HS SNP与瘤胃球菌2属之间均不存在任何反向因果关系（比值比分别为1.01、0.99，95%置信区间分别为0.97~1.06、0.96~1.04， P>0.05）。 结论:人肠道菌群与HS之间存在一定因果关系，肠杆菌属及瘤胃球菌2属对抑制HS有一定作用。

断点回归设计是一种严格的类实验方法，基于样本在临界值周围分布的随机性，使用两侧数据分别进行回归，从而避免潜在混杂因素的干扰，为因果关联的推断提供强有力的证据。本研究对断点回归方法以及其研究进展进行了介绍，并结合案例对其应用方法进行了详细阐述，以期为断点回归在流行病学研究中的使用提供参考。

基于变量调查（或测量）的共时性、统计学关联及幸存者偏倚等原因，横断面研究被认为不能进行因果推断，这是当前的"共识"。本文基于因果思维，借助因果图，首先明确定义真实截面和测量截面，并识别截面概念仅存在于理论的特性。实际横断面研究中，测量变量的共时性并不存在，而是无一例外地表现为非共时性时序，其实质上相当于测量变量间互为独立性假设，或不存在有差别错分偏倚。类似于累积病例对照研究和历史性队列研究，横断面研究均为暴露和结局已存在或发生后进行的测量，这种测量相当于对变量值的历史重建或"考古"。这类研究进行因果推断的共性前提条件之一是，测量变量与其历史变量间必须存在着因果律。测量变量均为真实变量的替代者，测量变量间的时序在因果推断上并不重要。应加强对横断面研究分析性角色的认识。

非酒精性脂肪性肝病(NAFLD)是以肝内脂质大量堆积为特点的肝脂质代谢异常类疾病,是目前世界范围内最常见的肝脏疾病.孟德尔随机化(MR)将基因组数据纳入传统的流行病学研究设计中,以推断暴露因素与疾病风险之间的因果关系.近年来,MR在NAFLD的病因推断研究中得到了较为广泛的应用,本文将系统总结MR在NAFLD研究中的应用进展,为认识疾病本质和科学化干预提供新思路.

目的:研究卒中后非痴呆型认知障碍(PS-CIND)患者的认知损害与同型半胱氨酸(HCY),脑白质疏松(LA)及P300检测的相关性.方法:对200例PS-CIND患者进行MoCA 评分,HCY检测、LA评分,P300潜伏期、波幅测定,进一步分析MoCA评分及其子项评分与 HCY水平,LA评分及P300潜伏期、波幅的相关性.结果:PS-CIND组与卒中后无认知障碍组(NCI)组相比在M oCA总分、及子项视空间/执行功能、延迟回忆两子项得分上,有统计学差异(P<0.05),其HCY水平显著增高,LA评分明显严重,并存在明显的P300潜伏期延长和波幅下降,均具有统计学差异(P<0.05或 P<0.01).HCY水平、LA 评分与MoCA评分均呈负相关,相关系数 r分别为-0.652(P<0.05)、-0.714(P<0.05),M oCA评分与P300潜伏期、波幅分别成负相关和正相关,相关系数 r分别为-0.693(P<0.05),r=0.442(P<0.05).HCY水平与M oCA子项中的执行/视空间功能,注意力,延迟回忆,命名四项均呈负相关,相关系数 r分别为 -0.704,-4.21,-6. 44,-5.98,差异均具有统计学意义(P<0.05);LA评分与MoCA子项中的执行/视空间功能、延迟回忆呈负相关,相关系数 r分别为-0.651(P<0.05),-0.538(P<0.05);P300潜伏期与M oCA子项中的执行/视空间功能、注意评分呈负相关,相关系数 r分别为-0.814(P<0.05),-0.795(P<0.05),P300波幅与M oCA子项中的执行/视空间功能、延迟记忆呈正相关,相关系数 r分别为0.601(P<0.05),0.554(P<0.05).结论:PS-CIND的认知损害主要表现在视空间/执行功能的损害及延迟回忆的损害,P300潜伏期、波幅,HCY水平及LA评分可作为反应PS-CIND患者总体认知损害程度和选择性认知损害领域程度的客观的电生理,生化和影像学指标.