目的 考察大学生自杀意念的动态发展轨迹及其相关因素.方法 在2019-2021年间,采用方便取样的方法,每年10月对广东省某高校1438名学生进行3次纵向追踪.采用自杀意念自评量表(SIOSS)测量自杀意念水平,使用潜变量混合增长模型(LGMM)分析大学生自杀意念发展轨迹,建立多元logistic回归模型,分析自杀意念各潜在类别变化的相关因素.结果 2019-2021年间,大学生自杀意念总分分别为3.79±2.87、3.71±2.90、3.32±2.87,随时间呈现一定的下降趋势;大学生自杀意念的发展轨迹可划分为3个类别:C1自杀意念加重组(11.8％)、C2持续低自杀意念组(79.6％)、C3自杀意念缓解组(8.6％);相对于C1组,童年生活一般(OR=2.52,95％CI:1.31～4.86)或童年生活幸福(OR=3.95,95％CI:2.02～7.74)、无自杀未遂史(OR=9.03,95％CI:3.49～23.36)、外向(OR=1.03,95％CI:1.01～1.05)的大学生更容易进入C2组,而生源地为城市(OR=0.54,95％CI:0.31～0.93)、神经质得分高(OR=0.86,95％CI:0.84～0.88)的大学生更不容易进入C2组;童年生活幸福(OR=3.29,95％CI:1.40～7.74)的大学生更容易进入C3组,而生源地为城市(OR=0.37,95％CI:0.18～0.78)的大学生更不容易进入C3组.结论 大学生自杀意念的发展轨迹可划分为3个类别:自杀意念加重组、持续低自杀意念组、自杀意念缓解组,自杀意念的发展轨迹与生源地、童年生活、自杀史、人格特质等因素相关.

目的:探讨腹腔镜结直肠癌根治术患者术后1年内创伤后应激障碍(posttraumatic stress disorder,PTSD)的纵向变化轨迹,分析心理维度、社会支持维度及炎症因子对轨迹类别的预测作用.方法:本研究为前瞻性观察性纵向研究.选取2021年2月至2022年8月我院结直肠癌行腹腔镜根治术的患者为研究对象,于术后6周(T1)、术后3个月(T3)、术后6个月(T4)、术后9个月及术后12个月(T5)采用创伤后应激障碍平民版量表调查患者5个时间点的PTSD水平,采用潜类别增长模型(latent class growth model,LCGM)识别轨迹类别,采用Logistics回归分析识别轨迹类别的预测因素.结果:采用LCGM方法将193例腹腔镜结直肠癌根治术患者术后1年内PTSD的变化轨迹分为3个潜在类别即PTSD低水平组(占41％)、PTSD下降组(占33％)、PTSD进展组(占26％);选择PTSD低水平组作为参考类别,采用逐步进入法筛选变量,Logistic回归分析结果显示,高中及以上文化程度、下尿路不适及炎症因子[肿瘤坏死因子-α(tumor necrosis factor-α,TNF-α)、白细胞介素-6(interleukin-6,IL-6)]水平越高的患者更易进入PTSD下降组,炎症因子(IL-6、TNF)、恐惧疾病进展及侵入性认知水平越高的患者更易进入PTSD进展组.结论:腹腔镜结直肠癌根治术患者术后1年内的PTSD存在异质性,有3种不同的变化轨迹,尤其应重视初始分值较低的PTSD进展组.临床医护人员应关注患者术后早期的身体不适症状,同时重视心理因素如恐惧疾病进展、侵入性认知及炎症因子对PTSD轨迹的预测作用,制定针对性的措施预防术后PTSD慢性进展.

目的 探讨头颈部恶性肿瘤放疗患者护理未满足需求发展轨迹及其相关因素.方法 选取2020年1月至2022年7月西安交通大学医学院第一附属医院接受放疗的头颈部恶性肿瘤患者230例作为研究对象,采用癌症患者未满足需求简明量表(SF-SUNS)评价放疗前(T1)、放疗结束时(T2)、放疗结束后3个月(T3)、6个月(T4)患者护理需求未满足情况.基于潜类别增长模型分析头颈部恶性肿瘤放疗患者护理未满足需求发展轨迹,采用多分类Logistic回归模型筛选其影响因素.结果 T1～T4时间点患者SF-SUNS评分分别为(70.23±5.20)分、(78.10 士6.15)分、(72.65±6.84)分及(72.20±6.21)分,各时间点SF-SUNS评分比较差异均有统计学意义(F=66.330,P＜0.05).潜类别增长模型分析结果显示,当潜变量类别数为3时,Akaike信息标准(AIC)、经过校正Bayesian信息标准(aBIC)值最小,3类别模型为最佳模型,按未满足需求量表评分＞60分为高需求未满足为参照,1类别在T1时SF-SUNS评分高于此水平,然后T4时间点降低,命名为高需求降低组,共69例(31.36％);2类别起始SF-SUNS评分较高后缓慢降低,命名为高需求稳定组,共77例(35.00％);起始SF-SUNS评分较低后升高,命名为低需求上升组,共74例(33.64％).3组年龄、受教育年限、家庭月收入、社会支持量表(SSRS)评分比较差异有统计学意义(P＜0.05);多分类Logistic回归分析显示:SSRS评分、受教育年限为头颈部恶性肿瘤放疗患者护理未满足需求上升的危险因素,年龄、家庭月收入＞5 000元为未满足需求上升的保护因素(P＜0.05).结论 头颈部恶性肿瘤放疗患者护理未满足需求发展轨迹存在个体差异,应重点关注患者年龄、受教育年限、家庭经济情况及社会支持度,针对头颈部恶性肿瘤放疗患者需求未满足发展轨迹类型及影响因素制定个性化干预方案.

目的探讨冠状动脉旁路移植术(coronary artery bypass grafting,CABG)患者术后疲劳发展轨迹及其影响因素.方法采用便利抽样法,选取2022年10月-2023年8月在我院心脏外科行CABG的145例患者为研究,随访过程中剔除18例,最终共纳入124例患者.分别在出院前1 d及出院后1、2、4个月采用一般资料调查表、中文版多维疲劳量表(multidimensional fatigue inventory-20,MFI-20)及Barthel指数评定量表对患者术后疲劳水平进行追踪调查,利用潜类别增长模型(latent class growth model,LCGM)识别潜在的术后疲劳轨迹类别,采用多分类logistic回归分析患者术后疲劳轨迹类别的影响因素.结果患者术后MFI-20得分随时间呈下降趋势(F=310.039,P＜0.001),得分从出院前1 d(58.44±8.65)分降至出院后4个月(36.92±9.89)分;最终识别出3种术后疲劳轨迹类别,根据其特征分别命名为中至轻度明显下降组、中度缓慢下降组、重至中度缓慢下降组,分别占比53.2％、32.3％、14.5％.单因素分析显示:有无心肌梗死病史、有无脑卒中病史、NYHA心功能分级、术后ICU停留天数、术后左室射血分数(LVEF)、红细胞计数、血红蛋白、总胆固醇、低密度脂蛋白、日常生活活动能力与术后疲劳轨迹类别有关(均P＜0.05);多因素logistic回归分析显示:术后ICU停留天数、有无心肌梗死病史、NYHA心功能分级、有无脑卒中病史、日常生活活动能力是CABG患者术后疲劳轨迹类别的主要影响因素(均P＜0.05)o结论CABG患者术后疲劳发展轨迹存在异质性,医护人员应根据不同变化轨迹制定针对性护理措施,积极改善可控因素,缓解患者术后疲劳程度.

目的:应用潜分类增长模型（latent class growth model，LCGM）纵向分析孕期血压数据，探索妊娠期间血压潜分类轨迹与未来发生子痫前期（preeclampsia，PE）风险之间的联系。方法:该研究为前瞻性队列研究，选取2016年11月1日至2018年5月30日在天津市19家社区医院招募的孕早期女性受试者，采集孕28周前5个阶段（即12、16、20、24和28周）的孕检血压数据纳入分析，同时采集孕妇基本信息问卷、体格检查、实验室检查等数据。应用LCGM建立收缩压、舒张压的潜分类轨迹模型；分别以收缩压轨迹、舒张压轨迹作为预测变量，以妊娠28周后至分娩时是否发生PE作为结局变量，进行PE的预测研究。结果:共纳入5 560例已生产的单胎妊娠孕妇，其中随访期间诊断PE 128例，妊娠期高血压106例。单因素及多因素logistic回归分析显示基线收缩压与舒张压水平越高，PE发生的风险越高。采用LCGM拟合血压动态变化，确定了4条收缩压潜分类轨迹和4条舒张压潜分类轨迹。收缩压潜分类轨迹对应的血压水平越高，PE发生的风险越高；在校正了初检体重指数、是否初产妇、白细胞计数、血红蛋白水平、血小板计数和谷丙转氨酶水平等混杂因素后，收缩压潜分类轨迹_4（SBP_traj_4）的 OR值为4.023（95% CI：2.368~6.835， P<0.001），收缩压潜分类轨迹_3（SBP_traj_3）的 OR值为1.854（95% CI：1.223~2.811， P=0.004）；舒张压潜分类轨迹_4（DBP_traj_4）的 OR值为2.527（95% CI：1.534~4.162， P<0.001），舒张压潜分类轨迹_2（DBP_traj_2）的 OR值为2.238（95% CI：1.328~3.772， P=0.002）。受试者工作特征曲线显示基线血压水平联合血压轨迹对PE具有较好的预测能力。 结论:应用LCGM潜分类增长模型构建孕12 ~28周的血压潜分类轨迹可对孕28周后PE进行早期预测分析，提示孕期血压潜分类轨迹是一种可用于孕妇PE危险分层的新方法。

目的:应用贝叶斯网状Meta分析预测头孢菌素类抗菌药物联合甲硝唑经验性治疗腹腔感染的疗效和安全性。方法:以( "intraabdominal infections" [MeSH Terms])AND( "Cephalosporins*" [MeSH Terms])AND( "randomized controlled trial" [MeSH Terms])、腹腔感染、继发性腹膜炎、腹腔脓肿、头孢为检索词，检索PubMed、Embase、the Cochrane Library、中国知网、万方数据库、维普数据库。检索时间为1990年1月至2018年5月。纳入评价头孢菌素类抗菌药物联合甲硝唑对比其他类型抗菌药物治疗腹腔感染疗效和安全性的随机对照试验(RCT)。试验组患者给予头孢菌素类抗菌药物联合甲硝唑治疗腹腔感染，对照组患者给予其他类型抗菌药物治疗腹腔感染。结局指标：临床治愈率、微生物清除率和药物严重不良反应发生率。使用R 3.6.2软件随机贝叶斯模型进行网状Meta分析。采用马尔可夫链蒙特卡罗法直接评估和间接预测。采用轨迹图法、密度图法、杠杆图法对模型收敛性和稳定性进行评价。因各干预措施间没有闭环，不需要评估研究的一致性。结果:(1)文献检索结果：最终纳入符合标准的相关文献18篇。18篇文献均为RCT，文献累积样本量为6 792例，其中试验组3 402例，对照组3 390例。(2)贝叶斯网状Meta分析结果：①第3代头孢菌素＋甲硝唑、碳青霉烯类的临床治愈率显著低于第4代头孢菌素＋甲硝唑，差异均有统计学意义[优势比( OR)＝0.46，0.61，95%可信区间( CI)为0.26～0.81，0.38～0.97， P<0.05]。第5代头孢菌素＋甲硝唑的临床治愈率与碳青霉烯类比较，差异无统计学意义( OR＝1.03，95% CI为0.59～1.80， P>0.05)。②第5代头孢菌素＋甲硝唑、碳青霉烯类的微生物清除率显著低于第4代头孢菌素＋甲硝唑，差异均有统计学意义( OR＝0.84，0.41，95% CI为0.73～0.98，0.23～0.74， P<0.05)。第5代头孢菌素＋甲硝唑的微生物清除率与碳青霉烯类比较，差异无统计学意义( OR＝0.76，95% CI为0.27～1.80， P>0.05)。③第3代头孢菌素＋甲硝唑、第4代头孢菌素＋甲硝唑、第5代头孢菌素＋甲硝唑、碳青霉烯类、喹诺酮＋甲硝唑、替加环素的药物严重不良反应发生率均显著低于喹诺酮类，差异均有统计学意义( OR＝0.13，0.13，0.14，0.13，0.15，0.13，95% CI为0.03～0.50，0.02～0.98，0.02～0.75，0.02～0.59，0.02～0.78，0.02～0.57， P<0.05)。第3代头孢菌素＋甲硝唑、第4代头孢菌素＋甲硝唑、第5代头孢菌素＋甲硝唑的药物严重不良反应发生率与碳青霉烯类比较，差异均无统计学意义( OR＝0.96，1.00，1.10，95% CI为0.52～1.60，0.31～3.50，0.49～2.30， P>0.05)。(3)疗效与安全性等级排名：①临床治愈率排名：各类用药方案的临床治愈率从高到低排名顺序为喹诺酮类＋甲硝唑、第4代头孢菌素＋甲硝唑、合成青霉素类、第2代头孢菌素＋甲硝唑、第5代头孢菌素＋甲硝唑、碳青霉烯类、第3代头孢菌素＋甲硝唑、替加环素、喹诺酮类。各排名相应的概率为51.73%、35.72%、22.57%、31.37%、24.98%、32.82%、34.69%、29.05%、72.36%。②微生物清除率排名：各类用药方案的微生物清除率从高到低排名顺序为喹诺酮类＋甲硝唑、第4代头孢菌素＋甲硝唑、第2代头孢菌素＋甲硝唑、合成青霉素类、第5代头孢菌素＋甲硝唑、碳青霉烯类、第3代头孢菌素＋甲硝唑、替加环素、喹诺酮类。各排名相应的概率为：89.62%、77.01%、38.60%、20.94%、26.26%、26.39%、22.22%、20.19%、62.55%。③药物严重不良反应发生率排名：各类用药方案的药物严重不良反应发生率从高到低排名顺序为喹诺酮类、喹诺酮类＋甲硝唑、第5代头孢菌素＋甲硝唑、碳青霉烯类、第3代头孢菌素＋甲硝唑、替加环素、第4代头孢菌素＋甲硝唑。各排名相应的概率为：96.21%、30.46%、21.09%、25.27%、27.26%、19.45%、31.69%。 结论:中、低危腹腔感染时，推荐经验性使用头孢菌素＋甲硝唑替代经验性使用碳青霉烯类。

目的:探讨收缩压轨迹与恶性肿瘤发生的关系。方法:收集54 888名开滦(集团)有限责任公司职工2006—2007年、2008—2009年、2010—2011年3次健康体检的相关资料，及其恶性肿瘤新发情况。利用研究对象3次健康体检的血压测量值，进行收缩压轨迹分组。采用寿命表法计算恶性肿瘤发生率，并进行恶性肿瘤发病影响因素的多因素Cox比例风险模型回归分析。结果:根据收缩压轨迹将54 888人分为5组，其中低-稳定组14 326人，正常-稳定组25 630人，中-高组5 390人，高-中组6 438人，高-稳定组3 104人。在(4.95±0.53)年的随访期内，共新发恶性肿瘤1 070例。低-稳定组、正常-稳定组、中-高组、高-中组和高-稳定组的恶性肿瘤发生率分别为1.3%(177/14 326)、2.2%(491/25 360)、3.1%(147/5390)、2.7%(156/6 438)和3.8%(99/3 104)，差异有统计学意义( P<0.001)。多因素Cox回归分析显示，校正性别、年龄、吸烟、饮酒、体育锻炼、体质指数、空腹血糖、总胆固醇、服用降压药、降糖药、降脂药后，收缩压轨迹与恶性肿瘤发病有关，与低-稳定组相比，正常-稳定组、中-高组、高-中组、高-稳定组的 HR分别为1.413、1.731、1.557和1.907(均 P<0.001)。 结论:高收缩压轨迹是恶性肿瘤发生的危险因素。

目的:通过降低斑马鱼接触蛋白相关蛋白基因2（contact protein-related protein gene 2, cntnap2）表达，建立注意缺陷多动障碍（attention-deficit/hyperactivity disorder，ADHD）遗传模型。方法:以注射剪接抑制型吗啉代寡聚核糖核酸降低cntnap2表达的斑马鱼作为实验组（ n=48），注射随机寡核苷酸片段的斑马鱼作为对照组（ n=48)，通过原位杂交观察斑马鱼神经递质系统相关基因th、dbh、gad1b、glyt2、vglut2表达的变化，使用实时荧光定量PCR检测神经递质系统及神经元发育相关基因表达量变化。使用斑马鱼行为轨迹跟踪系统记录（Noldus行为仪）检测cntnap2表达降低后多动/冲动行为及托莫西汀对多动/冲动行为的改善作用。 结果:原位杂交实验结果显示实验组较对照组神经递质系统相关基因dbh表达减少，gad1b表达增加；th、glyt2、vglut2表达未见明显差异。实时荧光定量PCR显示实验组较对照组神经递质系统相关基因dbh表达降低（ t=5.772， P=0.005），snap25、dat、gad1b表达增高（ t=14.310、22.210、22.090,均 P<0.01）；神经元发育相关基因shha、epha4b、netrin1b、noi及神经递质系统相关基因drd4、th、vglut2、glyt2表达差异均无统计学意义。6 dpf时药物未处理实验组较药物未处理对照组幼鱼游动总距离长 ［143.12(98.56, 212.22) cm 与99.03(80.54, 133.96) cm， Z=-3.547， P<0.01］，平均游动速度快［0.050(0.041, 0.067) cm/s与0.033(0.025, 0.042) cm/s, Z=-5.495， P<0.01］；经托莫西汀处理后，实验组药物处理较药物未处理幼鱼游动总距离减少［106.59(95.28, 120.10) cm与143.12(98.56, 212.22) cm, Z=-3.591， P<0.01］，平均游动速度减慢［0.031(0.028, 0.037) cm/s与0.050(0.041, 0.067) cm/s， Z=-7.035， P<0.01］。 结论:cntnap2表达降低斑马鱼可通过改变dbh、snap25、dat、gad1b表达量而产生多动/冲动表型，托莫西汀可有效缓解多动/冲动表型，cntnap2表达降低斑马鱼可作为ADHD遗传模型。

随着临床研究的发展，纵向数据的应用越来越受到重视。实践中常见多次暴露测量及单次结局测量的纵向数据，临床医生在分析中往往难以充分利用其信息。轨迹模型是一种用于处理这类纵向数据的灵活工具，可根据个体发展轨迹识别总体中的异质亚群，常见方法包括增长混合模型和潜类别增长模型。本文从实践的角度出发，以R语言为例，并结合示例数据，详细介绍轨迹模型的操作步骤、软件实现和分析结果的展示。

目的:探讨不同喂养不耐受变化轨迹外科危重症患者预后及影响因素。方法:采用回顾性队列研究方法。收集2018年3月至2019年7月中国重症营养临床研究小组-NEED研究中69家医学中心收治的354例外科危重症患者的临床资料；男247例，女107例；年龄为58（46，68）岁。根据354例患者喂养不耐受变化轨迹模型，分为无喂养不耐受、喂养不耐受下降、喂养不耐受持续分别为164、49、141例。观察指标：（1）不同喂养不耐受变化轨迹患者一般情况。（2）不同喂养不耐受变化轨迹患者营养治疗情况。（3）不同喂养不耐受变化轨迹患者生存情况。（4）外科危重症患者预后影响因素分析。正态分布的计量资料以 x±s表示，组间比较采用单因素方差分析；偏态分布的计量资料以 M（ Q1， Q3）表示，组间比较采用Kruskal-Wallis秩和检验。计数资料以绝对数或百分比表示，组间比较采用 χ2检验。等级资料比较采用Kruskal-Wallis秩和检验。两两比较采用Bonferroni校正。根据Stata17.0统计软件中Traj插件构建组基轨迹模型，最佳轨迹模型采用贝叶斯信息准则、平均后验概率参数评价。采用Kaplan-Meier法绘制生存曲线并计算生存率，Log-Rank检验进行生存分析。采用COX比例风险回归模型进行单因素和多因素分析。 结果:（1）不同喂养不耐受变化轨迹患者一般情况。354例外科危重症患者中，行肠内营养治疗257例、行肠内营养+肠外营养治疗97例，急性生理与慢性健康（APACHE Ⅱ）评分为17（13，21）分，序贯性器官功能衰竭（SOFA）评分为6（5，8）分，改良危重症营养风险（mNUTRIC）评分为4（2，5）分，合并症数量为2（1，3）个，急性胃肠损伤分级为Ⅰ级、Ⅱ级、Ⅲ级分别为293、55、6例，行机械通气、连续肾替代治疗、血管活性药物使用分别为224、17、61例。354例患者喂养不耐受发生率呈先上升后下降趋势，第3天发生率达峰值为25.42%（90/354），7 d内发生率为53.67%（190/354）。354例外科危重症患者中，无喂养不耐受、喂养不耐受下降、喂养不耐受持续患者APACHE Ⅱ评分，mNUTRIC评分，合并症数量，急性胃肠损伤分级，机械通气分别为16（12，20）分、17（14，25）分、18（13，22）分，3（2，5）分、4（3，6）分、4（3，5）分，2（1，2）个、2（2，3）个、2（2，3）个，Ⅰ级（152例、27例、114例）、Ⅱ~Ⅲ级（12例、22例、27例），95、39、90例，三者上述指标比较，差异均有统计学意义（ H=6.14、13.11、28.05， χ2=37.96、7.65， P<0.05）；进一步分析，与无喂养不耐受患者比较，喂养不耐受下降、喂养不耐受持续患者均有更高的合并症数量和急性胃肠损伤分级（ Z=60.32、54.69， χ2=39.72、9.52， P<0.05）；与无喂养不耐受患者比较，喂养不耐受下降患者有更高的mNUTRIC评分和机械通气比例（ Z=53.41， χ2=7.59， P<0.05）。（2）不同喂养不耐受变化轨迹患者营养治疗情况。无喂养不耐受患者促胃动力药物使用，幽门后喂养分别为36例，13例，喂养不耐受下降患者上述指标分别为25例，10例，喂养不耐受持续患者上述指标分别为46例，19例，三者上述指标比较，差异均有统计学意义（ χ2=15.76，6.20， P<0.05）。进一步分析，与无喂养不耐受患者比较，喂养不耐受下降患者有更高的促胃动力药物使用和幽门后喂养比例（ χ2=15.60，6.10， P<0.05）。（3）不同喂养不耐受变化轨迹患者生存情况。无喂养不耐受、喂养不耐受下降、喂养不耐受持续患者28 d总生存率分别为96.96%、95.92%、87.94%，三者比较，差异有统计学意义（ χ2 =10.39， P<0.05）；进一步分析，无喂养不耐受与喂养不耐受持续患者比较，差异有统计学意义（ χ2 =9.19， P<0.05）。（4）外科危重症患者预后影响因素分析。多因素分析结果显示：喂养不耐受持续是外科危重症患者28 d死亡的独立危险因素（风险比=3.92，95%可信区间为1.43~10.79， P<0.05）。 结论:与无喂养不耐受患者比较，喂养不耐受持续的外科危重症患者28 d总生存率下降，喂养不耐受持续是外科危重症患者28 d死亡的独立危险因素。