问题标签 [neuroscience]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
neural-network - Keras 不使用类进行模型匹配
我是 Keras 的新手,我正在尝试建立一个神经网络来识别 38 个案例。我创建了这样一个模型,但它不起作用。我认为最后一层有问题。我检查了摘要,看起来最后一层的输出应该是 38。有人可以帮我让它工作吗?
我的代码是:
错误看起来像:
matlab - IAF 神经元模拟没有运行?[Matlab-未定义]
我想知道是否有人可以提供帮助,我正在尝试阅读以下书籍:
Cohen, Mike X. 面向大脑和认知科学家的 MATLAB(麻省理工学院出版社)(第 490 页)。麻省理工学院出版社。Kindle版。
这是为简单的集成和触发神经元建模而给出的代码示例的精确剪切和粘贴,但是每次我运行它时都会收到一个未定义的错误:(,根据以前的知识,我知道通常未定义与调用一个没有已明确声明。问题是,我不知道“ti”或“timevec”应该与什么相关?(我尝试将它们切换到我已经声明的变量,例如“时间”,但这只会导致更多问题!)也许这是我错过的语法更改并快速修复?我在 OSX 上使用 R2017b。
谢谢!
python - MNE 如何将带有事件的完整原始 EEG 信号保存为 matlab 文件?
这看起来很简单,但我迷路了。我在 edf 文件中有原始 eeg 信号,在 csv 文件中有事件。我设法通过这些事件创建纪元,然后用标记的事件绘制原始信号,但是我需要将其传输到可在 matlab 中访问的文件中。我猜我不能将 scipy.io.savemat 用于情节?我不能简单地将这些文件加载到 matlab 中,因为它以某种方式使时期不同步,可能是由于平均采样频率。MNE 不这样做,但必须在 matlab 中进行进一步分析。非常感谢您的帮助,我在 MNE 网站上找不到答案。
python-3.x - 我是否正在阅读泄漏的集成触发或突发神经元模型的重侧步函数正确性的符号?
作为一个充满激情的项目,我正在纽约大学的 XJ Wang 实验室重建一个神经元模型。论文是 Wei, W., & Wang, XJ (2016)。皮质 - 基底神经节 - 丘脑皮质环的抑制控制:复杂的调节和与记忆和决策过程的相互作用。神经元,92(5),1093-1105。
我遇到的主要问题是解释计算神经元膜电压微分的方程。他们包括一个用于基底神经节和丘脑底核细胞的爆裂神经元模型。这些区域中膜电压的微分方程包含导致爆发和强直尖峰的超极化反弹。该等式位于先前论文的第 2 页,该论文使用基本完全相同的模型。我已经链接到下面的论文,并且我还提供了指向确切段落的图像链接。 http://www.cns.nyu.edu/wanglab/publications/pdf/wei.jns2015.pdf
这是我无法阅读的方程式,不要担心 Isyn 它是来自突触的输入电流
方程取自这篇论文:https ://www.physiology.org/doi/pdf/10.1152/jn.2000.83.1.588
显然,需要对等式进行区分,以便我可以使用 numpy 运行它,但我暂时忽略它,因为这样做相对容易。所有 H 的中间学期给我带来了麻烦。据我了解,我应该运行执行以下操作的代码:
gt * h * H(V-Vh) * (V-Vt)
其中H(V-Vh)是重侧阶跃函数,V 是前一时间步长Vh = -60mV和Vt = 120mV处的膜电压。gt是纳米西门子中的电导效率常数。我认为为python解释这个的正确方法是......
gt * h * 重边(-60, 0.5)*(V-120)
但我不是 100% 确定我是否正确阅读了符号。有人可以确认我已经按预期阅读了吗?
其次, h是引起爆裂的停用术语,如 Smith et al., 2000(我链接到的第二个 pdf)第 2 页的最后一段中所述。我非常了解控制 h 演化的微分方程,但是 h 的值是多少?在史密斯等人。2000 年,作者说 h 以 20ms 的时间常数松弛为零,并以 100ms 的时间常数松弛为一。h放松的价值是什么?放松到统一意味着什么?
python - 有没有办法在无需重写整个代码的情况下将完整的 Python (psychopy) 程序实现到神经网络演示文稿中?
我找到了一个让 python 播放演示程序的选项,但不是相反。我有大量代码,将其重写为另一种语言将是一场噩梦。
macos - 将 libjpeg v8 安装到 /opt/local
如何在 OSX 上将 libjpeg 8 安装到 /opt/local?我问的原因是我想使用MNE-C 库制作大脑活动的小电影。
我执行了'mne_make_movie',导致以下错误:
已选择电影制作(来自数据)...正在扫描 /Users/user/mne/MNE-2.7.4-3420-MacOSX-x86_64/lib 的插件 找到 lqt_mjpeg.so...从模块获取编解码器信息 尝试加载 /Users /user/mne/MNE-2.7.4-3420-MacOSX-x86_64/lib/lqt_mjpeg.so ... /Users/user/mne/MNE-2.7.4-3420-MacOSX-x86_64/lib/lqt_mjpeg 的 dlopen 失败.so:dlopen(/Users/user/mne/MNE-2.7.4-3420-MacOSX-x86_64/lib/lqt_mjpeg.so,2):未加载库:/opt/local/lib/libjpeg.8.dylib 引用来自:/Users/user/mne/MNE-2.7.4-3420-MacOSX-x86_64/lib/lqt_mjpeg.so
原因:找不到图像 错误:lqt_find_video_codec 找不到编解码器!
由于所有内容都已编译,我认为我无法更改“mne_make_movie”正在寻找 libjpeg 的位置。所以我尝试使用 Mac Ports 安装 libjpeg,它只有 libjpeg-turbo 可用。但是,然后我得到另一个错误:
已选择电影制作(来自数据)...正在扫描 /Users/user/mne/MNE-2.7.4-3420-MacOSX-x86_64/lib 的插件 找到 lqt_mjpeg.so...从模块获取编解码器信息 尝试加载 /Users /user/mne/MNE-2.7.4-3420-MacOSX-x86_64/lib/lqt_mjpeg.so ... /Users/user/mne/MNE-2.7.4-3420-MacOSX-x86_64/lib/lqt_mjpeg 的 dlopen 失败.so:dlopen(/Users/user/mne/MNE-2.7.4-3420-MacOSX-x86_64/lib/lqt_mjpeg.so,2):未加载库:/opt/local/lib/libjpeg.8.dylib 引用来自:/Users/user/mne/MNE-2.7.4-3420-MacOSX-x86_64/lib/lqt_mjpeg.so
原因:库版本不兼容:lqt_mjpeg.so 需要 13.0.0 或更高版本,但 libjpeg.8.dylib 提供版本 10.0.0 错误:lqt_find_video_codec 找不到编解码器!
Mac Ports 上似乎没有更新的版本。
在 libjpeg 的官方网站上,我下载了 8d 版的源代码并按照安装说明进行操作。这会将所有内容安装到 /usr/local。我不熟悉makefile,也不知道我必须改变什么才能将所有东西正确安装到/opt/local。
那么如何将 libjpeg 8 版本 13.0.0 或更高版本安装到 /opt/local?
PS:Homebrew 似乎不适用,因为它没有安装到 /opt/local(并且不建议将其更改为 /opt/local),并且 libjpeg 9 似乎是唯一可用的版本。
r - 如何在 R 中在它们的直方图上绘制尖峰时间的指数分布?
所以我的问题是在我上一个问题之后的发展。我一直在努力将尖峰时间作为尖峰列车的光栅图。我采用了 100 的射击率并获得了 20 次试验的脉冲训练:代码是:
完成所有这些之后,我的下一个任务是获取 Inter-Spike 间隔的向量并获取它们的直方图。因为 ISI 的分布遵循指数分布,如果我用相同的数据绘制 ISI 的指数分布,它将匹配由直方图的高度制成的曲线。因此,为了首先获得间隔时间,我使用了:
然后,为了获得峰值间隔及其直方图的向量,我使用了以下命令行,
它给了我这个情节:
现在,我想要的是在直方图中绘制指数分布,以证明尖峰间隔的指数分布性质是合理的。我对使用哪些参数和使用哪个速率感到困惑。如果有人使用过 Interspike 区间绘图,请提供帮助。如果我的数据看起来不完整,我很抱歉,如果我遗漏了什么,请告诉我。
r - 如何在 R 中的多次试验中在特定时间获得随机观察点?
我正在研究 Spike Trains 和我的代码以获得这样的尖峰火车:
下面写了 20 次试验。该图像具有 5 次试验的代表性。
每个试验都有随机时间点出现的尖峰。现在我正在努力的工作是获得一个适用于所有 20 次试验的随机样本时间点,并且我想为每个试验获得由该点所在的间隔长度组成的向量。获取尖峰发生点的时间向量的代码是,
然后,您为希望查看 Interspike 区间向量的任何试验调用 ISI(i)。我想要的视觉表示是:
对于每次试验,我想获得一个向量,该向量具有该点所在的间隔长度。我也想弄清楚它的分布,但那是以后的事了。有人可以帮我弄清楚如何编写代码吗?任何帮助表示赞赏,即使它只是关于如何开始/从哪里看。
python - How to get the cluster to which a cell belongs, using AllenSDK?
On the "A Cellular Taxonomy of the Mouse Visual Cortex" page, "Explore the Data" section, the top row shows the list of cluster names of mouse visual cortex cells:
Question: Given a mouse visual cortex cell's specimen ID (e.g. 505808144), how do I find the cluster name to which that cell belongs?
I am able to get cell metadata via AllenSDK get_cell(id)
method, but the returned metadata does not seem to contain the name of the cluster.
r - 如何在R中列表的每个向量(不同长度)中附加第一个和最后一个元素,而不会使代码变慢?
我是 R 和 StackOverFlow 的新手。因此,如果我遗漏了什么,请告诉我。我正在模拟非同质过程以更好地理解神经行为。我的代码以我拥有的方式工作,例如每次20
持续2
几秒钟的试验(每个试验代表一个尖峰列车)。然后该列表SpikeTimes
为我提供了一个20
向量列表,其中每个向量对应于该特定试验中出现尖峰的时间戳。[例如。SpikeTimes[1]
看起来像这样,0.002250802 0.053934034...1.971574170 2.023090384
意味着在第一个尖峰列车中,尖峰发生在0.002250802
,0.053934034
依此类推。我不知道为什么它还会显示超出我的秒数限制的时间戳2
,但我稍后会处理]。我的代码看起来像这样——
我的问题是,对于列表中的每个向量SpikeTimes
;它给出了尖峰的时间戳,我还想包括尖峰列车的开始(即0
)和结束(即2
)。所以我想附加这个列表,让每个向量包含第一个条目 as0
和最后一个条目 as 2
。
然后我的SpikeTimes[1]
看起来0 0.002250802 0.053934034...1.971574170 2
和其他SpikeTimes[i]
看起来相似。我尝试在开头SpikeTimes <- c(0, SpikeTimes)
输入0
,但它只使列表有 21 个向量而不是 20 个向量0
作为第一个元素(我的意思是我明白为什么会发生这种情况)。我怎样才能以一种不会使我的代码变慢的方式来做到这一点?我是 R 的新手,在互联网上阅读对这个特殊问题没有帮助。我将不胜感激任何形式的投入。