基于猕猴Spike运动解码的不同解码方法性能对比

  1. 公开数据集中文版详细描述 参考前文:https://editor.csdn.net/md/?not_checkout=1&spm=1011.2124.3001.6192
  2. 神经元Spike信号分析 参考前文:https://blog.csdn.net/qq_43811536/article/details/134359566?spm=1001.2014.3001.5501
  3. 神经元运动调制分析 参考前文:https://blog.csdn.net/qq_43811536/article/details/134401004?spm=1001.2014.3001.5501
  4. 离散卡尔曼滤波器算法详解及重要参数(Q、R、P)基于公开数据集解码的讨论 参考前文:https://blog.csdn.net/qq_43811536/article/details/134403022?spm=1001.2014.3001.5501
  5. 基于Spike信号的运动解码分析不同运动参数对解码的影响 参考前文:https://blog.csdn.net/qq_43811536/article/details/134477651?spm=1001.2014.3001.5501

目录

  • 0. 公开数据集
  • 1. 不同解码方法的性能比较
    • 1.1 基于Spike的运动解码
    • 1.2 特征提取
    • 1.3 解码算法
    • 1.3 解码结果
    • 1.4 解码性能


0. 公开数据集

  • 网址:Nonhuman Primate Reaching with Multichannel Sensorimotor Cortex Electrophysiology 1
  • Session:
    • " indy_20170124_01 "

1. 不同解码方法的性能比较

从 公开数据集中文版详细描述 以及其 神经元Spike信号分析 中我们知道:猕猴的感觉运动皮层对运动的调制主要有六个参数——二维平面中x和y方向上的位移、速度和加速度。这六个参数也是我们在运动解码中感兴趣的部分。

1.1 基于Spike的运动解码

参见 基于Spike信号的运动解码分析不同运动参数对解码的影响 的1.1节。

1.2 特征提取

  • 预处理 s p i k e _ f e a t = [ 5900 , 253 ] spike\_feat=[5900, 253] spike_feat=[5900,253], m o t i o n _ f e a t = [ 5900 , 6 ] motion\_feat=[5900, 6] motion_feat=[5900,6]
  • Spike特征:取前后各6个样本的feat作为当前样本的联合特征,即单个样本的总特征数为: 13 ∗ 253 = 3289 13*253=3289 13253=3289维。(卡尔曼滤波器没有此步骤)
  • 数据集划分:训练集:测试集=9:1。
  • 归一化:采用均值-方差归一化。

1.3 解码算法

分别采用线性回归(Linear Regression, LR)、卡尔曼滤波器(Kalman Filter, KF)、DNN以及LSTM 2 进行解码性能比较。

  • KF:参考 离散卡尔曼滤波器算法详解

  • DNN, LSTM的网络结构如下:

    • 三层网络:输入-隐藏层-输出,隐藏层均为200个单元。
    • 激活函数:DNN为‘relu’,LSTM为‘tanh’。
    • 损失函数:均为MSE均方误差。
    • 优化器:DNN为‘Adam’,LSTM为‘RMSProp’。

最后,为了方便比较,DNN和LSTM都训练了固定的20个epochs

1.3 解码结果

Figure 1 呈现四种不同方法的部分解码结果。其中LR和KF在x、y两个不同方向的位置预测上比其他两个神经网络更精准,后者波动明显较大;但前者在速度和加速度的预测上明显弱于神经网络,后者可以捕捉到速度和加速度较大的波动,当然也正是因为这个原因导致后者预测的位置曲线出现了很多意料之外的毛刺。

在这里插入图片描述

Figure 1. 不同方法的部分解码结果. a. LR. b. KF. c. DNN. d. LSTM. 从左到右依次为十折交叉验证的第1、5、10折的前100个samples.

1.4 解码性能

Figure 2 展示了十折交叉验证的平均解码性能。Figure 2a 为Pearson相关性的柱状图,横轴为不同的方法,纵轴为相关性大小,同颜色的点代表不同折实验的性能,黑色竖线代表方差范围,反应解码方法的鲁棒性。Figure 2b 为不同方法在不同运动参数上的性能折线图,可以看到与1.3节的结论一致:LR和KF在位置预测上性能更好,而DNN和LSTM在速度和加速度的预测上更加精确,其中同色阴影为十折实验的误差范围。Figure 2c 为不同方法在不同运动参数上的 R 2 R^2 R2,结果表现与b图一致。

在这里插入图片描述

Figure 2. 不同方法的解码性能

Figure 3 绘制了DNN和LSTM在训练过程中训练集和验证集的性能曲线。可以看到DNN存在过拟合,而LSTM还处于欠拟合的状态,继续训练LSTM的性能可能会更高。

在这里插入图片描述

Figure 3. DNN与LSTM的训练过程曲线. a. MSE损失曲线. b. Pearson相关性曲线.


创作不易,麻烦点点赞和关注咯!

  1. Makin, J. G., O’Doherty, J. E., Cardoso, M. M. B. & Sabes, P. N. (2018). Superior arm-movement decoding from cortex with a new, unsupervised-learning algorithm. J Neural Eng. 15(2): 026010. doi:10.1088/1741-2552/aa9e95 ↩︎

  2. Glaser J I, enjamin A S, Chowdhury R H, et al. Machine learning for neural decoding[J]. Eneuro, 2020, 7(4). ↩︎

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/161105.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

2023 年戴森设计大奖得主是谁?给大楼降温、争取救援机会

2023 年戴森设计大奖得主是谁?给大楼降温、争取救援机会 ​编辑拉风的极客2023/11/22 摘要 当今社会除了持续不断对科技创新保持注目,还有很多年轻发明家为了实际场景的难题提供解决方案。 11 月 15 日,2023 年戴森设计大奖国际大奖名单正…

2023年危险化学品生产单位安全生产管理人员证模拟考试题库及危险化学品生产单位安全生产管理人员理论考试试题

题库来源:安全生产模拟考试一点通公众号小程序 2023年危险化学品生产单位安全生产管理人员证模拟考试题库及危险化学品生产单位安全生产管理人员理论考试试题是由安全生产模拟考试一点通提供,危险化学品生产单位安全生产管理人员证模拟考试题库是根据危…

文件的原理和应用

常识: 1 文件包括属性和内容 2 文件有打开和未打开文件, 3 本文先讨论谁打开的文件,以及如何管理已经打开的文件 一 回忆c接口 1 fopen 我们在test.c里面用一下fopen函数,不存在打开的文件会默认创建,那为什么默认新建在当前…

【Kettle实战】字符串处理及网络请求JSON格式处理

经过大量的kettle操作实践,我们会渐渐掌握一些技巧,大大减轻清洗的工作量。比如在哪里 处理字符串更方便,在哪儿处理更合理都是一个取舍问题。 字符串拼接 MySQL中使用concat(字段1,字段2),但是如果“字段2”为NULL,结…

如何在windows使用别名远程执行命令

需求背景 在开发中,需要在服务器执行脚本,需要如下几步操作: 1.打开xshell 2.登录服务器 3.进入命令脚本的路径 4.执行脚本 但是,作为懒人来说,操作太繁琐了,真麻烦,能不能一键就解决那么多操作?所以,开始研究windows有没有这个东西,而且不需要额外的软件就可以实现的.结…

ABeam Recruiting | ABeam旗下德硕管理咨询(上海)最新招聘岗位

诚聘英才 市场部经理 招聘背景 作为起源于亚洲的全球化咨询公司, ABeam立志成为中国、乃至全球的杰出品牌 本次招聘岗位 以加快ABeam中国本土各法人 在中国市场的品牌推广、提高企业知名度 强化Marketing部门的体制为目标 工作职责 在现有的PR业务基础上&#…

C++多线程学习(二):多线程通信和锁

参考引用 C11 14 17 20 多线程从原理到线程池实战代码运行环境:Visual Studio 2019 1. 多线程状态 1.1 线程状态说明 初始化 (lnit):该线程正在被创建就绪 (Ready):该线程在就绪列表中,等待 CPU 调度运行 (Running):…

xss-labs靶场6-10关

文章目录 前言一、靶场6-10关1、关卡62、关卡73、关卡84、关卡95、关卡10 总结 前言 此文章只用于学习和反思巩固xss攻击知识,禁止用于做非法攻击。注意靶场是可以练习的平台,不能随意去尚未授权的网站做渗透测试!!! …

在win10上安装pytorch-gpu版本2

安装anaconda即下载了python,还可以创建虚拟环境。 目录 1.1 anaconda安装 1.2 pytorch-gpu安装 1.1 Anaconda安装 anaconda的安装请看我之前发的tensoflow-gpu安装,里面有详细的安装过程,这里不做重复描述,传送门 1.2 pyt…

羊大师提示,羊奶都有哪些惊人功效?

羊奶不仅是一种美味的健康饮品,在近年来备受瞩目的的健康圈子里,羊奶还被赋予了更多的功效,成为一种备受推崇的保健品。羊奶不但富含营养,而且还有着非常多的益处,它能够用来美容、保健,甚至还可以治疗某些…

【Java】多线程-单例模式/volatile-指令重排序

单例模式即代码中只有一个实例的模式 适用场景:有些场景下,有的类只能有一个对象,不能有多个 要注意:在单例模式下,要保证不能产生多个实例 1、饿汉模式 class Singleton{private static Singleton instance new …

Mybatis plus 简介

简介 MyBatis-Plus (opens new window)(简称 MP)是一个 MyBatis (opens new window)的增强工具,在 MyBatis 的基础上只做增强不做改变,为简化开发、提高效率而生。 官网:https://baomidou.com/pages/24112f/ 特性 无侵入&…

《QT从基础到进阶·三十八》QWidget实现炫酷log日志打印界面

QWidget实现了log日志的打印功能,不仅可以在界面显示,还可以生成打印日志。先来看下效果,源码放在文章末尾: LogPlugin插件类管理log所有功能,它可以获取Log界面并能打印正常信息,警告信息和错误信息&…

runnergo全栈测试平台

一、全栈测试平台runnergo使用 官网 官方使用文档 二、单接口测试 三、性能测试 1.性能测试 2.性能测试报告 四、自动化测试(暂时不支持UI自动化,或许会上)

Jmeter 压测实战保姆级入门教程

1、Jmeter本地安装 1.1、下载安装 软件下载地址: https://mirrors.tuna.tsinghua.edu.cn/apache/jmeter/binaries/ 选择一个压缩包下载即可 然后解压缩后进入bin目录直接执行命令jmeter即可启动 1.2 修改语言 默认是英文的,修改中文,点击…

使用Java Servlet生成动态二维码

文章目录 引入ZXing库创建QRCodeServlet部署到Servlet容器拓展功能1. 动态生成二维码内容2. 调整二维码尺寸3. 错误修正级别4. 日志输出 结语 🎉欢迎来到Java学习路线专栏~探索Java中的静态变量与实例变量 ☆* o(≧▽≦)o *☆嗨~我是IT陈寒🍹✨博客主页&…

【追求卓越04】数据结构--栈与队列

引导 今天我们开始学习栈与队列的内容,我觉得栈并不难,所以篇幅也就不会那么多了。在虚拟空间中,栈是用户空间中的一种数据结构,它主要用于保存局部变量。那么问题来了,为什么用栈来保存局部变量,不用别的数…

Go 语言函数、参数和返回值详解

函数是一组语句,可以在程序中重复使用。函数不会在页面加载时自动执行。函数将通过调用函数来执行。 创建函数 要创建(通常称为声明)一个函数,请执行以下操作: 使用 func 关键字。指定函数的名称,后跟括…

10年开发工程师总结,8大主流程序员兼职平台,月入30k不是梦!

今年互联网行业陆续裁员减薪,许多人怨声载道的同时也开始另谋出路。而对于程序员更是应该提早做好准备,活跃在兼职接单的最前沿。 我们程序员是一门技术工种,与互联网其他行业相比薪水会相对高一点,不过钱也不是那么好赚的&#…

2023年【广东省安全员B证第四批(项目负责人)】报名考试及广东省安全员B证第四批(项目负责人)复审考试

题库来源:安全生产模拟考试一点通公众号小程序 广东省安全员B证第四批(项目负责人)报名考试是安全生产模拟考试一点通总题库中生成的一套广东省安全员B证第四批(项目负责人)复审考试,安全生产模拟考试一点…