简介：摘要降秩回归是一种用于降维的多变量线性回归，近年来被越来越多的应用于营养流行病学领域来提取膳食模式。但Stata软件作为主流的统计分析工具之一，目前尚无进行降秩回归的程序。本文介绍了降秩回归方法的基本思想、实现降秩回归的基本原理，并开发和公开发布了功能全面的Stata程序（rrr），最后结合中国慢性病前瞻性研究数据进行示例分析。rrr程序的研发和公开，为该统计方法在流行病学和公共卫生领域的广泛应用提供了实用工具。

简介：摘要利用生物学年龄衡量个体的衰老进程，可以弥补实足年龄难以解释同龄个体间健康水平差异的缺陷。近年来，涌现出多种基于临床或表型、分子生物学指标或复合型指标的生物学年龄测量方法。本文对一些常见的生物学年龄评价指标进行概述，并比较这些评价指标的效度。

简介：摘要安全有效的疫苗是应对新型冠状病毒（新冠病毒）肺炎大流行的核心措施之一，从疫苗上市前安全性、保护效力的实验性研究到上市后安全性监测和人群保护效果的观察性研究，流行病学发挥了重要作用。本文通过列举国内外研究案例，系统介绍了临床试验、队列研究和病例对照研究设计在疫苗上市前后的应用及设计要素，以期丰富流行病学教学案例。另外，在新冠病毒疫苗上市后，国外多个国家在短时间内组织开展了规模数十万到千万人的人群观察性研究，得益于这些国家完善的公共卫生、医疗、医保等信息系统，通过不同数据系统的互联互通，实现整合分析。这对我国健康医疗大数据的发展具有重要启示。

简介：摘要目的描述中国成年人的共病流行情况及其长期变化，以及常见的共病模式。方法本研究纳入参与中国慢性病前瞻性研究第二次重复调查的25 033人，利用研究对象参加2004-2008年基线调查和2013-2014年第二次重复调查时采集的信息进行分析。根据自报疾病史、现场体格检查及血液样本检测等信息定义了13种慢性病或健康问题，将共病定义为同时患有≥2种慢性病或健康问题。采用系统聚类分析法描述共病模式。结果研究对象在基线调查时的年龄为（51.5±10.1）岁，第二次重复调查时为（59.5±10.2）岁。2次调查平均间隔时间（8.0±0.8）年，共病率由33.5%上升至58.1%，人均患病数由1.15个增加至1.82个，平均每5年增长0.42个患病数。年长者、城市人群、文化程度低者的共病率较高，且患病数量随年龄的增速更快。吸烟及过量饮酒者的患病数量随年龄增速也更快。该人群中最常见的共病组合为：超重肥胖、高血压、糖尿病、中风和冠心病。结论我国成年人群共病率较高，且随年龄增长而增加。共病情况在不同地区、文化程度水平和不同生活方式人群中存在差异。

简介：摘要目的比较Fried表型以及由不同数量的疾病缺陷构成的衰弱指数（FI）对衰弱状态评价的一致性以及与死亡风险的前瞻性关联。方法利用中国慢性病前瞻性研究（CKB）第二次重复调查的23 615名研究对象的数据，采用5种表型指标构建Fried表型，并分别纳入28个和40个疾病缺陷构建FI-28和FI-40。计算加权Kappa系数比较3种指标对衰弱状态分类的一致性。采用Cox比例风险模型分析衰弱指标与死亡风险的关联。结果采用Fried表型、FI-28和FI-40计算的衰弱率分别为5.4%、7.9%和4.0%。Fried表型与FI-28和FI-40的Kappa系数分别为0.357和0.408，FI-28与FI-40的Kappa系数为0.712。经过（3.9±0.5）年的随访，死亡755人。当采用Fried表型时，与无衰弱组相比，衰弱前期和衰弱组的死亡风险均增加，多因素调整后的风险比（HR）（95%CI）分别为1.60（1.32~1.94）和2.90（2.25~3.73）；采用FI-28时，衰弱前期和衰弱组的死亡HR值分别为1.71（1.39~2.11）和2.52（1.95~3.27）；采用FI-40时，衰弱前期和衰弱组的死亡HR值分别为1.98（1.60~2.44）和3.71（2.80~4.91）。衰弱状态与死亡风险的关联在不同年龄组间存在差异，在低年龄组中的关联强度高于高年龄组。结论Fried表型和基于不同数量的变量构建的FI表现出较好的一致性，都能较好地预测死亡风险。

简介：摘要目的了解苏州市居民体力活动水平，探索社区建成环境与居民体力活动的关联。方法2017年采用多阶段整群随机抽样方法抽取苏州市25~64岁常住人群进行面对面调查，采用国际体力活动量表长卷（IPAQ-L）评估居民体力活动水平，采用居民环境步行量表简表（NEWS-A）评价社区建成环境主观感知。结果苏州市居民过去1周总体力活动水平M=3 610.42 MET-min/w，以工作相关体力活动水平为主，交通、家务及休闲相关体力活动水平较低。控制社会人口学因素后，公共服务可及性与社区居民的总体力活动水平呈负相关（OR=0.522，95%CI：0.329~0.830），场所设施多样性与工作相关体力活动水平呈负相关（OR=0.701，95%CI：0.492~0.999），步行和自行车道设施与工作相关体力活动水平呈正相关（OR=1.603，95%CI：1.004~2.559）；交通安全与交通相关体力活动水平呈负相关（OR=0.642，95%CI：0.416~0.990）；住宅密度与休闲相关体力活动水平呈正相关（OR=1.001，95%CI：1.000~1.002）；此外，社区美观与舒适的主观感知程度越高，工作、交通、家务及总体力活动水平越高（OR=1.889，95%CI：1.176~3.033；OR=1.671，95%CI：1.120~2.495；OR=1.775，95%CI：1.143~2.756；OR=1.593，95%CI：1.079~2.350）。结论完善步行道和自行车道设施、提高社区的美观和舒适程度对于增加居民体力活动有重要作用。

简介：摘要目的分析身高降低与跟骨骨密度的关联。方法研究对象为参加中国慢性病前瞻性研究第二次重复调查者，此次调查时采用定量超声骨密度仪测定跟骨骨密度。本研究纳入分析24 231人。采用线性回归分析基线至重复调查的身高变化与重复调查时骨密度指标之间的关联。结果基线至重复调查平均间隔时间为（8.0±0.8）年，身高降低≥1.0 cm者占33.0%，≥3.0 cm者占3.7%。控制潜在混杂因素后，身高降低与骨密度水平的关联存在线性趋势（均P<0.001）。身高降低1.0 cm，宽带超声衰减、超声声速和强度指数分别变化-0.79（-0.95～-0.63）、-2.74（-3.35～-2.13）和-1.29（-1.54～-1.04）。与身高稳定者相比，身高降低≥3.0 cm者的宽带超声衰减、超声声速和强度指数的β值（95%CI）依次为-3.29（-4.08～-2.50）、-10.70（-13.66～-7.73）和-5.16（-6.36～-3.96）。亚组分析显示，上述关联在女性、年龄≥55岁和低体力活动者中更为显著。结论身高降低幅度越大，骨密度指标越低。身高作为容易测量的体格指标，应该定期关注其变化情况，以早期发现和预防骨质疏松症等骨骼健康问题。

简介：摘要目的描述中国慢性病前瞻性研究（CKB）项目10个地区非吸烟人群被动吸烟行为特征的地区和人群分布差异。方法分析10个地区317 486名30~79岁非吸烟者在基线调查时的被动吸烟行为信息，描述其被动吸烟行为特征的地区和人群分布差异。结果按照全国人口普查结果进行标化，我国成年人被动吸烟率为56.7%，与吸烟者共同居住率为66.5%，其中农村高于城市。周被动吸烟频率、周被动吸烟累计时长、日均被动吸烟时长具有明显的地区差异，城市地区的周被动吸烟累计时长随周被动吸烟频率的增加而增加，女性中湖南省的周被动吸烟频率较高，但周被动吸烟累计时长较低，日均被动吸烟时长最低，河南省却与之相反；与吸烟者共同居住者的被动吸烟率是从未共同居住者的2.27倍（95%CI：2.24~2.29），且在女性中关联更强（OR=2.61，95%CI：2.58~2.64），而在男性中二者无关（OR=1.01，95%CI：0.95~1.06）。女性除日均被动吸烟时长小于男性外，其余指标均大于男性；且年龄较小、文化程度较低者的多数被动吸烟指标较高；家庭年收入较低者的被动吸烟率和与吸烟者共同居住率较低，但日均被动吸烟时长较高；女性在婚者的被动吸烟率、与吸烟者共同居住率较高，男性与之相反。结论CKB项目10个地区非吸烟人群的被动吸烟率、周被动吸烟频率、周被动吸烟累计时长和日均被动吸烟时长以及与吸烟者共同居住率、共同居住年限均存在明显的地区和人群分布差异。

简介：摘要目的分析骨折住院率的三间分布。方法纳入参与中国慢性病前瞻性研究基线调查者，剔除2009年以前失访或死亡者，共纳入506 004人。采用负二项回归模型分析2009-2016年10个项目地区任意骨折和5个部位骨折（上肢、下肢、椎骨、骨盆和髋部骨折）住院率的时间、地区和人群分布特征。结果在（7.7±1.2）年（累积3 899 814人年）的随访期间，共记录到17 118例任意骨折住院事件，任意骨折的粗住院率为4.39/1 000人年。在控制了固定队列人群年龄增长后，2009-2016年间，不同部位的骨折住院率均呈上升趋势，任意骨折年增长率（95%CI）为9.1%（8.3%~9.9%）。除髋部骨折外，农村骨折住院率高于城市（P<0.05），农村和城市的任意骨折住院率分别为5.42/1 000人年和3.24/1 000人年。骨折住院率随年龄增加。<50岁时，除骨盆骨折外，男性骨折住院率高于女性；≥50岁时，女性骨折住院率均高于男性。结论骨折住院率随年龄增加的同时，也表现出长期增长趋势。我国正在迈入老龄化社会，骨折的疾病负担将会越来越重，预防骨质疏松和防止伤害相关骨折的发生，从而减少骨折发病率具有重要意义。

简介：摘要目的分析中老年人群衰弱状态随时间的变化以及衰弱恶化的影响因素。方法利用美兆健康体检中心的人群体检数据，纳入2008-2019年至少参加过2次体检且间隔>3年者，共13 689人。采用28个变量（包括疾病、症状、体格测量指标等）构建衰弱指数，衰弱指数≥0.25定义为衰弱，衰弱指数>0.10且<0.25为衰弱前期。按照首次体检时的衰弱状态进行分层分析，采用logistic回归模型探究社会人口学因素和生活方式特征对衰弱状态恶化的影响。结果研究对象首次和末次体检时的年龄分别为（42.3±9.2）岁和（47.9±9.3）岁，首末次体检间隔时间为（5.7±1.9）年。首次体检时人群衰弱率和衰弱前期率分别为2.5%和50.3%，末次体检时分别上升至3.9%和55.4%。67.3%的研究对象在首末次体检时的衰弱状态保持不变，21.2%恶化，12.5%改善。在首次体检衰弱状态为健康的人群中，年龄大、女性、低文化程度、戒烟和每天吸烟、BMI定义的肥胖、以及腰臀比定义的中心性肥胖是衰弱状态恶化的危险因素。在首次体检为衰弱前期的人群中，年龄大、女性、BMI定义的肥胖、以及腰臀比定义的中心性肥胖是衰弱状态恶化的危险因素。结论低文化程度、吸烟和肥胖等会增加个体衰弱恶化的风险。

简介：摘要目的探索我国成年人主要慢性病与肌肉重量、力量和质量的关系。方法主要利用中国慢性病前瞻性研究项目的第二次重复调查，采用logistic回归分析主要慢性病种类、个数及患病年限与肌肉重量、力量和质量的相关性。结果在24 242名研究对象中，糖尿病、冠心病、脑卒中和慢性阻塞性肺疾病（慢阻肺）的患病率分别为9.6%、5.8%、3.2%和26.8%，患以上任一慢性病者占38.8%。患上述任一慢性病与低手握力和低上肢肌肉质量（AMQ）呈正相关，且其患病年限越长，风险越高（线性趋势P值均<0.001），患任一慢性病≥10年者低手握力和低AMQ的OR值（95%CI）分别为1.64（1.42~1.90）和1.83（1.60~2.10）；患1、2、≥3种慢性病者低手握力的OR值（95%CI）分别为1.26（1.17~1.37）、1.42（1.23~1.64）和2.27（1.55~3.32），低AMQ的OR值（95%CI）分别为1.28（1.18~1.38）、1.67（1.46~1.92）和2.41（1.69~3.45），且关联均存在线性趋势（线性趋势P值均<0.001）。糖尿病、冠心病和脑卒中均与低手握力和低AMQ呈正相关；慢阻肺患者低四肢肌肉重量指数（ASMI）、低全身肌肉重量指数（TSMI）和低手握力的风险高于非慢阻肺患者，且风险与患病年限呈正相关（线性趋势P值均<0.05），患慢阻肺≥10年者低ASMI和低TSMI的OR值（95%CI）分别为2.39（1.36~4.20）和2.41（1.37~4.26）。结论慢性病患者肌肉力量和质量更低，且随患慢性病个数及年限延长，低手握力、低AMQ的比例更高。

简介：摘要目的探究我国成年男性和女性中慢性阻塞性肺疾病（COPD）与骨密度的关系。方法利用中国慢性病前瞻性研究第二次重复调查的数据，采用线性回归分析COPD及慢性阻塞性肺疾病全球倡议（GOLD）分级与跟骨骨密度指标（宽带超声衰减、超声声速、强度指数）之间的相关性。结果共纳入23 876名研究对象，男性COPD患病率为34.3%（3 130/9 125），女性COPD患病率为23.2%（3 416/14 751）。以健康男、女性为参照，男、女性COPD患者的宽带超声衰减、超声声速、强度指数均降低，男性对应的β值（95%CI）分别为-0.87（-1.37~-0.36）、-3.42（-5.42~-1.43）、-1.53（-2.34~-0.71），女性对应的β值（95%CI）分别为-0.66（-1.09~-0.23）、-2.24（-3.92~-0.55）、-1.06（-1.71~-0.40），男性中下降幅度更大（交互作用均P<0.05）。GOLD分级与宽带超声衰减、超声声速、强度指数呈负相关，且存在线性趋势（线性趋势检验均P<0.001）。亚组分析显示，男性年龄≥60岁、男性低体力活动水平者、男、女性非超重/肥胖者中COPD患者骨密度下降幅度更大。结论COPD与骨密度呈负相关，应加强对COPD患者骨密度的关注，尤其应关注年龄大、体力活动水平低、非超重/肥胖的COPD患者。

简介：摘要目的探索我国成年人体力活动、休闲静坐时间与低肌肉重量、力量和质量的关系。方法利用中国慢性病前瞻性研究（CKB）项目的第二次重复调查，采用logistic回归分析体力活动水平、休闲静坐时间与低四肢肌肉重量指数（ASMI）、低全身肌肉重量指数（TSMI）、低肌肉力量和低肌肉质量（AMQ）的相关性。结果研究纳入24 245名研究对象，体力活动水平为（18.3±13.8）MET-h/d，休闲静坐时间为（4.4±1.9）h。以体力活动最低组为参照组，高水平的体力活动与较低的低ASMI、低TSMI、低手握力和低AMQ有关，对应的OR值（95%CI）分别为0.68（0.60~0.77）、0.66（0.58~0.75）、0.82（0.72~0.94）和0.84（0.74~0.95）。工作、交通、家务和休闲相关体力活动亚型亦与上述指标呈现不同程度的负相关。与休闲静坐时间最短组相比，静坐时间最长组低TSMI的危险性增加（OR=1.13，95%CI：0.99~1.30）。结论体力活动与低肌肉重量、力量和质量呈负相关，而休闲静坐时间与低肌肉重量呈正相关。

简介：摘要目的探索我国成年人睡眠时长及失眠症状与低肌肉重量、力量和质量的关系。方法主要利用中国慢性病前瞻性研究（CKB）项目的第二次重复调查，采用logistic回归分析睡眠时长、有无失眠症状及其年限与低肌肉重量、力量和质量的相关性。结果研究人群每天睡眠时长为（7.4±1.5）h，9.3%睡眠时间过短（<6 h），17.4%睡眠时间过长（≥9 h），29.9%自报有失眠症状。以睡眠时间7~h为参照组，睡眠时间过长与低肌肉重量、力量和质量呈正相关，低四肢肌肉重量指数（ASMI）、低全身肌肉重量指数（TSMI）、低手握力和低上肢肌肉质量（AMQ）的OR值（95%CI）分别为1.32（1.18~1.48）、1.26（1.13~1.41）、1.33（1.18~1.49）和1.16（1.03~1.30）。与无失眠症状者相比，自报有失眠症状者低ASMI和低TSMI的危险性均增加，且其失眠症状年限越长，低ASMI和低TSMI的风险越高（线性趋势均P<0.001）。与睡眠时间7~h且无失眠症状者相比，睡眠时间过短且有失眠症状者低ASMI和低TSMI的风险明显增加，OR值（95%CI）分别为1.26（1.08~1.47）和1.25（1.07~1.46）。结论睡眠时间过长者肌肉重量、力量和质量更低，有失眠症状者肌肉重量更低，且随失眠症状年限增加，低ASMI和低TSMI的比例越高。

简介：摘要目的探索我国成年人主要膳食模式与低肌肉重量、力量和质量的关系。方法基于中国慢性病前瞻性研究（CKB）项目第二次重复调查，对20个食物组进行因子分析提取膳食模式。将四肢肌肉重量指数（ASMI）、全身肌肉重量指数（TSMI）、手握力和上肢肌肉质量（AMQ）分性别的最低五分位数作为判断低肌肉重量、力量和质量的截断值。采用logistic回归分析膳食模式与低肌肉重量、力量和质量的相关性。结果提取出2种主要膳食模式，平衡膳食模式以摄入多种食物为特征，米肉膳食模式以摄入大米、肉类、家禽和鱼类较多为特点。平衡膳食模式得分高者，低TSMI、低手握力和低AMQ的危险性降低，对应OR值（95%CI）分别为0.83（0.74~0.95）、0.64（0.56~0.74）和0.82（0.72~0.93），且关联具有线性趋势（P<0.05）；米肉膳食模式得分高者，低ASMI、低TSMI和低手握力的危险性更低，OR值（95%CI）分别为0.67（0.55~0.82）、0.69（0.56~0.85）和0.74（0.60~0.91），且关联具有线性趋势（P<0.05）。结论平衡膳食模式和米肉膳食模式得分较高者的肌肉重量、力量和质量更高。

简介：摘要目的比较分析5种甲基化年龄加速指标与生活方式和心血管代谢因素间的相关性。方法基于中国慢性病前瞻性研究中有基线外周血全基因组甲基化检测数据的研究对象980人，计算Horvath时钟、Hannum时钟、DNAm PhenoAge、GrimAge和Li时钟5种甲基化年龄。根据甲基化年龄对实足年龄回归的残差值计算甲基化年龄加速。研究的生活方式因素包括吸烟、饮酒、饮食习惯、体力活动、经BMI和腰围联合定义的体型。心血管代谢因素包括血压、血糖和TC。利用一般线性模型分析生活方式和心血管代谢因素与各甲基化年龄加速指标的相关性［β值（95%CI）］。结果GrimAge加速指标与吸烟、饮酒、体力活动水平及BMI存在关联。与不吸烟、不饮酒或BMI为18.5~23.9 kg/m2 者相比，吸烟者（每天吸烟1~14、15~24、≥25支者对应的β值依次为0.71（95%CI：0.57~0.86）、0.88（95%CI：0.73~1.03）、0.99（95%CI：0.81~1.18）、重度饮酒者［每日纯乙醇量≥60 g：0.33（95%CI：0.11~0.55）］、BMI<18.5 kg/m2者［0.23（95%CI：0.03~0.43）］表现为加速衰老；与体力活动水平低者相比，中、较高、高体力活动者表现为减速衰老［β值依次为-0.14（95%CI：-0.27~-0.00）、-0.12（95%CI：-0.26~0.02）、-0.13（95%CI：-0.26~0.01）］。随着健康生活方式数量的增加，GrimAge加速指标呈现下降趋势（P<0.001）；与具有0~1个健康生活方式者相比，具有2、3、或4~5个健康生活方式者的β值依次为-0.30（95%CI：-0.47~-0.12）、-0.47（95%CI：-0.65~-0.30）、-0.72（95%CI：-0.90~-0.53）。其他4个指标与多数生活方式因素不存在有统计学显著性的关联。5种甲基化年龄加速指标与血压、血糖和TC均不存在关联。结论生活方式不健康者表现出表观遗传年龄的加速，即DNA甲基化预测年龄老于实足年龄。