当前软件版本及合法许可证：KLFromRecordingDays

当前软件版本指向此版本可合法软件许可证计算平台/操作系统安装要求依赖性如果有用户手册的链接-如果正式出版，请在参考列表中包括问题支持电子邮件1.1https://github.com/soderstromk/KLFromRecordingDays/tree/master/cx_freeze/KLFromRecordingDaysMITMicrosoft Windows，64位需要Sound Analysis Pro 2011，http://soundanalysispro.com/如果使用Sound Analysis Recorder生成 *.wav文件，则该实用程序ParseSAPRecorderWaves非常有用，https://github.com/soderstromk/KLFromRecordingDays/tree/master/cx_freeze/ParseSAPRecorderWaves文档和手册包含在手稿中。安装说明包含在GitHub存储库README.mdblob/master/README.mdhttps://github.com/soderstromk/KLFromRecordingDays/文件中soderstromk@ecu.eduSoftwareX 6（2017）271原始软件出版物用于客观比较数周音频记录中的声音声学特征的软件：KLFromRecordingDaysKen Soderstrom*，Ali Alalawi药理学和毒理学系，Brody医学院，东卡罗来纳大学，Greenville，NC 27834，美国ar t i cl e i nf o文章历史记录：2017年6月27日收到2017年9月11日收到修订版2017年10月9日接受保留字：声乐学习声学分析Kullback–Leiblera b st ra ctKLFromRecordingDays允许测量人声声学特征的2D概率分布之间的Kullback-Leibler（KL）距离。更大的KL距离的措施反映了增加语音分歧的发声比较。该软件已被用于比较由声音分析记录器2011制作的鸣禽发声前和后药物和手术操作的 *.wav文件记录。来自单个动物的 *.wav格式的录音首先按录音日期组织到子目录中，然后使用Sound Analysis Pro 2011（SAP）将其分割为发出的单个音节和这些音节的声学特征。KLFromRecordingDays使用SAP输出到MySQL表的音节声学特征数据来生成和比较“模板"（通常是治疗前）和”目标"（通常是治疗后）概率分布。这些分布是每个音节的持续时间（作为x轴）与SAP测量的13个其他声学特征中的每一个的一系列虚拟2D图对于该音节（作为y轴）。“模板"和”目标“概率分布之间的差异通过计算KL距离来确定每个声学特征，KL距离是目标2D分布图案与模板的2D分布图案的发散度的度量。计算每个记录日的KL距离和所有声学特征的平均KL距离，并输出到Excel电子表格。然后，可以汇总各治疗组的个体受试者的结果数据，并以图形方式进行总结，用于统计学比较。由于SAP生成的MySQL文件是直接访问的，因此避免了与电子表格输出相关的数据限制，并且可以一次性客观地分析数周内的语音输出总量。该软件可用于测量药物对鸣禽发声的影响，以及评估发声运动皮层区域受损后的恢复情况。它可能在使用其他物种的研究中有用，并作为语言治疗的一部分，跟踪孤立地产生不同语音的进展。©2017作者。由爱思唯尔公司出版这是CC BY许可下的开放获取文章（http://creativecommons.org/licenses/by/4.0/）中找到。软件元数据通讯作者。电子邮件地址：soderstromk@ecu.edu（K. Soderstrom）。缩略语：SAR，Sound Analysis Recorder 2011; SAP，Sound Analysis Pro 2011http://dx.doi.org/10.1016/j.softx.2017.10.0032352-7110/©2017作者。由爱思唯尔公司出版这是CC BY许可下的开放获取文章（http://creativecommons.org/licenses/by/4.0/）。可在ScienceDirect上获得目录列表SoftwareX期刊主页：www.elsevier.com/locate/softx*272K. Soderstrom，A.Alalawi/SoftwareX 6（2017）271当前代码版本此代码版本使用的代码/存储库的永久链接使用的代码版本系统使用的软件代码语言编译要求、操作环境依赖性（如有）开发人员文档/手册问题支持电子邮件1.1https://github.com/ElsevierSoftwareX/SOFTX-D-17-00050MITgitPython 3.5需要Microsoft Windows，64位文件和手册被纳入手稿soderstromk@ecu.edu代码元数据1. 介绍大约4%的儿童在6岁时会受到言语障碍的困扰，并且通常继发于其他更严重的发育障碍，包括难治性儿童癫痫和自闭症[1]。这些疾病通常用言语语言物理治疗来治疗。有效的药物干预仍然不可用，可能是由于缺乏适当的临床前动物模型。很少有动物会学习以类似于人类的方式进行语音交流。鸣禽是这一小群声乐学习者中的一员。因此，我们已经开始开发一种鸣禽，斑胸草雀，作为实验室动物模型，以评估药物对改善声音学习和从CNS损伤中恢复的作用。该模型依赖于方法的发展，以客观地比较几周内治疗前和治疗后的发声质量。存在用于记录和分析动物发声的声学特征的优秀软件（例如Sound Analysis Recorder和Sound Analysis Pro 2011 [SAP，[2]]），并且已经描述了通过KL距离测量比较发声语音质量的良好方法[3]。开发当前软件的目标是以一种允许有效分析持续数周的实验中产生的全部发声的2. 问题和背景2.1. 解决问题所描述的软件解决了将量化音位学的KL距离方法（由[3]开发）应用于由SAP计算的数百万个声学特征测量的问题，这些声学特征测量是从动物在多周实验期间记录的每个音节中得出的从音频记录到KL距离测量的最终输出，分析都是客观完成的。2.2. 背景SoundAnalysis Pro（SAP）是一款功能强大的免费开源软件，用于研究鸣禽的发声发展[2]。它包括允许以 *.wav格式记录和管理录音的模块，以及分析软件，该分析软件自动将声音分割成发出的单个音节，并通过计算14个声学特征（音节持续时间、平均振幅、平均音高、平均FM、平均AM 2、平均熵、音高的平均优度、平均频率、音高方差、FM方差、熵方差、音高方差的优度、平均频率方差、AM方差）来分析其频谱结构。SAP使用开源MySQL关系数据库系统来管理数据。SAP提供以原生MySQL表或Excel电子表格格式导出这些数据。其他人已经开发出一种很好的方法来使用SAP电子表格输出进行统计，使用KL距离测量比较鸣禽语音学随时间的变化[3，4]。这些方法已经成功地应用于研究对发声学习和运动产生重要的大脑区域的手术操作的效果[4在使用鸣禽作为临床前动物模型来筛选药物对发声的影响时，我们最初的目的是应用[4]已经开发并免费提供的优秀SongSeq软件（https://www.math.fsu.edu/~bertram/software/birdsong/JNM_12/，最后访问时间为2017年9月5日）。该软件以类似于我们开发的KLFromRecording- Days程序的方式计算KL距离。此外，它还提供了一种强大的方法来区分单个音节，并分析它们在鸣禽发声中以一致的顺序产生的程度。我们继续在我们的项目中使用SongSeq的序列分析功能。在SongSeq的KL距离测量能力的情况下有两个问题使该功能对我们的工作不切实际。第一个问题与SongSeq的KL距离的产生有关，每次测量一个概率分布。也就是说，当使用音节持续时间作为X轴时，剩余的13个声学参数中的仅一个可以被用作y轴以生成跨天比较的分布。额外的参数需要单独的、顺序的分析，由于我们项目的规模，这些分析很耗时，并且在组织和汇总结果方面产生了问题。我们在应用SongSeq生成KL距离度量时遇到的第二个也是更重要的问题是，它使用SAP生成的电子表格进行数据输入。这些电子表格是较旧的 *.xls文件格式，最多限制为65，536行。由于SAP处理的每个音节的声学测量占据了一个电子表格行，并且在一天的记录期间产生数十万个音节并不罕见，为了客观地分析每个音节发出的数据，我们需要建立KL距离测量软件，该软件直接访问SAP采用的MySQL数据库系统。MySQL数据库表的大小仅受所用计算机的内存资源限制。鉴于Python的开源性质和PyMySQL和openpyxl包分别与MySQL和Ex- cel接口，我们选择这种编程语言（版本3.5）来开发我们的应用程序。该软件被开发为在64位Windows系统上运行。3. 软件框架3.1. 软件构架该软件由一个图形用户界面控制，该界面来自Tkinter Python软件包（见图1）。①的人。该界面收集来自用户的信息，包括：用于存储KL距离测量结果的Excel电子表格的目录;用于分析的SAP生成的音节表;用作“模板”的记录日（剩余K. Soderstrom，A.Alalawi/SoftwareX 6（2017）271273Fig. 1. KLFromRecordingDays GUI。Tkinter-derived窗口允许选择用于存储Excel格式输出的目录;要分析的SAP生成的MySQL音节表;以及用作与目标记录天数进行比较的模板的天数。控制分析各方面的参数可在窗口的下半部分获得，并在第3.2.1节用户控制变量中详细解释。before *.wav文件用于分析;用于空单元的平滑因子（下面解释）;估计音节是否是基元的一部分的音节间隔截止（ms）;是否用估计的音节类型修改SAP音节表以及是否从分析中排除非基元音节。该程序的流程图如图所示。 二、该程序使用PyMySQL查询SAP数据库以获得音节表列表。一旦选择音节表进行分析，就创建在所选音节表中表示的所有记录日的此列表显示在一个选择框中，用户可从中指定哪一天用作模板日通常是治疗前或操作前基线记录。选择框中列出的未指定包含在“模板”中的每一天都用作“目标"天，用于与模板进行输入信息后，用户激活界面上的这将启动对表和模板条目的检查，在出现错误时将控制权返回给GUI如果图二、 流程图说明了KLFrom-RecordingDays的体系结构和执行顺序。不是表格或模板选择遗漏，“”MainKL“”然后调用“”ModifyTableWithDay“”，将新列添加到选定的音节数据库中，并为表中的每个音节输入记录日期信息。如果用户选择了为了仅分析基元类型音节（滤出非基元音节）或为了计算音节类型但使用全部用于分析（计算基元类型但使用全部），调用函数来计算每个音节的前音节间隔和后音节间隔，并将该信息添加到音节表（Syll_Duration）。然后调用函数来创建新的模板和目标MySQL 表 ， 这 些 表 包 含 来 自 相 应 记 录 日 期 的 音 节 信 息（MakeTemplateTable和MakeTargetTable函数）。然后，目标表用于为每个目标记录日创建一系列子表（TargetDayTables）。将这些目标日特定表中的每一个与模板表和为每个声学参数生成的KL距离进行比较。一旦创建了模板和一系列记录特定于日期的目标MySQL表，就为每个目标日期输入一个for循环。进入第二个嵌套的for循环，该循环迭代将用作2D数组中的y轴的13个声学参数中的每一个（上面列出的，不包括总是用作x轴的音节持续时间计算的每个声学参数的值由SAP为每个音节从模板和目标访问使用传递模板或目标日期表名称以及要访问的声学参数的名称的“TableParameters "函数创建日期表。然后，该函数使用PyMySQL查询每个音节的参数值的表，并将 这 些 值 作 为 一 个 列 表 存 储 在 由 参 数 名 称 键 控 的 字 典 中 。“”TableParameters“” 返回此参数名称键字典。最终生成了四个“TableParameters”字典：模板和目标日期表各两个;二维图的x轴;另一个为y轴。参数键控274K. Soderstrom，A.Alalawi/SoftwareX 6（2017）271==-=然后将模板和目标日2D图的字典数据传递给类对象，以计算该参数“GenerateKL”的KL距离。GenerateKL为每个声学参数创建一个类对象。它使用上述四个“TableParameter”字典、轴划分值和平滑因子进行初始化。从然后，这些最小值和最大值与模板或目标参数值以及轴分割值一起传递给Python numpy包函数“histogram2d”。Histogram2d使用等宽条创建一个虚拟的二维单元阵列。默认的轴划分值为15，这些单元格152225然后，Histogram2d将参数值数据中的每个音节解析为与适当的音节持续时间和参数值相对应的2D单元。落入虚拟2D阵列的每个单元内的音节的总和由histogram2D作为关联阵列返回。为模板和目标参数值生成音节的2D关联阵列。用户可选择的平滑因子（默认值为10−6）添加到每个像元以避免分裂由空细胞引起的零。模板和目标概率数组是通过将每个单元格值除以所有单元格的总和来生成的结果数组1的所有单元格的总和。因此，每个单元表示模板或目标音节将被分配给它的分数概率。生成模板和目标概率数组后该函数采用由[3]描述的KL距离计算公式 。 请 注 意 ， KLCalc 函 数 的 输 出 （ 下 面 的 代 码 片 段 ） 与scipy.stats.entropy函数的输出等效KLCalc函数比较目标概率数组的2D模式与模板的2D模式的发散度。值越大表示分歧越大。该函数返回在所评估的参数的模板和目标分布之间计算的KL距离，并且该值被分配给以参数名称为键的一旦计算出所有参数的KL距离值并将其添加到字典中，则该记录日的KL距离测量完成，并且将该日然后执行迭代到下一个记录日。处理完所有记录日的数据后，将按记录日键入的完整结果字典传递给一个函数，以将结果写入Excel电子表格：“WriteXL”函数使用openpyxl包，该包允许通过Python操作Excel电子表格对象。将工作目录更改为GUI中指定的目录，创建新的Excel工作簿，并按行写入每个目标日的所有13个声学参数的KL距离结果还记录了关于记录的动物的信息、GUI设置以及针对模板和目标日分析的音节数（参见图2）。电子表格文件由动物ID、模板使用的日期和目标天数命名，并保存在工作目录中。3.2. 软件功能该软件提供了一个度量标准（KL距离），说明动物发声的光谱特征在一天中如何变化。它将这些值存储在用户指定目录中的Excel电子表格中。这些措施是有用的研究药物的语音效果，手术或其他操作在一系列天或几周。3.2.1. 用户控制变量该软件由GUI控制，该GUI收集来自用户的信息，如上文第3.1节“体系结构”所述。除了用于模板的分析和记录天数的表之外，用户还可以指定x轴和y轴将被分隔成的条格或分区的数量。这些箱用于形成2D概率阵列的单元，每个音节基于音节持续时间（x轴）和被迭代的特定参数（在mainKL函数中，y轴）被分配到该单元。轴划分的数量越大，模板和目标概率阵列之间的发散程度可能越高，这将倾向于产生更大的KL距离值。当细胞的数量接近1时，KL距离接近0。因此，这个参数应该优化不同的实验，产生不同的轴范围和音节数分析。使用斑胸草雀作为动物模型，我们和其他人[6]发现，15轴分割的默认值对运动皮层的操作产生良好的敏感性，同时保持一致的治疗前基线值。GUI还允许指定录音中接受分析的最小音节数。音节很少的录音通常不是歌曲的爆发，而是呼叫，笼子的砰砰声，翅膀的颤动或其他外来的声音。我们已经发现，默认设置为每个记录最少五个音节，这消除了许多非歌曲的回合记录，但不排除那些包含完整主题的记录。母音是鸣禽按固定顺序发出的一系列音节。主调音节通常是在短时间的沉默中产生的。因此，作为基元的一部分产生的音节通常可以使用基元间隔截止来可靠地估计并与呼叫、介绍性注释和其他噪声区分开。此参数可由用户控制，默认值为35 msec。个别动物在音节之间的最大间隔方面有所不同，因此应针对每个科目优化该参数。使用基元间隔截止，用户可以选择仅包括估计为用于分析的基元的一部分的音节。或者，所有音节都可以用于分析，在这种情况下，由SAP分析的每个音节将用于计算KL距离值，包括呼叫和介绍性注释。注意，选择使用所有音节避免了触发估计基元类型和修改音节表的功能，从而提高了信息分析速度。为了满足用户可能希望使用所有音节进行KL距离分析，但仍具有预测的基序类型并添加到音节表的情况，GUI上提供了“计算基序类型但使用全部”选项（图11）。①的人。3.3. 示例代码段代码片段1呈现类对象（GenerateKL），其计算单个记录日的特定声学参数的KL距离这个类的操作在上面的体系结构部分3.1中有详细描述4. 执行情况和实证结果4.1. 实现细节SAP 应 在 KLFromRecordingDays 之 前 安 装 这 也 将 确 保 安 装MySQLSAP可从以下网址获得：http：//soundanalysispro.com/ （ 最 后访 问 日期 ： 2017 年 6 月 27 日 ） 。Sound Analysis Pro Recorder（SAR）软件包也可在同一网站上找到 。 KLFromRecordingDays 和 ParseSAPRecorderWaves 可 作 为Python生成的可执行文件使用包年龄cx_freeze对的GitHub仓库https://github.com/K. Soderstrom，A.Alalawi/SoftwareX 6（2017）271275soderstromk/KLFromRecordingDays.git（last accessed 6/27/2017）.要安装，请按照README.md文件中的说明进行这两个可执行文件都应该在任何64位Windows计算机上运行。源代码也可以在GitHub存储库中获得。使用KLFromRecordingDays的顺序步骤如图3的流程图所示。为了适应我们实验室中使用的记录系统文件格式，我们首先使用一个实用程序处理SAR生成的记录，该实用程序根据276K. Soderstrom，A.Alalawi/SoftwareX 6（2017）271\==××图三. 分析步骤总结。录音是在几天内与声音分析专业2011录音机。这些记录在 *.wav文件格式处理到一系列的子目录，根据记录的日子，并根据指定的文件名格式重命名。然后使用Sound Analysis Pro 2011的批处理功能对记录的每个音节进行分段并测量第2.2节背景中描述的一系列声学特征。每个音节的这些声学特征度量存储在MySQL音节表中，然后KLFromRecordingDays和template以及一系列日常目标表访问该表被创造出来。每个目标日的每个声学参数用于创建2D概率图，该2D概率图用于与针对模板生成的可比图进行比较。KL距离，模板和目标2D阵列偏离的程度计算并存储在Excel电子表格中。见图4。图形用户界面的配套程序ParseSAPRecorderWaves重命名文件根据我们的约定（在图中描述。3，第二个面板），并按记录日期将它们移动到子目录中。此实用程序具有将重命名的文件存储在不同目录中（默认情况下，在源目录中创建记录日子目录）和复制文件而不是移动文件（默认情况下，文件，因为它们通常很大，并且它们的数量消耗大量的磁盘空间）。输入到KLFromRecordingDays中的文件必须根据因为动物ID和记录日期是从文件名访问我们的常规格式：*-##-#.wav，其中 * =动物ID，##=记录日期，#=*.wav 文 件 编 号 。 这 个 实 用 程 序 软 件 被 称 为“ParseSAPRecorderWaves"，可以在cx_freezeParseSAPRecorderWaves GitHub 中 找 到 （ 见 图 1 ） 。 4 ） .KLFromRecordingDays依赖于SAP从该文件格式的 *.wav文件中确定音节声学特征，以读取动物ID和记录日期信息。还提供了此实用程序的源代码，以便可以进行更改以适应其他4.2. 实证结果根据正在处理的MySQL表的大小，KLFromRecordingDays可能需要几个小时来完成分析。进度在GUI底部以红色字体显示，完成后将更改为 典型的KLFromRecording-Excel电子表格格式的天数输出如图所示。 五、5. 说明性实例开发该软件的项目涉及确定实验性新药可以在多大程度上改善斑胸草雀皮层发声前运动区的分数损伤（约10%的微损伤）的恢复这种微损伤程序最初由[6]开发，导致典型动物的发声模式暂时中断约7天。能够减轻这些病变影响的药物应减少恢复时间和/或声音模式破坏的程度。与溶剂对照组相比，药物处理动物的KL距离测量值降低，应可明显在最初的将这些剂量组进一步分为三个手术条件：无微损伤、假微损伤和微损伤组（每个药物剂量n将 动物圈养在置于记录室（302416“塑料盒，装有灯和麦克风，悬挂在顶板中心）。通过计算机连续监测麦克风，并使用声音分析记录器将符合歌曲标准的声音以 *.wav格式保存到硬盘驱动器中（见图1）。分析步骤请参见http://soundanalysispro.com/manual-1/chapter-7-recording-vocal-sounds/recorder-overview，最后访问日期：2017年9月11日）。基线记录进行了三天。每日药物治疗在外科手术前6天开始在基线后第7天进行外科手术。恢复期再随访10天。可以使用基线记录作为“模板"和随后的六个术前药物治疗日作为”目标“来评估手术程序之前药物治疗的效果可以使用在作为“模板”的六个手术前药物治疗日和作为“目标”的十个手术后日期间进行的记录来评估微损伤之后的声音模式的恢复将 *.wav格式的记录复制到USB硬盘或云存储中进行分析。使用实用程序Pars-eSAPRecorderWaves（与KLFromRecordingDays在同一GitHub存储库中提供）按记录日期将 *.wav文件解析为子目录，并使用SAP将记录的声音分割为不同的音节。通过SAP分析记录的 总 体 情 况 。 使 用 SAP 生 成 的 MySQL 表 输 入KLFromRecordingDays，计算13个非持续时间声学参数中每一个的KL距离值，并以Excel电子表格格式存储。KLFromRecordingDays使用SAP分割的所有音节，每个 *.wav文件至少有五个音节（在GUI中指定）。因此，分析是完全客观的，处理了所有记录，并且所有*.wav文件由用于KL距离计算的至少五个分段音节组成。总结了剂量和手术条件组的结果值，并通过使用Holm-Sidak事后比较的双因素ANOVA进行统计学比较初步分析表明，剂量和手术条件之间存在显著差异。在检查结果从K. Soderstrom，A.Alalawi/SoftwareX 6（2017）271277图五、 Excel电子表格格式的典型KLFromRecordingDays结果。目标记录天数从第五行开始按行列出计算KL距离的13个声学特征被写入列列P包含跨所有声学特征的平均KL距离测量列Q示出了用于该记录日的分析的单个音节的数量标题包含动物ID、GUI设置（包括模板使用的天数地区（见图）6），并应增加鸣禽作为临床前模型的价值。我们还开始应用该软件来评估药物治疗在鸣禽发育过程中的语音效果，以及在鸣禽受到约束和慢性轻度不可预测压力除了我们将该软件应用于斑胸草雀的歌声之外，它还可用于评估其他物种随时间的语音变化，这些物种的发声符合SAP分割成音节。这包括其他发声学习者，如蝙蝠和鲸目动物，以及大量产生类似叫声的动物。考虑到将人类语言分割成音节的挑战和语音的不规则性，软件在分析大多数自然产生的人类发声的能力方面可能受到限制。然而，对这些类型的人类发声进行分析仍然是可能的（包括笑声和不同音素和音素的产生）。因此，该软件可能有用的小，但潜在的重要见图6。使用KLFromRecordingDays生成的实验结果显示，当比较损伤前记录（作为模板）与损伤后记录天数（目标）时，药物对语音分歧的影响。较低的KL距离测量值表明病变诱导的变化程度较低。因此，药物治疗似乎显著减少了病变诱导的语音中断（* =与0 mg/kg对照组的显著差异，双因素ANOVA，一 组 接 受 微 损 伤 的动物（ 图 6 ， 使 用 GraphPad Prism 软 件https://www.graphpad.com/，最后一次访问日期为2017年6月27日）在恢复期第1天，100 mg/kg组和10 mg/kg组在恢复期第2、3、6、7天，所有声学参数的平均KL距离测量值均较0 mg/kg对照组显著降低（*p <0.05，Holm-Sidak事后检验）。这些初步结果表明，10和100 mg/kg剂量的实验药物均降低了损伤效应的程度（在各剂量组中未观察到梗死面积的差异，数据未显示）。至微损伤后7天，10 mg/kg药物组的KL距离测量值显著低于0 mg/kg对照组，表明歌曲音韵恢复时间缩短6. 结论KLFromRecordingDays的开发使我们能够利用SAP产生的大量音节声学数据的全部来客观测量语音随时间的变化。它的可用性使我们能够使用斑胸草雀来研究药物对声音运动脑现场学习笑声它也可能是临床实用的评估进展的语音治疗，其中不同的语音重复练习孤立。该软件目前设置为每天的分析，应该是适合的，因为是典型的鸣禽实验的基础上，随着时间的推移语音变化为了适应不同的实验设计，源代码在存储库中可用，使用Python进行修改应该是直接的。致谢这项工作得到了GW Pharmaceuticals的资助我们非常感谢OferTchernichovski领导的团队，他们创建并免费提供了我们工作所基于的Sound Analysis Pro 2011软件（http：//soundanalysispro.com/credits，last accessed 6/27/2017）.引用[1] Shriberg LD，Tomblin JB，McSweeny JL. J.演讲Lang.听. Res. JSLHR1999;42：1461-81.[2] Tchernichovski O，Nottebohm F，Ho CE，Pesaran B，Mitra PP. 动画。行为举止。2000;59：1167-76.[3] 吴W，汤普森JA，伯特伦R，约翰逊F. J. Neurosci. 方法2008;174：147 - 54.[4] 放大图片作者：DaouA，Johnson F，Wu W，Bertram R. J. Neurosci. 方法2012;210：147-60。[5] [10]杨文辉，李文辉. 神经科学杂志关闭. J. Soc（英国足球协会）Neurosci.2011;31：322-30.[6] 杨伟，王伟，王伟. 神经科学杂志关闭.神经科学学会杂志2007;27：12308-20.