膝骨关节炎是临床上最常见的一种骨关节炎，主要临床表现为膝关节疼痛和关节活动障碍^[1]。膝骨关节炎发病率和治疗费用较高，对中老年患者的生活质量和经济负担造成了较大影响^[2]。膝骨关节炎是中医骨伤科及推拿科的优势病种之一，中医药治疗膝骨关节炎疗效确切、费用合理、不良反应少，近年来受到了学术界的持续关注^[3]。本研究应用陈超美博士基于Java语言开发的CiteSpace可视化分析软件^[4]，对近10年推拿治疗膝骨关节炎相关文献进行结构性研究和分析，为了解该领域的发展趋势和推进未来研究的开展提供参考。

在中国知网数据库中以主题词“推拿”或“按摩”或“手法”并“膝骨关节炎”进行检索，检索时间为2010年1月1日至2020年5月11日，共检索出文献438篇。剔除会议报道和重复文献，一共收入有效文献426篇。

在中国知网高级检索页面，勾选符合条件的文献并以Refwork格式导出，将文件名保存并转换为“download_XXX.txt”格式。在CiteSpace 5.6 R3软件中进行知识图谱绘制，按照年发文量、作者、发表机构、关键词等节点分型，统计年发文量，进行作者、机构和关键词共现分析。网络中节点的重要性一般用中心度（centrality）来表示，选取中心度≥0.1且频次＞10的关键词作为高频关键词。对426篇文献的关键词进行聚类分析，选用对数似然比（log-likehood ratio，LLR）算法对聚类标签逐一归纳，得到目前推拿治疗膝骨关节炎主要研究内容的分类。设置切片时间为1年，选取其中被引频次最高的50篇文献，将年份和聚类作为参考值绘制时区视图（time zone view），观察各聚类近10年的发展演变。使用pathfinder、pruning sliced network、pruning the merged network等方法对图谱进行处理^[5]。用突现词（Burst term）检测工具（Burstdetecting）对关键词进行提取，以2年为最小区间（minimum duration）进行探测，获得突现词。

对纳入的426篇文献按照发表时间段进行统计，2010年至2020年推拿治疗膝骨关节炎的相关文献发文量分别为15、19、27、23、54、44、51、45、65、71、12篇，总体呈增长趋势。其中，2014年的发文量（54篇）是2013年的2.3倍，之后3年略有回落，2018年发文量再次升高（65篇），2019年发文量为近10年之最（71篇）。

426篇文献共有302名作者，作者合作网络图谱见图 1A。从连线节点来看，最主要的作者合作群为上海中医药大学附属岳阳中西医结合医院的龚利、房敏、严隽陶、孙武权等的合作群，其中龚利发文14篇，房敏发文7篇，孙武权、樊远志、吴耀持、严隽陶各发文5篇；其他作者合作群中，曲崇正和刘姣各发文4篇。

图 1 推拿治疗膝骨关节炎研究作者合作网络（A）与机构合作网络（B）分析

推拿治疗膝骨关节炎相关研究机构较为分散（density＝0.031 7），主要发表机构为上海中医药大学附属岳阳中西医结合医院（19篇）和广州中医药大学第三附属医院（9篇），其次为上海中医药大学（7篇）、天津中医药大学第一附属医院（6篇）和上海交通大学附属第六人民医院（5篇）。机构之间合作以同一地区或同一系统内的合作较为紧密，而异地的各机构之间合作较少（图 1B）。

通过对关键词的共现分析，得到频次＞10且中心度≥0.1的关键词13个（表 1）。对426篇文献的关键词进行聚类分析，选取被引频次最高的50篇文献，得到的知识图谱由92个节点和290条连线组成（图 2）。其模块性聚类数Q为0.456 4（＞0.3），轮廓性聚类数S为0.588 9（＞0.5），表明该次聚类分析所划分的社团结构显著，聚类合理。本次统计共出现7个聚类，聚类标签的规模按照数字递增的方式对应值递减。对这7个聚类标签逐一归纳，得到目前推拿治疗膝骨关节炎主要研究内容的分类（表 2）。将年份和聚类作为参考值，得到推拿治疗膝骨关节炎的时区视图，从图中可以观察各聚类近10年的发展演变（图 3）。

表 1 (Tab 1)

表 1 推拿治疗膝骨关节炎的高频关键词 排序 频次 中心度 年份 关键词 1 151 0.15 2010 推拿 2 149 0.30 2010 膝骨关节炎 3 103 0.21 2010 膝骨性关节炎 4 89 0.21 2011 膝关节骨性关节炎 5 30 0.19 2011 骨关节炎 6 25 0.10 2015 综述 7 24 0.12 2010 推拿手法 8 23 0.11 2010 推拿疗法 9 22 0.16 2012 针灸 10 20 0.36 2013 温针灸 11 19 0.30 2012 临床研究 12 17 0.16 2014 针刺 13 13 0.24 2013 临床疗效

表 1 推拿治疗膝骨关节炎的高频关键词

图 2 推拿治疗膝骨关节炎关键词共现分析图谱

表 2 (Tab 2)

表 2 基于LLR算法的推拿治疗膝骨关节炎关键词聚类分析 聚类号 规模 轮廓值 年份 LLR 0 13 0 2015 手法推拿 1 12 0 2015 小针刀 2 11 0 2016 中医药 3 10 0 2015 痹症 4 10 0 2013 按摩 5 10 0 2015 Myd88 6 9 0 2014 中医药治疗 LLR：对数似然比；Myd88：髓样分化因子88

表 2 基于LLR算法的推拿治疗膝骨关节炎关键词聚类分析

图 3 推拿治疗膝骨关节炎的时区视图

用突现词检测工具对关键词进行提取，以2年为最小区间探测，共得到4个突现词，分别为“推拿手法”“骨关节炎”“综述”“疗效”。其中，热度持续3年的有“骨关节炎”（2011年至2013年）和“疗效”（2018年至2020年）。

本研究通过检索中国知网数据库，利用CiteSpace软件分析了近10年推拿治疗膝骨关节炎的相关文献，结果显示近10年有2个文献发表的高峰期，分别为2014年和2019年。纵观近10年的文献发表量，以完整的自然年来计算，2019年的发文量为2010年的4.7倍。上述结果提示随着社会老龄化，中医推拿治疗膝骨关节炎的研究热度在不断升高；研究成果的产出具有一定的周期性，符合科学研究的发展规律。

根据文献作者和研究机构共现分析，推拿治疗膝骨关节炎研究发文量最多的机构为上海中医药大学附属岳阳中西医结合医院，其中以龚利、房敏、孙武权、樊远志等研究者发表的文献数量最多。龚利（通信作者）等^[6]在丁氏推拿流派特色手法的基础上，结合现代康复医学关节松动理念创造出“坐位调膝”推拿手法，该手法治疗膝骨关节炎不但可以缓解疼痛，还对膝关节周围肌群平衡具有恢复作用。此外，根据机构网络图谱分析发现，推拿治疗膝骨关节炎文献发表机构主要集中于我国东南部地区，尤其以上海、广东为主，且各机构之间离散度较高，缺乏合作。这与《膝骨关节炎阶梯治疗专家共识（2018年版）》^[7]中指出的西南、西北地区膝骨关节炎患病率高于东南部沿海地区的调查结果恰好相反。这可能与我国各地区医疗科研水平发展不均衡有关^[8]，同时也提示作者及各研究主体之间应加强合作，为推拿治疗膝骨关节炎相关研究的进一步开展提供更有利的条件。

在研究趋势和热点预测方面，从时区视图中可以发现早期治疗手段主要以推拿手法为主，后逐步发展为联合多种中医特色疗法（如内服汤药、针刺、小针刀、中药熏蒸）^[9-15]或现代康复疗法^[16-18]相结合；观察指标从疼痛积分、关节功能评估、步态分析向更多元评估方向发展，如关节液脂肪因子（visfatin及chemerin）含量、髓样分化因子88（myeloid differentiation factor 88，Myd88）等^[19-21]被引入机制研究。根据突现词探测分析发现，近10年研究推拿治疗膝骨关节炎的热点仍然局限于不同手法的临床疗效评价和文献综述方面，而推拿治疗膝骨关节炎的机制研究、动物实验和基于大数据平台或新型科学方法的文献计量学研究较少。由此可见，推拿治疗膝骨关节炎目前的研究过于保守，需引起重视。