目的 研制并评价老年人社会支持评价量表简化版（the social support appraisals scale, SSA）。方法 采用整群随机抽样法，在东莞市选取1130名符合纳入标准的老年人进行一对一访问调查，随机抽取其中565名被试对象作为训练样本，用经典测量理论和项目反应理论的方法对训练样本进行条目筛选；剩下565名被试对象作为考核样本，用于评价简化版量表的信效度。结果 每个维度保留了4个较优条目（朋友维度保留了5个条目）,形成了13个条目的简化版量表SSA-13。SSA-13 Cronbach′s α系数为0.903,各维度的Cronbach′s α系数介于0.755～0.866之间。验证性因子分析拟合度指标为：卡方自由度比值（χ²/df）为4.93,比较适配优度指数（CFI）为0.94,增值适配指数（IFI）为0.94,比较适配度指数（GFI）为0.92,非标准适配指数（TLI）为0.92,渐进残差均方和平方根为0.08,残差均方和平方根（RMR）为0.02,各项拟合度指标均达到理想指标。项目反应理论分析结果显示，各条目的区分度介于2.34～3.10之间，除条目SS-A18难度参数（-5.56～1.23）外，其他条目难度参数介于在-3.76～2.25,平均信息量在0.764～1.311。与原量表相比，总量表的平均信息量减少了16.06%。结论 SSA-13具有良好的信度和效度，可以作为老年人社会支持评价研究中简单易行、可靠有效的测量工具。

目的 分析脑血管病死因对我国城市不同性别、年龄人口期望寿命的影响。方法 应用蒋庆琅法，用excel软件计算和编制2019年全国城市人口分性别、年龄简略寿命表，去脑血管病死因寿命表，再运用最终死于某死因概率计算脑血管病对期望寿命的影响程度。结果 2019年我国城市人口出生时期望寿命为84.81岁(男性81.15岁，女性89.10岁),因脑血管病损失的预期寿命为4.9岁(男性4.22岁，女性5.7岁),男性尚存人数曲线下降较女性更陡峭，死亡高峰较女性更提前。男性因脑血管病死亡概率最高的为80～岁年龄组(0.237395),女性最高为70～岁年龄组(0.211967),各年龄组比较脑血管病占男性全死因比例均高于女性，且死亡人群集中年龄较女性更年轻。结论 因脑血管病损失的期望寿命在不同性别、年龄人群中呈现差异性，在男性中老年人群中影响更明显且呈年轻化趋势，卫生部门应针对此制定相应健康策略，以期延长人口寿命。

围绕医学影像的人工智能医疗器械，介绍其临床试验的统计学设计要点，主要涉及临床试验设计类型确定、对照设置、比较类型选择、主要结局指标确定、样本量估计方法、随机化策略、盲法实施、缺失值处理及敏感性分析等内容。本共识旨在为从事该类医疗器械临床试验的相关人员提供统计学设计参考，以保证临床试验结果的科学性、可靠性及规范性。

目的 分析1990—2019亚洲主要国家及全球宫颈癌疾病负担，为制定宫颈癌防控策略提供科学依据。方法 收集1990—2019年全球疾病负担研究(global burden of disease, GBD)数据，利用发病、死亡、伤残调整寿命年(disability-adjusted life years, DALY)等指标定量评估宫颈癌的疾病负担情况，采用估计百分百变化(estimated annual percent change, EAPC)分析疾病负担时间变化趋势。结果 1990—2019年，全球宫颈癌标化发病率从7.64/10万下降至6.81/10万，标化死亡率从4.46/10万降低至3.40/10万，标化DALY率从139.98/10万下降至107.20/10万，年平均降低分别为0.39%、0.96%、0.94%(EAPC<0,P<0.05)。亚洲主要国家中，中国的标化发病率从4.20/10万上升至5.53/10万，平均每年增长1.63%(EAPC>0,P>0.05)。中国宫颈癌的死亡数从2.64万增加到5.34万，相对增长1.02倍，全球及亚洲主要国家标化死亡率、标化DALY率均呈逐年下降趋势，而中国未出现下降趋势。相关分析显示，发病率与社会人口指数(socio-demographic index, SDI)无明显相关性(ρ=-0.13,P=0.11),死亡率与SDI呈负相关(ρ=-0.74,P<0.001),DALY与SDI也呈负相关(ρ=-0.77,P<0.001)。结论 中国宫颈癌的标化发病率、死亡率逐年增加，疾病负担呈上升趋势，应积极采取科学有效的防控措施，降低宫颈癌疾病负担。

目的 调查新冠疫情防控校园封闭管理大背景下医学生焦虑的影响因素及领悟社会支持的中介效应。方法 2022年4月至5月，采用自主学习能力测量量表、焦虑量表、领悟社会支持量表对沈阳市某高校1554名医学生展开问卷调查。结果 医学生焦虑检出率为20.8%,居住地、健康状况、参与志愿活动是医学生焦虑的影响因素(P<0.05)。自主学习能力与焦虑呈显著负相关(r=-0.115,P<0.01),与领悟社会支持呈正相关(r=0.137,P<0.01),焦虑与领悟社会支持呈负相关(r=-0.478,P<0.01)。医学生的自主学习能力直接影响焦虑(β=-0.127,P<0.01),领悟社会支持在自主学习能力与焦虑的关系中发挥部分中介作用(95%CI-0.090～-0.039)。结论 学校和家庭应及时关注医学生的心理健康状况，给予全方位的支持，同时注重自主学习能力的提升。

目的 对基于分组的轨迹模型(group-based trajectory modeling, GBTM)的基本原理、适用性及SAS软件程序进行介绍。方法 使用SAS软件中的PROC TRAJ过程步，构建不同数据类型的GBTM模型。结果 选取多项前瞻性数据，根据数据类型构建了CNORM、ZIP等单轨迹及双轨迹模型，在此基础上将协变量纳入模型，观察影响发展轨迹的因素。结论 GBTM为解决纵向数据复杂性和结果可解释性之间的平衡问题提供了有价值的工具，同时SAS软件可以很好地实现轨迹模型的构建，并为模型的不断拓展提供了较为方便的实现形式。

目的 对2012—2021年江苏省13个市的村卫生室医疗资源配置效率进行评价。方法 运用数据包络分析(data envelopment analysis, DEA)模型和Malmquist指数分析江苏省村卫生室医疗资源配置效率。结果 2021年江苏省13个市中南京市、无锡市等11个市为DEA有效，淮安市和扬州市为DEA无效。2012—2021年江苏省村卫生室医疗资源Malmquist指数均值为1.006,10个市(无锡市、徐州市、常州市、苏州市、南通市、连云港市、淮安市、盐城市、扬州市、镇江市)全要素生产率指数大于1,其中5个市(无锡市、徐州市、常州市、苏州市、南通市)全要素生产率指数高于平均水平，3个市(南京市、泰州市、宿迁市)全要素生产率指数小于1。江苏省村卫生室医疗资源全要素生产率主要受技术进步指数的影响。结论 江苏省村卫生室医疗资源存在地区差异，部分地区尚需优化，村卫生室人力资源配置趋于合理，可以通过加大推动技术进步的投入、合理配置各级医疗资源、升级村卫生室医疗设备、加强卫生信息化建设、提高村卫生室医疗资源的利用效率。

目的 评估希望区域法应用于适应性设计样本量再估计的准确性及稳健性，为该方法的适用条件提供理论参考。方法 以二分类资料为例，以包含期中分析的两阶段适应性试验为框架，采用Monte Carlo方法，在相同的组间响应率率差、同合并率等参数的模拟场景下，比较固定样本量设计、成组序贯设计和希望区域法进行样本量再估计的准确性和稳健性。结果 模拟研究显示希望区域法在期中检验效能落入希望区域内进行适应性样本量调整时，Ⅰ类错误与同场景下固定设计和成组序贯设计近似。初始估计样本量被低估时，希望区域法比固定设计的检验效能平均提高约5%,比成组序贯设计提高约8.8%,样本量比固定设计多耗费0.18倍，比成组序贯设计多0.38倍。高估初始估计样本量时，三种设计的检验效能相差仅在1%左右，但希望区域法比固定设计多消耗7.5%、比成组序贯设计多耗费48%的样本量。结论 相比于固定样本量及成组序贯设计，希望区域法对于试验的整体检验效能提高不大，其平均样本量均大于相同设定场景下的固定设计及成组序贯设计，使用该设计方法时应谨慎权衡试验收益。

目的 探讨老年食管癌患者术后医院感染危险因素并构建术后医院感染风险评分预测模型，为预防老年食管癌患者术后感染提供参考依据。方法 选取2019年5月至2021年5月行手术治疗老年食管癌患者236例，根据术后医院感染情况分为感染组、未感染组，通过多因素logistic回归分析老年食管癌术后感染危险因素并根据回归系数对所得危险因素赋值建立术后医院感染风险预测评分模型，使用受试者工作特征曲线(receiver operating characteristic, ROC)、Hosmer-Lemeshow检验评估该模型的预测效能(区分度和拟合度)。结果 236例老年食管癌患者术后医院感染50例，术后医院感染发生率为21.19%。糖尿病、术前中性粒细胞水平≥75%、留置胸腔引流管时间≥8 d为老年食管癌患者术后医院感染危险因素，风险预测评分模型ROC曲线下面积为0.929,Hosmer-Lemeshow检验P=0.236。结论 根据老年食管癌患者术后医院感染危险因素建立的预测评分模型对术后医院感染有较高的预测风险价值，便于快速识别术后医院感染高风险人群。

目的 基于潜变量增长混合模型(latent growth mixture model, LGMM)和中国健康与养老追踪调查(China Health and Retirement Longitudinal Survey, CHARLS)3期追踪数据，考察我国老年人认知水平的发展轨迹及影响因素。方法 资料来源于CHARLS 2013年、2015年、2018年数据，共包含4662例60岁及以上老年人。采用MMSE量表对老年人的认知水平进行评估。利用Mplus工具构建适用于队列数据的LGMM模型探讨不同老年人群体的认知水平变化轨迹，考察时间恒定因素(性别、年龄、受教育程度、婚姻状况、吸烟、饮酒)和随时间变化因素(日常生活活动能力、睡眠时长、主观生活满意度)对老年人认知水平的影响。结果 根据相关指标可将老年人认知水平的发展趋势划分为3个类型，分别是C1缓慢下降组(35.2%)、C2快速下降组(20.6%)和C3稳定组(44.2%);老年人的认知水平随着时间推移呈整体下降趋势，3个组别的认知水平发展轨迹均有显著的时间效应，C2快速下降组的认知水平随时间下降较快，C1缓慢下降组的认知水平随时间缓慢降低，C3稳定组的认知水平变化不明显，呈稳定态。结论 老年人认知水平的发展轨迹具有显著的群体异质性和个体差异性，存在不同类别的发展趋势。

<正>2024年7月19日至21日,2024年中国生物统计学术会议在广东广州召开。会议由中华预防医学会生物统计分会、中国卫生信息与健康医疗大数据学会统计理论与方法专业委员会、中国卫生信息与健康医疗大数据学会卫生统计学教育专业委员会、中国现场统计研究会医药与生物统计分会联合主办,中山大学公共卫生学院、北京大学公众健康与重大疫情防控战略研究中心联合承办,中国统计教育学会生物医学统计分会、中国医药教育协会医药统计专业委员会、广东省卫生统计学会协办。

<正>组学数据在分子层面上加深了研究者对表型相似疾病的辨别能力及对认识不足疾病的理解，对组学数据的分析能够在疾病的诊断与治疗中发挥至关重要的作用。组学分子间存在复杂的调控关系，过去的研究表明细胞的功能是模块化的^[1],在疾病的实际进展中，关键的调控网络出现异常往往表现为网络中所有分子表达水平都发生变化^[2],因此表达模式相似的分子很可能存在共调控或功能相关的情况^[3]。

目的 大量环境流行病学研究表明PM_2.5浓度短期变化与人群每日死亡相关，但是大多数关于PM_2.5对每日死亡的急性效应的研究是通过回归来控制混杂因素，而通常可用的人群资料里仅包含少量测量的混杂因素，这就面临着大量未观测的混杂因素未纳入模型的问题从而导致估计有偏，而工具变量法可以较好地解决未观测混杂带来的效应估计问题，本文通过工具变量法估计PM_2.5对每日死亡的急性效应。方法 收集中国某市2016—2019年PM_2.5日均浓度、气象数据及每日非意外死亡人数，采用边界层高度和风速作为工具变量，分析该地PM_2.5浓度对每日非意外死亡的影响；采用阴性暴露对照法检验工具变量假设；采用时间序列的bootstrap方法估计置信区间。并与广义相加模型的效应估计进行比较。结果 工具变量法得出PM_2.5浓度与人群每日非意外死亡存在相关，PM_2.5浓度每升高10μg/m³,人群每日非意外死亡增加0.94%（95%CI:0.39%～1.55%）;阴性暴露对照结果显示阴性暴露与人群日非意外死亡不相关（P=0.19）,说明前述工具变量模型不受未测量且未控制混杂的影响。传统广义相加模型估计PM_2.5浓度每升高10μg/m³,人群每日非意外死亡增加0.24%（95%CI:0.01%～0.47%）。结论 经工具变量法估计，该地PM_2.5浓度与居民每日非意外死亡存在相关，边界层高度和风速可以作为工具变量估计PM_2.5浓度对人群每日非意外死亡的急性效应。

目的 描述辽宁省城市脑肿瘤患者10年生存率并分析其影响因素，为评价辽宁省肿瘤防治工作提供参考依据。方法 该研究对辽宁省肿瘤发病报告数据库中2000—2002年确诊的292例脑肿瘤患者进行主动随访及被动随访，使用寿命表法计算观察生存率（OSR）,EdererⅡ计算期望生存率（ESR）,相对生存率等于观察生存率除以期望生存率。使用Kaplan-Meier法绘制生存曲线，应用Log Rank Test法进行单因素分析，Cox比例风险模型进行多因素分析。结果 辽宁省城市脑肿瘤患者10年OSR、RSR和T₅₀分别为25%、28.6%和7.37个月。女性10年OSR、RSR、T₅₀（30.5%、34.6%、15.1个月）,高于男性10年OSR、RSR、T₅₀（20.5%、23.6%、14.2个月）;年龄在64岁及以下的脑肿瘤患者10年OSR、RSR、T₅₀（31.2%、32.9%、20.3个月）显著高于65岁以上患者组10年OSR、RSR、T₅₀（11.1%、19.2%、8.1个月）;Ⅰ～Ⅱ期患者10年OSR、RSR、T₅₀（41.1%、46.6%、52.5个月）显著高于Ⅲ～Ⅳ期患者10年OSR、RSR、T₅₀（6.45%、7.5%、6.0个月）,手术组10年OSR、RSR、T₅₀（35.3%、38.9%、31.8个月）显著高于非手术组10年OSR、RSR、T₅₀（10.2%、13.1%、5.3个月）,多因素结果显示，患者的确诊年龄（HR=1.693,95%CI:1.236～2.320）、临床分期（HR=3.535,95%CI:2.430～5.143）、确诊后采取的治疗方式（HR=2.322,95%CI:1.706～3.161）是影响辽宁省脑肿瘤患者生存率预后的独立因素。结论 辽宁省城市脑肿瘤患者的10年生存率在各国家地区处于中等水平，与发达国家相比还存在一定差距，提高脑肿瘤早诊早治率，加快研究脑肿瘤疾病的早期筛检技术，重视以手术治疗为核心，多学科综合治疗的治疗方式是未来提高脑肿瘤患者生存率的方向。

目的 探索构建门急诊医生服务评价量表，为提升患者满意度和三级公立医院绩效考核提供参考。方法 基于SERVQUAL模型及国内外研究量表进行编制，并使用德尔菲法、折半信度、Cronbach′s α系数、探索性及验证性分析法等方法评价量表的信度和效度，予以修订形成最终版量表。结果 量表包含有形性、响应性、可靠性、保证性和移情性五个维度共计24个条目。该量表折半信度为0.937,内部一致性为0.983,量表总分、量表各条目得分与医生服务总体满意度及医院总体满意度得分的Spearman相关系数较高。提取的5个公因子可解释总变异的92.427%,且各成分条目与量表理论构想基本一致，参照探索性因子分析结果调整量表后结果显示χ²/df=1.991,RMSEA=0.087,CFI=0.941,IFI=0.942,TLI=0.931,表明量表信度、效度较好。结论 该量表的信度与效度均较高，适合在大型医院中推广使用。

目的 建立季节性自回归移动平均模型（seasonal auto-regressive integrated moving average, SARIMA）,并探讨其对山西省百日咳发病的预测效果。方法 用R 4.1.3软件对2010—2019年山西省百日咳月报告发病率进行SARIMA模型构建与拟合，用所拟合的最优模型对2020—2022年山西省百日咳月发病率进行预测并与实际值比较，分析该模型的预测效果。结果 2010—2021年山西省共报告百日咳病例1956例，年平均发病率为46.50/1000万。发病率总体上呈现上升趋势，2021年时达到峰值140.67/1000万，2010—2019年百日咳报告发病主要集中在夏秋两季，2020年百日咳报告发病集中于上半年，2021年集中于下半年。2020年发病率与2019年相比减少了72.98%。SARIMA（5,2,2）(0,1,2)₁₂模型对于短期百日咳疫情有较好的预测性能。结论 SARIMA（5,2,2）(0,1,2)₁₂模型对于山西省短期百日咳发病有较好的预测性能，对于指导百日咳疫情防控有一定积极意义。

目的 了解南京市水痘的发病趋势，探究防控的重点地区，比较SaTScan与FleXScan法在空间聚集性研究中的差异和应用效果。方法 利用R 4.2.0开展乡镇/街道水痘报告发病率全局空间自相关分析，计算Moran′s I。分别利用SaTScan 9.4.1和FleXScan 3.1.2软件进行空间聚集性扫描，并对扫描结果进行比较分析。结果 自2017年纳入丙类管理之后，南京市水痘报告发病水平快速上升，2019年报告18 635例，无死亡病例报告，报告发病率为220.89/10万。全局空间自相关分析Moran′s I=0.20(Z=3.48,P<0.05)。SaTScan法和FleXScan法分别探测到10个(25个街道)和7个(24个街道)空间聚集区域(P均<0.05),发生聚集的街道范围基本一致，除南部高淳区的淳溪街道、古柏和溧水区的永阳街道外，其他聚集区域都分布于主城区周围的城郊地区，但聚集区域形状和每个聚集区域覆盖的街道有较大差异，SaTScan法聚集区覆盖了部分实际发病水平较低的街道，还将每个形状不规则的聚集区域探测为2个聚集区域。结论 建议将水痘纳入法定传染病进行报告管理。南京市水痘存在空间聚集性，两种方法均可在实践中选择使用，但FleXScan法更合理准确。

目的 分析1990—2019年中国人群精神障碍的流行状况、疾病负担变化趋势及危险因素。方法 从2019年全球疾病负担研究(global burden of disease, GBD)中获取中国人群精神障碍疾病负担指标，分析1990—2019年中国人群精神障碍患病率、伤残调整生命年(disability adjusted of life years, DALY)和归因危险因素DALY,采用Joinpoint 4.7估计精神障碍DALY率的年度变化百分比(annual percent change, APC)和平均年度变化百分比(average annual percent change, AAPC),分析精神障碍疾病负担的变化趋势。结果 2019年中国人群精神障碍患病率为11260.42/10万，抑郁症(3519.26/10万)和焦虑症(3363.60/10万)患病率最高。2019年女性DALY率(1564.22/10万)高于男性(1294.51/10万)。抑郁症(531.65/10万)、焦虑症(321.60/10万)和精神分裂症(43.15/10万)的DALY率最高，分别占2019年中国人群精神障碍总DALY率的37.26%、22.54%和17.59%。1990—2019年中国人群精神障碍DALY率由1321.15/10万增长至1426.78/10万，年龄标化DALY率呈下降趋势，平均每年下降0.19%(95%CI:-0.22%～-0.15%),但进食障碍年龄标化DALY率的增幅明显(AAPC=1.94%,95%CI:1.85%～2.02%)。2019年中国人群抑郁症归因于亲密伴侣暴力的DALY率最高(26.72/10万),1990—2019年男性抑郁症归因于童年性虐待的DALY率呈明显增长趋势(AAPC=0.68%,95%CI:0.26%～1.09%)。结论 应重视精神障碍疾病的防治，开展重点人群筛查并采取干预措施，针对危险因素开展病因预防，从而降低精神障碍导致的疾病负担。

目的 探讨适应性富集设计中的效应点估计方法，为富集设计下的疗效估计提供参考。方法 通过模拟适应性富集设计中的疗效数据，比较极大似然方法、第二阶段估计方法、Lindley估计方法、Bootstrap参数方法和条件无偏估计方法这五种方法在估计药物疗效方面的统计学性能。结果 极大似然估计方法往往产生较大偏倚，且容易受参数影响；Lindley估计方法和Bootstrap参数方法偏倚较小，但Lindley估计方法可能低估药物疗效；第二阶段估计方法和条件无偏估计方法的估计几乎无偏。除了第二阶段估计方法外，其他方法的均方误差均较小。结论 在适应性富集设计中，条件无偏估计方法的偏倚较小，均方根误差在可接受范围内，因此推荐使用该方法来估计药物疗效。

目的 已有研究证实体力活动及久坐与健康状态密切相关。但个体体力活动水平和久坐时间在生命历程中并非一成不变，常常随年龄增加发生变化。本研究旨在识别中国成年人体力活动(PA)与久坐时间(ST)轨迹并探索与全因死亡率之间的关系。方法 使用2004年至2015年中国健康与营养调查的数据，采用多变量组基轨迹模型中删失正态模型识别随生命历程变化的体力活动与久坐时间轨迹模式，使用多变量Cox比例风险回归模型评估体力活动与久坐时间轨迹分组与全因死亡风险的关联。结果 共纳入10576名具有2轮及以上完整随访数据的研究对象进行分析，本研究识别出4种体力活动与久坐时间轨迹。34.9%表现出PA水平下降，ST维持中等水平的特点，其余为低PA低ST(13.5%)、中PA高ST(32.8%)、高PA高ST(13.8%)。中PA高ST人群(HR=0.54,95%CI:0.37～0.78,P=0.001)、高PA高ST人群(HR=0.27,95%CI:0.11～0.65,P=0.002)死亡风险明显低于低PA低ST者，降低PA中ST组死亡风险差异没有统计学意义(HR=0.82,95%CI:0.59～1.16,P=0.257)。结论 本研究识别出了4种中国成年人PA与ST轨迹，并强调长期保持较高PA水平对降低死亡风险的重要性。

目的 利用TCPA及UCSC Xena数据库获取样本信息，构建胃腺癌预后相关蛋白的预测模型，并检验其应用价值，以期为临床治疗与相关机制研究提供参考。方法 利用R语言软件，通过单因素Cox分析和LASSO-Cox回归分析筛选与胃腺癌预后有关的蛋白质，并计算风险评分，根据风险评分结果将患者划分为高风险组和低风险组开展生存分析。同时绘制风险评分的ROC曲线、蛋白质热图等，用以评估这些蛋白质的预测能力。然后，将风险评分与患者其他临床表型数据结合，进行多因素回归分析，建立预后模型，并绘制列线图。通过C-index指标、ROC曲线、决策曲线等评估预后模型建立的准确性。进行GO富集分析与KEGG富集分析，探索相关蛋白所参与的生物学过程。结果 通过单因素Cox分析及LASSO-Cox分析筛选出15个差异蛋白质(CKIT、CHK1、CLAUDIN7、COLLAGENVI、DVL3、EGFR＿pY1173、ERALPHA＿pS11、NFKBP65＿pS536、SYK、ETS1、MYOSINIIA＿pS1943、P21、P90RSK、RAPTOR、XBP1),用于计算蛋白质风险评分。通过多因素回归分析，筛选出年龄、性别、M分期、N分期、风险评分这5种因素，构建预后模型。生存分析结果认为，蛋白质风险评分划分的高低风险组之间的总生存率有显著差异(P<0.0001)。列线图可以用来预测患者3年、5年的生存率，并认为具有较好的预测价值(C-index=0.7257)。ROC曲线显示模型在预测中表现出稳定的灵敏度和特异度(AUC=0.77和0.86)。结论 胃腺癌预后相关蛋白质的预测模型可用来预测患者的临床预后，研究结果可为个性化治疗和管理提供参考。

目的 应用一种可以同时解决少数类和多数类类间和类内不平衡问题的类别不平衡处理方法，并将其与随机森林(random forest, RF)分类器结合实现对弥漫大B细胞淋巴瘤(diffuse large B-cell lymphoma, DLBCL)患者早期复发的预测，为DLBLC患者的治疗提供参考。方法 首先使用一种基于高斯混合模型双向聚类重采样的类别不平衡处理方法(Gaussian mixture model, GMM-GMM)处理数据，并与随机过采样(random over sampling, ROS)、合成少数类过采样技术(synthetic minority over-sampling technique, SMOTE)、Borderline-1 SMOTE、Borderline-2 SMOTE、GMM上采样、GMM下采样、SMOTE+RUS、SMOTE+GMM和GMM+RUS进行比较，然后以RF作为分类器验证10种类别不平衡方法的性能，之后为验证RF的性能，在处理后的数据集上使用logistic回归和决策树(decision tree, DT)作为对照，最后从区分度和校准度两方面对模型进行评价。结果 在本文所有模型中，采用GMM-GMM的RF模型取得了相对最优的分类性能(accuracy=0.79, AUC=0.87, sensitivity=0.71, specificity=0.87, G-means=0.79, MSE=0.21)。结论 GMM-GMM优于其他传统的重采样方法，结合RF用于DLBCL患者早期复发的预测取得了相对较好的分类结果，可以很好地实现对DLBCL患者早期复发的预测。

目的 使用SAS以及操作系统的“任务计划”功能实现医学研究数据实时核查。方法 通过SAS直接访问数据集，对制定好的逻辑核查类型编写宏程序。结果 采用自行编制的SAS宏程序及其调用实例、SAS设置邮件发送和操作系统自动化运行SAS,完成缺失值以及异常值的核查。结论 当前研究所提供的SAS宏程序和系统设置，可帮助医学研究实现简单的数据核查功能，保证研究数据质量。

<正>信息技术的突飞猛进，带动着教学手段的改革。国家精品课程、精品资源共享课、MOOC等教学手段颠覆了学生只能在课堂学习的教学模式。以不同方式进行的混合式教学正在逐渐替代传统课堂教学。混合式教学主张将传统的面对面教学的优势与数字化教学的优势结合起来，二者优势互补，从而取得更佳的教学效果，成为国内外高等教育发展的重要趋势^[1]。

目的 研究真实世界研究中连续型结局变量下多个协变量缺失的数据模拟方法，以期为缺失值填补方法的模拟研究提供参考。方法 基于多元线性回归模型和三种缺失机制的定义，模拟了存在混杂因素情况下连续型结局变量的完整数据，设置了完全随机缺失、随机缺失、非随机缺失机制下3个协变量均缺失一定比例(如2%)的情境，并采用SAS软件实现缺失机制的模拟过程。结果 构建了不同缺失机制下多元协变量缺失的模拟数据集。结论 提供的数据模拟方法和SAS代码实现了三种不同缺失机制的多元缺失数据模拟，可为相关模拟研究提供参考。

目的 调查护士工作幸福感现状，并分析其影响因素。方法 采用组织支持感量表、工作幸福感量表、情感支持量表、工作旺盛感量表对1127名护士进行调查。结果 相关分析表明组织支持感与情感支持、工作旺盛感以及工作幸福感两两之间呈现显著正相关(P<0.001);结构方程模型表明组织支持感可以通过情感支持(95%CI:0.055～0.171)与工作旺盛感(95%CI:0.174～0.293)的中介作用影响护士的工作幸福感，也可以通过情感支持经由工作旺盛感的链式中介作用(95%CI:0.072～0.146)影响护士的工作幸福感。结论 情感支持与工作旺盛感在组织支持感与护士工作幸福感之间起链式中介作用。

目的 探讨智能体模型(agent-based model, ABM)在艾滋病流行病学中的应用价值。方法 介绍ABM的基本原理和模型特点，基于NetLogo构建以异质性性行为为主体的人类免疫缺陷病毒(human immunodeficiency virus, HIV)传播模型，在不同HIV检测频率的条件下估算男男性行为(men who have sex with men, MSM)中的HIV流行参数。结果 模拟结果显示随着HIV检测频率的增加，年新增感染数和诊断比例明显改善，且当HIV检测频率为1次/年时，感染HIV并知晓病情的比例为91%。结论 智能体模型适用于具有复杂传播网络的健康体系，可用于预测疾病传播趋势和比较不同政策组合的长期健康收益。

目的 利用带有缺失基因型观测的家族关联数据（correlated family data, CFD）和全基因组关联研究（genome-wide association study, GWAS）探索阿尔茨海默病的潜在致病基因及关联强弱。方法 研究人群来自华盛顿高地-英伍德哥伦比亚老龄化项目（Washington Heights-Inwood Columbia Aging Project, WHICAP）,该项目收集了先证者的基因型信息，并通过调查访谈收集了先证者及其家庭成员的人口统计学信息。本研究纳入352名先证者和820名关联家庭成员共1172人，首先利用家族结构信息和孟德尔遗传定律估计家庭成员缺失基因型的概率分布，然后应用混合效应logistic回归模型，并利用极大似然估计和EM算法估计基因效应值。最后，将分析结果分别与仅用先证者信息的logistic回归模型和使用主成分校正的模型进行对比。结果 该GWAS+CFD研究新发现了7个显著的单核苷酸多态性（single nucleotide polymorphisms, SNP）位点，其中4个SNPs对应已知的基因位点，分别为rs7918428(DNAJC12,OR=2.362,P=1.82×10^-9),rs6135509(MACROD2,OR=2.238,P=7.40×10^-9),rs4750496(FRMD4A,OR=2.454,P=1.12×10^-8), rs4721323(MAD1L1,OR=1.593,P=2.04×10^-8),另外3个为新发现的SNP位点：rs764009(OR=2.321,P=2.23×10^-8),rs7593443(OR=1.745,P=2.83×10^-8)和rs2170560(OR=2.603,P=3.11×10^-8)。相比于其他方法，本研究提出的方法能获得最小的基因膨胀系数（λ=1.007）。结论 本研究提出了一种新型混合效应logistic回归模型来进行带有缺失基因型观测的家族关联数据全基因组关联分析，通过加入多层随机效应有效控制了混杂因素，能显著提高检出遗传变异的能力。该方法发现了多个阿尔茨海默病的潜在致病风险位点，将有助于后续的疾病通路探索，也为疾病检测和治疗药物的研发提供了更多可能性。

目的 筛选并分析我国人均卫生费用影响因素。方法 采用LASSO方法和分位数回归模型对2018年全国31个省、直辖市、自治区的人均卫生总费用的影响因素进行分析。结果 LASSO方法从45个协变量中筛选出8个对我国各省人均卫生费用最重要的影响变量。在最终筛选出的8个变量中，孕产妇死亡率和每千人口卫生技术人员数显著促进各省人均卫生费用，且系数估计值随人均卫生费用的上升呈现先增加后减小的特征。居民消费支出的系数估计值为正，系数估计值除了70%分位点没有统计学意义，在其他四个分位点均有统计学意义。死亡率、公立医院中政府办医院数、三级医院数、门诊健康检查人数、居民人均可支配收入的系数估计值在不同分位数有正有负。结论 不同分位的人均卫生费用水平，各影响因素效应各异。对于人均卫生费用水平较低的省份，应关注居民消费支出的影响；对于人均卫生费用水平中等的省份，应关注孕产妇死亡率和每千人口卫生技术人员的影响。在增加卫生总费用投入的同时，适当控制卫生总费用的增长速度，实现卫生事业与国民经济均衡发展。

2024年7月19—21日,中国生物统计学术年会在广州市召开,来自海内外的1000余名生物统计专业人士参加了会议。中山大学国家教学名师方积乾教授作为生物统计学领域前辈,应大会特邀在开幕式上对青年学者寄语。方教授提出“三学习”——数学、现代计算技术、服务对象,“三传统”——尊师爱徒、合作共赢、同行厚道。85岁高龄的方教授激情洋溢的演讲和谆谆教诲赢得了全场参会人员的热烈掌声。应广大参会者和读者要求,我们将方教授的讲话内容整理成稿,以飨读者。

<正>第七次全国人口普查结果显示， 我国60岁及以上的老年人口总数达到2.64亿人， 占总人口的18.75%^[1], 已进入中度老龄化社会。随着我国人口老龄化的不断加剧， 老年人群疾病经济负担日益沉重， 特别是作为慢性病高发人群， 老年人的患病率是年轻人的3倍^[2]。 此外， 退行性疾病及损伤类疾病也成为老年人的多发性疾病， 不仅影响其身体健康也在一定程度给家庭带来了严重的养老负担^[3]。

目的 使用关联规则挖掘上海市≥35岁常住居民慢性病相关行为危险因素的关联性及关联强度。方法 利用上海市2017年慢性病及其危险因素监测数据，描述≥35岁常住居民现在吸烟、危险或有害饮酒、蔬菜和水果摄入不足、红肉摄入过量、身体活动不足和每天睡眠时长不合格6项行为危险因素共存情况，采用R 4.1.2软件arules包Apriori算法对6项行为危险因素在不同性别和年龄组人群的关联情况进行关联规则分析。结果 19578名居民中至少存在2项行为危险因素的占比为39.4%。支持度超过5%、置信度超过50%且提升度最高的关联规则为：男性35～59岁组{危险或有害饮酒}→{现在吸烟}(C=82.38%)、男性≥60岁组{危险或有害饮酒}→{现在吸烟}(C=63.78%)、女性≥60岁组{蔬菜和水果摄入不足}→{睡眠时长不合格}(C=51.35%)。结论 上海市≥35岁常住居民慢性病行为危险因素共存现象常见且相互间存在关联，不同特征人群的共存关联特点不一，提示在今后的行为危险因素干预项目中，需提高行为危险因素综合干预的意识，根据不同特征人群特点实施有针对性的多项行为危险因素综合干预措施。

目的 探讨江苏省成年人体质指数纵向变化轨迹及分类。方法 基于中国营养健康调查数据，使用线性混合模型树探究江苏省18～65岁人群体质指数的变化轨迹及其分类情况。结果 线性混合模型树生成了13个节点，树的深度为6,分类节点为基线BMI、平均摄入卡路里、基线年龄。结论 线性混合模型树可以识别体质指数的变化轨迹，拓展了纵向数据动态变化的研究方法。

目的 探讨基于积极心理学理论的心理危机干预对急性脑卒中患者的影响。方法 选取2018年5月-2020年5月我院113例急性脑卒中患者进行前瞻性研究，按照建档顺序分为对照组(n=56)和实验组(n=57),分别给予常规心理护理、常规心理护理+基于积极心理学理论的心理危机干预，连续干预4周。观察并对比护理前后患者的心理状态(type D personality scale-14, DS-14)、创伤后成长、创伤后应激障碍(PTSD self-rating scale, PTSD-SS)评分、神经功能(stroke neurological deficit score, CSS)、日常生活能力(activities of daily living, ADL)情况。结果 干预4周后，实验组社交抑制、消极情感等维度DS-14评分低于对照组(P<0.05);实验组主观评定、回避症状、社会功能、反复重现体验和警觉性增高各维度PTSD-SS评分低于对照组(P<0.05);实验组新的可能性、与他人关系、精神改变、个人力量和欣赏生活等维度PTGI评分高于对照组(P<0.05);干预4周后，两组神经功能CSS评分低于干预前，日常生活能力ADL评分高于干预前(P<0.05),但两组CSS评分、 ADL评分相比差异无统计学意义(P>0.05)。结论 基于积极心理学理论的心理危机干预能明显改善急性脑卒中患者心理状态，减轻创伤后应激障碍，促进创伤后成长。

目的 研究某三甲综合医院冠心病患者住院费用的构成及影响因素，分析不同病例组合，为冠心病患者住院费用提供参考标准。方法 以某三甲医院2016年1月1日至2020年12月31日住院病案首页中主要诊断为冠心病(I25.1)的10129例住院患者为研究对象，采用单因素分析、 logistic回归分析获得住院费用的影响因素，采用基于CHAID算法的决策树模型进行病例组合分析。结果 冠心病患者平均住院费用为14741.49元，中位数为8351.34元，其中诊断费、药品费和材料费构成较高。住院费用的影响因素包括医疗付费方式、入院时情况、实际住院天数、是否转科、是否抢救、其他诊断个数、手术个数等。将具有统计学意义的影响因素纳入决策树模型，共形成22个病例组合。结论 医疗机构应实施科学、规范的临床路径，加强对高值耗材使用的精细化管理，充分利用医疗资源，提高服务效率，缩短住院时间，减轻患者经济负担。通过决策树模型获得的冠心病DRGs分组合理，住院费用标准及病种权重具有一定参考价值。

目的 在高体重指数(body mass index, BMI)率快速增长的情况下，分析1990-2019年中国归因于高BMI的2型糖尿病疾病负担。方法 基于2019年全球疾病负担研究数据，分性别、年龄组分析1990-2019年中国因高BMI造成的2型糖尿病疾病负担，并利用Joinpoint回归模型分析标化后的归因死亡率和归因伤残调整寿命年(disability-adjusted life year, DALY)率的变化趋势。结果 1990-2019年间，我国2型糖尿病患病率从2928.78/10万增至6328.79/10万，归因高BMI的2型糖尿病死亡人数从1.05万增至4.75万，标化死亡率从1.25/10万增至2.39/10万；归因DALY从77.18万人年增至373.76万人年，标化后的归因DALY率从80.21/10万增至181.54/10万，早死损失寿命年(years of life lost, YLL)和伤残损失寿命年(year lived with disability, YLD)及其标化率也呈大幅增长。1990-2019年间，我国归因于高BMI的2型糖尿病标化死亡率和标化DALY率的年平均变化百分比分别为2.28%和2.81%,均有统计学意义(P<0.05),且男性均高于女性。男性标化后的归因DALY率和归因死亡率分别在2010年和2014年超过女性。年龄分层结果显示，50岁以上人群高BMI导致的2型糖尿病疾病负担更加严重，15～49岁人群归因于高BMI的YLD率增幅最大，达到323.99%。结论 1990-2019年间，中国归因于高BMI的2型糖尿病疾病负担显著增加。需加大防控力度，有效管理人群体重指数，并对高危人群采取重点干预，以减轻2型糖尿病疾病负担。

目的 探究老年患者胃肠手术后手术部位感染(surgical site infection, SSI)发生水平及相关高危因素，为SSI防控提供理论依据。方法 对沧州市中心医院1 704例老年胃肠手术患者的病历资料进行回顾性分析。SSI包括手术切口、深部器官或腔隙的感染。对手术后发生SSI的危险因素进行单因素和多因素分析。结果 1 704例老年胃肠手术患者SSI发生率为11.5%(196例)。按手术类型分层后，胆囊切除术、结肠、胃、直肠、小肠手术发生SSI的老年患者中术中输血、手术时间长及手术切口Ⅲ/Ⅳ类的患者比例较高(P<0.05)。多因素分析结果示，常见胃肠手术患者SSI的主要危险因素依次是手术切口Ⅲ/Ⅳ类、术中输血、类固醇类药物、抗糖尿病药物、急诊手术、ASA分级≥三级及年龄(P<0.001)。结论 胃肠手术老年患者SSI的发生率相对较低，影响因素较多。预防SSI发生应针对5种以上手术主要危险因素进行重点干预。

目的 分析2016-2021年广州市南沙区道路交通伤害的疾病负担现状及变化趋势，为防控政策提供理论依据。方法 以2016-2021年广州市南沙区道路交通伤害死亡人群为研究对象，通过计算死亡率及标化死亡率、早死所致寿命损失年(years of life lost, YLL)及标化早死所致寿命损失年、早死所致寿命损失年率及标化早死所致寿命损失年率等指标对道路交通伤害所致疾病负担进行描述，并分析其年龄分布及性别差异，并使用年度变化百分比分析其变化趋势。结果 2016-2021年共有死亡例数329例；粗死亡率为12.10/10万，标化死亡率12.20/10万；YLL为11950.76人年，标化YLL为12861.22人年，YLL率为4.39‰,标化YLL率为4.73‰。除2021年外，2016-2020年不同性别道路交通伤害所致死亡率及标化死亡率差异均有统计学意义。男性死亡人数、粗死亡率及标化死亡率有逐年下降的趋势，差异有统计学意义。对不同性别的年龄分布进行检验，差异有统计学意义。总体中位死亡年龄为50(31～62)岁，男性中位死亡年龄为45(29～58)岁，女性中位死亡年龄为57.5(43～69)岁。不同性别发生道路交通伤害时使用交通工具分布差异有统计学意义，其中男性使用摩托车发生道路交通伤害例数最多。不同性别疾病负担不同；男性疾病负担有时间变化趋势，女性疾病负担稳定无明显趋势。结论 不同性别道路交通伤害疾病负担不同；男性疾病负担有随时间下降的变化趋势。男性发生道路交通伤害死亡年龄更年轻；使用摩托车是发生道路交通伤害的危险因素。

目的 对现有样本代表性评估方法进行全面比较和探讨，为临床研究样本代表性评估方法选择提供参考。方法 结合国内肺癌患者特征的分布以及国内临床研究样本筛选的实际情况，模拟肺癌患者目标人群，抽取不同样本量和不同偏离程度的样本，使用现有样本代表性评估方法计算样本代表性，同时计算疗效估计偏差（bias）,通过建立各方法代表性测量值与bias之间的相关性模型，分析各方法评估代表性的准确性和稳定性。结果 整体结构差异率（rate of overall struction variation, RV）RV1和RV2及基于倾向评分的C统计量、基尼集中比求和（sum Gini concentration ratio, SGCR）及K-S距离（kolmogorov-smirnov distance, KSD）均能较好地测量不同样本的偏离程度。在不同样本量下，RV2和RV1与bias相关模型的R²值均大于0.90,C统计量、SGCR及K-S距离的R²大于0.80。结论 因考虑了特征权重，整体结构差异率更为准确、稳定，尤其是RV2能更好地测量不同偏离程度样本的代表性、准确反映估计偏差；在难以获得特征重要性信息时，SGCR及利用倾向评分的方法中的C统计量和K-S距离测量代表性的可靠性也可以接受。

目的 初步编制国家卫生城市主观建成环境满意度评价量表，并进行信度和效度评价。方法 基于政策梳理、专家咨询确定量表的初始条目，通过居民个人访谈结合统计分析方法确定量表的最终条目，利用探索性因子分析主成分分析确定量表的维度，最后使用克朗巴赫系数和分半系数以及验证性因子分析对量表进行信度和效度评价。结果 本研究建立了包括20个条目的4维度量表，分别代表城市生活方式、城市基础功能、城市环境卫生和城市治理四个方面的评价。问卷的克朗巴赫系数为0.910,分半系数为0.851;验证性因子分析结果显示各检验量均处于标准范围。结论 本研究初步编制的国家卫生城市建成环境满意度评价量表信度和效度良好，应开展量表在人群中的实证研究，进一步验证量表的实际应用价值。