该研究针对血清尿酸(UA)水平与肌肉减少症之间的关系这一关键知识空白,特别是在中国人群中,此类关联尚未得到全面探讨。
该研究针对血清尿酸(UA)水平与肌肉减少症之间的关系这一关键知识空白,特别是在中国人群中,此类关联尚未得到全面探讨。
既往研究关于血清尿酸(UA)水平与肌少症之间关系的发现尚不一致。然而,目前尚无研究全面探讨这一关系在中国更广泛人群中的表现。本研究旨在探讨45岁及以上中国成年人中血清尿酸水平与肌少症之间的关系,并重点关注与年龄相关的差异。本研究纳入了2015年中国健康与养老追踪调查(CHARLS)中的10,938名参与者。采用加权逻辑回归和加权线性回归模型评估肌少症(包括其各组分)与血清尿酸水平之间的关联。在按年龄组别对参与者进行分类后,通过亚组分析进一步深入考察了与年龄相关的变化。参与者根据尿酸水平被分为四个四分位组。经校正的分析结果显示,在65岁及以上的个体中,较高的血清尿酸水平仅在此人群中与肌少症呈负相关。加权线性回归分析结果表明,血清尿酸水平与握力(HGS)和骨骼肌指数(SMI)之间均存在具有统计学意义的正相关关系。此外,Q4组(≥5.70 mg/dL)在所有年龄组中均保持了这种正相关性。结果表明,在45岁及以上的中国人中,较高的UA水平与较高的SMI和HGS显著相关。升高的血尿酸水平可能对肌少症的发生具有潜在的保护作用,尤其是在65岁及以上的个体中。
肌少症是一种影响骨骼肌的退行性疾病,导致肌肉质量逐渐减少。该病与跌倒、功能衰退、虚弱和死亡等不良后果相关1,2,3,已成为一个重大的全球公共卫生问题。过去,肌少症常被视为一种老年疾病,但近期研究表明,由于长期不良的生活习惯和机体老化,部分中年人群的肌肉质量已显著下降,其肌少症发病率不容忽视4。全球肌少症的患病率因诊断标准和目标人群不同而存在显著差异,介于10%至27%之间5。既往研究显示,老年亚洲人群中肌少症的患病率在6.8%至25.7%之间波动6,7,8,9。中国肌少症的总体患病率约为16.4%-18%,且患病率随年龄增长呈上升趋势10。肌少症的发生不仅影响个体和家庭,还加重了社会医疗系统的负担。随着全球人口老龄化加剧,识别可干预的肌少症发生与进展危险因素,对于预防或延缓其发展至关重要。
肌少症是一种复杂的疾病,与多种危险因素相关,例如体力活动减少、热量摄入降低、慢性炎症、神经肌肉接头退化以及氧化应激11,12。作为肌少症的重要致病因素之一,减轻氧化应激所造成的损伤可能是预防肌少症的有效途径。尿酸(UA)是一种重要的内源性抗氧化剂,能够清除活性氧(ROS),从而缓解氧化应激的影响。然而,尽管尿酸具有显著的抗氧化特性,但也表现出矛盾的促氧化活性13,14,因此尿酸对肌少症的影响仍存在争议。现有研究表明,尿酸(UA)可能通过以下机制影响肌少症:一方面,尿酸的抗氧化特性可减少活性氧(ROS)对肌细胞的损伤,维持肌细胞的正常结构与功能,并促进肌肉蛋白质合成,从而对肌肉质量和力量产生积极影响15;另一方面,尿酸的促氧化活性可能诱发炎症反应,干扰细胞内信号通路,并抑制肌肉蛋白质合成,加速肌肉蛋白质降解,导致肌肉质量和力量下降16。既往研究提示,在中国人群中,尿酸水平与肌肉质量呈正相关17,而在日本和美国人群中,尿酸水平与肌肉力量呈正相关18,19。相反,也有研究指出,尿酸水平升高与肌肉质量减少及肌肉力量下降相关20,21。
关于血清尿酸与肌少症之间关联的现有研究存在若干局限性,包括样本量较小、数据质量较差以及尿酸检测精度不足。这些问题使得在不同人群中就尿酸水平与肌少症之间的关系得出可靠结论变得困难。此外,探讨血清尿酸与肌少症在不同年龄阶段差异的研究十分有限。本研究旨在利用中国健康与养老追踪调查(CHARLS)的全国代表性数据进行横断面分析,以探讨45岁及以上中国成年人中尿酸与肌少症之间的相关性,并同时探索潜在的年龄相关差异。
中国健康与养老追踪调查(CHARLS)已获得北京大学机构审查委员会(IRB)的批准,批准编号为 IRB00001052-11014 和 IRB00001052-11015。所有参与者在参加研究前均签署了知情同意书,符合伦理标准和《赫尔辛基宣言》的要求。这确保了所使用的 CHARLS 数据遵循伦理研究准则,包括保密性和知情参与的要求。
研究人群
本研究的数据来源于北京大学国家发展研究院于2015年开展的中国健康与养老追踪调查(CHARLS)。CHARLS始于2011年,是一项针对中国45岁及以上成年人的纵向研究。该调查通过使用标准化问卷进行结构化访谈,系统收集了与老龄化相关的社会经济因素及其影响的全面数据集。我们采用多阶段分层概率比例抽样方法,从17,708名参与者中收集数据。这些参与者居住在中国28个省份的150个区县、450个村庄中的10,257户家庭中22。数据包含个体加权变量,以确保调查样本具有全国代表性。有关CHARLS研究设计的其他信息已有前期记录6,更多关于CHARLS研究的内容可从其官方网站获取(https://charls.pku.edu.cn/)。
纳入标准如下:(a) 2015年时年龄为45岁或以上的参与者,(b) 有关于其肌少症状态的可用数据的参与者,以及 (c) 接受过尿酸(UA)血液检测的参与者。考虑到这些疾病对血清尿酸水平的影响,患有肾脏疾病、肿瘤和痛风的参与者被排除在外。最终,分析共纳入10,938名参与者。图1展示了排除过程的详细情况。本研究遵循加强流行病学中观察性研究报告规范(STROBE)报告指南进行。
肌肉减少症状态的评估
本研究采用2019年亚洲肌肉减少症工作组(Asian Working Group for Sarcopenia, AWGS)制定的诊断标准23。大量既往研究已充分证实AWGS诊断标准在中国中老年人群中的适用性和可靠性24,25。根据2019年AWGS标准,肌肉减少症的特征是肌肉量减少,同时伴有肌肉力量下降或身体功能受损。
肌肉力量
我们使用握力计26测量双侧手部的握力(HGS)。每位受试者以双侧手各测量两次,测量时将握力计保持在直角(90°)位置。分析时采用双侧手测量值的最大值的平均值,单位为千克(kg)。握力不足的判定标准为女性低于18 kg,男性低于28 kg23。
身体功能测试
在此测试中,要求受试者坐在椅子上,双臂交叉抱于胸前。随后,他们需尽可能快速地连续完成五次从坐位站起再坐下的动作。使用秒表记录完成全部五次重复动作所需的时间(单位为秒)。椅子应紧靠稳定表面放置,受试者仅能使用腿部完成动作,不得借助手臂力量27。若个体在12秒内无法完成五次坐-立测试,则被认为存在身体功能下降23。
骨骼肌质量指数(SMI)
本研究中,四肢骨骼肌质量(ASM)采用2019年亚洲肌少症工作组(AWGS)推荐的测量公式进行计算23。该公式已在本研究中经过严格验证,并在中国人群样本中表现出对ASM估算的高预测准确性28,29。具体公式如下:
ASM = 0.193 × 体重 + 0.107 × 身高 - 4.157 × 性别 - 0.037 × 年龄 - 2.631
使用身高计和数字体重秤分别测量参与者的身高和体重,精确至0.1 cm和0.1 kg。骨骼肌质量指数(SMI)通过将四肢骨骼肌质量(ASM)除以身高的平方(ASM/身高²)计算得出。因此,女性ASM/身高²值<5.4 kg/m2、男性<7.0 kg/m2提示肌肉质量偏低23。
尿酸
由中国疾病预防控制中心(中国疾控中心)的医学培训人员按照标准操作规程采集参与者的血液样本。血液样本采集后立即在 4 °C 下保存,并运送至当地疾控中心实验室或卫生中心进行处理。随后将样本分离为血浆和白膜层组分,血浆分装至 0.5 mL 冻存管中,在 -20 °C 冷冻后运往位于北京的中国疾控中心。样本到达后于 -80 °C 保存,直至在首都医科大学实验室进行分析30。尿酸(UA)水平采用 UA Plus 方法进行定量,结果以毫克/分升(mg/dL)表示22。
协变量
协变量包括社会人口学信息和健康相关特征31,32。社会人口学变量包括年龄和性别,而健康相关协变量包括体重指数(BMI)、吸烟和饮酒史、慢性疾病以及血液检测值。通过问卷收集参与者的健康状况数据,包括高血压、血脂异常和糖尿病等情况。血液检测结果涵盖多种生物标志物,即血红蛋白、肌酐、血尿素氮(BUN)、高密度脂蛋白(HDL)、低密度脂蛋白(LDL)、总胆固醇(TC)、甘油三酯(Tg)和C反应蛋白(CRP)。
统计分析
鉴于CHARLS复杂的抽样设计,所有分析均使用R软件(版本4.2.2)和EmpowerStats 2.0进行。根据血清尿酸浓度将参与者分为四分位数组,最低浓度组作为参照组(四分位数组Q1)。此外,还按年龄进行亚组分析。由于85岁及以上个体数量较少,故将其划分为四个年龄组: <55、55–64、65–75 和 ≥75。本研究中协变量的缺失值均低于5%,采用常见的填补方法(中位数填补或均值填补)进行填补。分析内容包括连续变量的均值、中位数、标准差、范围和四分位数,以及分类变量的频数表。对于非正态分布的连续变量,采用中位数和四分位距(IQR)进行描述。所有分析均使用生物标志物子样本权重,该权重包含对家庭和个体应答的调整,以生成代表中国人群的估计值。为此,使用 R 中的 survey 软件包,调用相关函数如 svymean(),svyvar(),以及 svyquantile() 计算均值、标准差、中位数和四分位数。为检验组间变量特征的差异,对连续变量采用t检验和调查加权Wilcoxon秩和检验,对分类变量的加权百分比采用Rao-Scott卡方检验,以全面描述总体人群特征。采用加权逻辑回归模型分析血清尿酸与肌少症之间的关联。使用加权多变量线性回归模型探讨血清尿酸与握力(HGS)、起立-行走计时测试(FTSST)以及骨骼肌质量指数(SMI)之间的关联。回归模型的构建通过survey软件包中的svyglm()函数实现。模型对多种协变量进行了调整,包括年龄、性别、吸烟史、饮酒史、血液参数(肌酐[Crea]、血尿素氮[BUN]、高密度脂蛋白[HDL]、低密度脂蛋白[LDL]、总胆固醇、甘油三酯、C反应蛋白[CRP])以及慢性疾病(高血压、血脂异常、糖尿病或糖代谢异常)。双侧P < 0.05 被视为具有统计学意义的阈值。
研究对象特征
本研究共纳入10,938名年龄在45至94岁之间的中国参与者,其中肌肉减少症的加权患病率为6.68%。根据血清尿酸(SUA)水平的四分位数,研究队列被分为四组:Q1:≤3.90 mg/dL,Q2:4.00–4.70 mg/dL,Q3:4.80–5.60 mg/dL,Q4:≥5.70 mg/dL。表1展示了按血清尿酸四分位数分层的加权基线特征。结果显示,随着SUA四分位数的升高,加权平均年龄和体重指数(BMI)呈逐步上升趋势(P < 0.001)。此外,尿酸水平的升高与男性比例增加及女性比例减少相关。血液学分析显示,随着SUA四分位数的升高,加权后的肌酐、尿素氮(BUN)、低密度脂蛋白胆固醇(LDL)、总胆固醇(TC)、甘油三酯(TG)和C反应蛋白(CRP)水平均逐步升高(P < 0.001)。相反,加权高密度脂蛋白胆固醇(HDL)水平则随SUA升高而逐步下降(P < 0.001)。在SUA四分位数之间,尿酸水平的升高与吸烟者比例、高血压和血脂异常的加权患病率上升相关,而饮酒者比例则下降。人体测量指标显示,随着SUA四分位数的升高,加权握力(HGS)和骨骼肌质量指数(SMI)均呈逐步上升趋势(P < 0.001),而五次起坐-站立测试(FTSST)时间则逐步缩短(P < 0.001)。此外,肌肉减少症的加权患病率随尿酸水平升高而下降(P < 0.001)。
血清尿酸水平与肌肉减少症
表2展示了血清尿酸(UA)水平四分位数与肌肉减少症之间的关联。以Q1作为参照组。在未调整模型中,总体人群中肌肉减少症的患病率与UA水平呈负相关(Q2:OR 0.660(0.531–0.821),P < 0.001;Q3:OR 0.583(0.464–0.733),P < 0.001;Q4:OR 0.479(0.380–0.603),P < 0.001)。在年龄亚组分析中,Q4四分位组在所有年龄组中均持续表现出负相关关系(45–54岁:OR 0.096(0.021–0.433),P = 0.002;55–64岁:OR 0.370(0.218–0.629),P < 0.001;65–74岁:OR 0.256(0.178–0.369),P < 0.001;≥75岁:OR 0.389(0.241–0.628),P < 0.001)。在对年龄、性别、吸烟与饮酒史、血液参数(肌酐、尿素氮、高密度脂蛋白、低密度脂蛋白、总胆固醇、甘油三酯、C反应蛋白)以及慢性疾病(高血压、血脂异常、糖尿病或糖代谢异常)进行综合调整的模型中,肌肉减少症总体患病率与UA水平之间的负相关关系仍然存在(Q2:OR 0.687(0.528–0.896),P = 0.005;Q3:OR 0.606(0.463–0.793),P < 0.001;Q4:OR 0.477(0.352–0.647),P < 0.001)。此外,在65岁及以上的个体中,肌肉减少症与Q4四分位组之间也表现出显著的负相关关系(65–74岁:OR 0.357(0.230–0.553),P < 0.001;≥75岁:OR 0.485(0.275–0.856),P = 0.013)。
血清尿酸水平与握力/HGS、起立行走测试时间/FTSST及肌肉质量指数/SMI
表3展示了血清尿酸水平与握力/HGS、起立行走测试时间/FTSST及肌肉质量指数/SMI之间经调整的加权多变量线性回归模型分析结果。第一四分位组(Q1)作为基线参照组。第二四分位组(Q2)的血清尿酸水平显著升高(β 1.474(0.787–2.161),P < 0.001);第三四分位组(Q3)(β 2.913(2.328–3.498),P < 0.001)和第四四分位组(Q4)(β 4.210(3.561–4.859),P < 0.001)也均显著升高,且与握力/HGS呈正相关,表现出剂量效应关系。在年龄亚组分析中,Q4组在所有年龄组中均保持正相关关系(45–54岁:β 4.913(3.605–6.220),P < 0.001;55–64岁:β 4.170(3.163–5.176),P < 0.001;65–74岁:β 3.253(2.070–4.435),P < 0.001;≥75岁:β 2.825(0.978–4.671),P = 0.003)。血清尿酸水平在Q3组(β -0.448,95% CI -0.684–(-0.213),P < 0.001)和Q4组(β -0.410,95% CI -0.654–(-0.166),P < 0.001)与起立行走测试时间/FTSST相比差异具有显著性。
在年龄亚组分析中,Q4组在45-54岁(β -0.721,95% CI -1.086 至 -0.357,P < 0.001)和65-75岁(β -0.692,95% CI -1.187 至 -0.196,P = 0.006)均显示出显著的负相关。血清UA水平在Q2组(β 0.089(0.026–0.152),P = 0.006)、Q3组(β 0.159(0.094–0.224),P < 0.001)和Q4组(β 0.270(0.199–0.342),P < 0.001)均较高,且与SMI呈正相关,并存在剂量效应。在年龄亚组分析中,Q4组在所有年龄组中均保持正相关(45–54岁:β 0.294(0.160–0.428),P<0.001;55–64岁:β 0.223(0.111–0.336),P < 0.001;65–74岁:β 0.315(0.200–0.430),P < 0.001;≥75岁:β 0.342(0.151–0.533),P < 0.001)。
数据可用性:
本研究分析的数据来自公开可用的数据集,可通过以下链接访问: https://charls.pku.edu.cn/。

图 1:本研究所用样本的流程图。 请点击此处查看此图的放大版本。
表1:基于血清尿酸四分位数的研究人群加权特征 请点击此处下载该表格。
表2:在调整相关混杂因素后评估年龄分层的肌肉减少症与尿酸四分位数的加权逻辑回归模型。 缩写:CI,置信区间;OR,比值比。注:调整后模型包括年龄、性别、吸烟史、饮酒史、血液参数(肌酐、尿素氮、高密度脂蛋白、低密度脂蛋白、总胆固醇、甘油三酯、C反应蛋白)以及慢性疾病(高血压、血脂异常、糖尿病或糖代谢异常)请点击此处下载该表格。
表3:在调整相关混杂因素后,测试不同年龄组中尿酸(UA)四分位数与握力(HGS)、五次起坐试验(FTSST)和骨骼肌指数(SMI)之间关系的加权多变量线性回归模型。缩写:HGS,握力;FTSST,五次起坐试验;SMI,骨骼肌指数;注:调整模型包括年龄、性别、吸烟史、饮酒史、血液参数(肌酐、尿素氮、高密度脂蛋白、低密度脂蛋白、总胆固醇、甘油三酯、C反应蛋白)以及慢性疾病(高血压、血脂异常、糖尿病或血糖异常)。请点击此处下载该表格。
本研究利用具有全国代表性的数据进行横断面研究,旨在评估中老年人群中血清尿酸(UA)水平与依据亚洲肌少症工作组(AWGS)2019标准诊断的肌少症之间的关联。研究对象为年龄在45至94岁之间的中国成年人,其肌少症患病率为6.68%。在调整了混杂因素后发现,血清UA水平与肌少症呈负相关。血清UA水平的四分位分组与握力(HGS)及骨骼肌质量指数(SMI)之间存在显著正相关。这些结果提示,UA可能对肌肉质量和力量具有保护作用。结合既往研究结果18,33,34,35,36,37,本分析结果与一项在中国西部开展的研究一致,表明在50岁及以上人群中,UA水平升高与肌肉质量和握力均呈显著正相关33。此外,一项基于社区的研究发现,在40至75岁的中国人群中,血清UA水平升高与肌肉质量增加之间存在相关性34。另有针对936名老年人的横断面回顾性研究显示,高尿酸血症与肌少症呈负相关,而与肌肉质量和肌肉力量呈正相关37。上述研究提示,血清UA水平升高可能在维持骨骼肌质量和力量方面发挥关键作用,从而延缓肌少症的进展。
在年龄亚组分析中,观察到尿酸(UA)与肌肉减少症及其组分(握力(HGS)、五次起坐测试(FTSST)、骨骼肌指数(SMI))之间的关系在不同年龄组中存在差异。本研究显示,在45-54岁和55-64岁年龄组中,UA与HGS呈显著正相关,并表现出剂量-反应效应。在65-75岁和≥75岁年龄组中,仅Q4组(≥5.70)显示血清尿酸水平与HGS之间存在正相关。该结果与现有文献的结论不同。根据一项在韩国开展的研究,血清UA浓度与青中年人群的HGS之间无统计学意义,但在60岁及以上人群中,这两者之间存在正相关38。相反,另一项研究记录了UA对较年轻人群HGS的不利影响39。造成这种差异的可能原因在于不同研究之间基线血清UA水平的差异以及年龄分组标准的不同。未来的研究应重点关注年龄对UA与肌肉力量之间关系的影响。
此外,Q4组血清尿酸水平在65-75岁及75岁以上的老年人群中与肌少症呈显著负相关。相比之下,在45-54岁和55-64岁的中年人群中未观察到相关性。然而,Q4组的血清尿酸水平在45岁及以上的成年人中与握力(HGS)和骨骼肌指数(SMI)均呈显著正相关。综上所述,较高的血清尿酸水平似乎对45岁及以上成年人的握力和骨骼肌指数具有保护作用,但在45至64岁的中年人群中未观察到对肌少症的显著预防作用。一种合理的解释是,较高的血清尿酸水平可能在45岁以上个体中保护肌肉质量和力量。然而,由于中年人通常处于体力巅峰阶段,肌肉质量和力量较高,仅有极少数人符合肌少症的诊断标准,导致高血清尿酸水平与肌少症之间的关系不明显。随着年龄增长,机体代谢水平逐渐下降,各种功能障碍和慢性炎症的发生率上升,导致肌肉质量和力量进行性下降。在此阶段,较高血清尿酸水平对肌少症的保护作用变得显著。一项基于英国12个普通人群队列、涉及49,964名参与者的研究所显示,肌肉力量在成年早期达到峰值,维持至中年,之后从中年开始逐渐下降40。这可以解释不同年龄组中较高血清尿酸水平与肌少症之间关系的差异,但仍需更多纵向研究加以验证。
血清尿酸(UA)水平、握力(HGS)与肌肉量之间存在正相关关系的具体机制尚不明确。一个可能的原因在于尿酸强大的抗氧化能力。已有研究表明,骨骼肌产生的活性氧和活性氮物质(RONS)会加剧多种变性性肌肉疾病的氧化损伤41。此外,在与年龄相关的肌肉量减少背景下,氧化还原稳态的紊乱会影响多种转录因子和激酶的功能,包括FOXO、NF-κB以及MAPK家族成员,从而阻碍肌肉再生过程42。研究人员发现,能够降低活性氧水平的细胞内环境有利于肌生成43,提示抗氧化剂可能在缓解肌肉衰老和功能衰退方面具有治疗潜力。尿酸是主要的水溶性抗氧化剂,约占血浆抗氧化能力的60%44。尿酸可通过清除自由基,并与过渡金属离子(尤其是铁离子)形成稳定的复合物来发挥作用。这种双重机制使UA能有效抑制抗坏血酸的氧化以及脂质的过氧化,减轻由氧化应激引起的损伤45。因此,尿酸的存在可增强机体的抗氧化能力,从而保护肌细胞免受氧化损伤,维持骨骼肌的功能。然而,尿酸水平升高与肌肉量及肌肉功能之间正相关的具体机制仍未完全阐明,有待进一步研究。
本研究利用了来自CHARLS数据库的带有样本权重的血液检测数据,确保研究结果具有对中国总体人群的代表性。使用调查权重可校正复杂的抽样设计,从而提高研究结果的准确性和可推广性46。采用CHARLS进行数据分析的一个重要优势在于其全国代表性强、数据质量高以及样本量大47。这些因素增强了研究的稳健性,使得所得结论能够可靠地反映中国广泛而多样化人群的健康状况。大样本量还提高了统计效能,支持在中国不同年龄组中开展细致的亚组分析。展望未来,CHARLS为后续纵向研究提供了宝贵机会,有望进一步揭示各类健康因素(如血清尿酸水平)对衰老相关疾病长期影响的深入见解。然而,本研究也存在一定的局限性。首先,CHARLS缺乏关于影响尿酸水平的药物使用信息,以及饮食习惯和结构的数据,因此无法对部分相关混杂因素进行调整。其次,在肌少症的测量方面,肌肉质量是依据亚洲肌少症工作组(AWGS)推荐的公式计算得出,属于间接测量方法,而非双能X线吸收测定法这一金标准。最后,由于本研究为横断面研究,无法阐明变量之间的因果关系,仅能展示特定时间点上各因素之间的关联,因而难以推断尿酸水平与肌少症之间的因果方向。
综上所述,本研究表明,在45岁及以上的中国成年人群中,血清尿酸水平升高与肌肉质量和肌肉力量之间存在显著的正相关关系。此外,维持较高的血清尿酸水平可能有助于延缓肌少症的进展,尤其是在65岁及以上的个体中。鉴于肌少症在全球范围内普遍存在且老年人口数量不断增加,本研究结果对于肌少症的预防及促进健康老龄化具有重要意义。
作者声明无利益冲突。
本研究基于中国健康与养老追踪调查(China Health and Retirement Longitudinal Study, CHARLS)的基线数据。我们诚挚感谢辛勤工作的CHARLS团队所做出的贡献。
本研究由以下项目资助:中国博士后科学基金(资助号:2022TQ0398)、重庆医科大学附属第一医院基金(PYJJ2022-01)、重庆市重点学科及区域重点医学学科发展项目(0201 [2023] 160 202412)、国家自然科学基金(资助号:82401836)以及重庆医科大学附属第一医院研究生导师团队项目(CYYY-DSTDXM-202407)。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 中国健康与养老追踪调查(CHARLS) | 北京大学 | https://charls.pku.edu.cn/. | |
| EmpowerStats | https://www.empowerstats.net/en/ | 版本 2.0 | |
| R | https://www.r-project.org/ | 版本 4.2.2 |