E07 · Neurophysiological recordings from parietal areas of macaque brain during an instructed-delay reaching task
Diomedi, S., Vaccari, F.E., Gamberini, M., De Vitis, M., Filippini, M., & Fattori, P. (2024). Scientific Data, 11, 624. doi:10.1038/s41597-024-03479-7
这篇数据论文公开两只食蟹猴在延迟注视伸手任务中的682个顶叶单元，覆盖V6A、PEc和PE。主要贡献是带事件时间戳与元数据的HDF5资源及读取示例，而不是新提出一种群体控制机制。
核对版本：Zotero VR 集合中的本地 PDF（2026-09-18）；Zotero 条目 · 打开原文。图下页码均为该 PDF 的文件页序，从 1 开始，可能不同于期刊印刷页码。
背景
公开电生理数据要可复用，不仅需要时间戳，也需要任务、脑区、条件和事件定义。本文整理延迟伸手的顶叶单元记录，降低读取与理解成本。
研究思路
六个HDF5文件按动物和脑区组织，包含单元、条件、试次、spike时间及事件标记，并提供Python和MATLAB示例。任务时程和ISI分析是数据质量参考；不包含宽带信号意味着无法据此重新分拣spike或恢复完整LFP。
方法
两只雄性食蟹猴（Macaca fascicularis）：MonkeyS 4-6 岁、实验期体重 3.9-4.9 kg，MonkeyF 5-8 岁、3.7-4.0 kg。任务是完全黑暗中的 instructed-delay foveated 伸手：猴子按住嵌在猴椅中央的 Home Button（直径 2.5 cm，体前 4 cm），用对侧手臂向 9 个 LED 靶（直径 6 mm，绿/红双色）之一伸手。靶位排在眼水平的一个水平面板上，按眼坐标而非笛卡尔坐标设计：方向 version 为 −15°/0°/+15°，深度对应 vergence 6.9°、11.4°、17.1°，离 HB 距离 15.8-27.9 cm——这样排是为了照顾 PPC 对眼位置的强依赖，是这间实验室的传统设计。眼动用 ISCAN ETL200（100 Hz）经 45° 热镜（hot mirror）监测，4°×4° 电子窗判注视，每 session 前做 10 靶标定。记录用 5 通道 MiniMatrix 微电极系统（Thomas Recording），信号增益 10000、0.5-5 kHz 带通，Multi Spike Detector（Alpha Omega）做波形鉴别，锋电位采样 100 kHz；行为由 LabVIEW 软件控制。手术植入钛合金头柱（ML 0、AP +19 mm）和记录 chamber：MonkeyS 在 AP −13 mm、后倾 29°，MonkeyF 在 AP −25 mm、后倾 37°（图 1 上标注直径 18 mm），组织学重建确认电极道落在 V6A/PEc/PE 内。任务时序为 FREE 1000 ms → 眼动反应（RT SACC）→ 注视延迟 DELAY（1850 或 2350 ms 两档随机）→ go 信号（LED 变红）→ 外伸手（反应时超 1 s 或运动超 1 s 则 abort）→ HOLD（800 或 1200 ms 两档）→ 回程伸手 → WAIT END 200 ms；一个 run 为每靶 10 个正确 trial（共 90 trial），伪随机呈现，错误 trial 后重复同一靶。数据集共 682 个已分离神经元：V6A 285 个（S 157 + F 128）、PEc 231 个（134+97）、PE 166 个（62+104）。
主要结果
成果就是数据集本身。六个 HDF5 文件按"猴×脑区"命名（如 MonkeyS_V6A_reach9pos.h5），托管在 G-Node，总共约 500 MB。每个神经元每个 trial 存两样东西：spike 时间戳（ms）和 10 个行为事件标记的时机，全部数据统一对齐在 Move out on（伸手起始）上，其他时间全是相对值，用户可以自行重对齐到任何事件。16 类 attribute 把复用所需信息都挂进文件内部，比如想筛选"MonkeyS 在 V6Ad 记录的单元"或"某个靶的三维坐标"，读属性就行，不用回论文里翻。
技术验证做得比多数数据论文认真。epoch 时长全部对得上任务设计：FREE设定约1000 ms，实测1007±0.18 ms，装置控制时长的 DELAY（两档随机）平均 2130±249 ms，HOLD 平均 1013±200 ms；依赖猴子反应的时段也符合生理常识——saccade 反应时 159±133/137±118 ms，与猴类文献的 120-200 ms 区间一致；外伸手运动时间 388±86/378±84 ms，回程显著快于去程（222±82/252±44 ms），作者归因于重力帮忙或 HB 这个大目标好找。ISI 检验同样过关：全部 6 个数据集里短于 1 ms 的 ISI 只占 1.9%，最小 ISI 0.38 ms，说明单细胞分离质量高。
可用性方面附了 MATLAB（R2022a）和 Python（3.7.16）示例脚本：H5_raster 画 raster 图，H5_ISI 画 ISI 直方图并统计 violation 比例，改一个字符串参数就能在 unit、condition、trial 任意层级取数。作者还在一台 i7-1360P 笔记本上实测了读取速度：提取一个神经元的全部数据 MATLAB 要 1.42 s、Python 0.18 s，整个 MonkeyS_V6A 文件是 260 s 对 33 s。这些数字本身不算成果，但透露出作者确实站在第一次打开文件的用户角度想过问题。
逐图解读
来源：原文 Fig. 1，PDF 第 2 页。
图 1 · V6A、PEc、PE 的解剖定位与记录 chamber 位置。 (a) 右半球后外侧观加左半球后内侧观，品红、绿、橙虚线分别圈出 SPL 上 V6A、PEc、PE 的范围，右半球做了部分剥离以露出藏在 parieto-occipital sulcus 里的区域；实色区域是两只猴共四个半球记录范围的叠加重建，也就是说位置信息出自自家组织学而非标准图谱。(b, c) 内侧面观的三维脑重建，标出 chamber 的直径（18 mm）和倾角——MonkeyS 后倾 29°、MonkeyF 37°；(d, e) 背面观给出 chamber 在颅顶的位置，黑色圆圈是对 chamber 位置和尺寸的估计。重建用 Caret 软件基于组织学切片完成。对使用者来说，这张图的价值在于判断"我关心的区域和他们的记录位点重合度有多高"，这是决定数据能不能复用的第一步。
来源：原文 Fig. 2，PDF 第 3 页。
图 2 · 实验装置全景。 (a) 侧视示意图，配色标出各部件：青色是带热镜（hot mirror）的眼动追踪系统，绿色是奖励系统（量筒、电磁阀加喂水管），红色是记录 chamber，紫色是猴椅，橙色是带 9 个靶的面板，Home Button 嵌在猴椅中央作为每次伸手的起点。(b) 顶视图给出 9 个靶的空间坐标——以 version（方向角，−15°/0°/+15°）和 vergence（深度角，6.9°/11.4°/17.1°）表示，从近左到远右共九个位置。读这张图的要点是靶位不按等欧氏距离排，而按眼坐标组织（给出的集合深度角数值并非严格等差），这直接服务于 PPC 眼位置依赖性的研究传统，也意味着用这份数据时深度和方向这两个运动学分量的定义与通常的 center-out 范式不一样。
来源：原文 Fig. 3，PDF 第 5 页。
图 3 · 任务时序图。 上排是 10 个事件标记，从左到右：Start（按下 HB）、Green on（靶亮绿）、Fix on（开始注视）、Green to red（go 信号）、Move out on（松开 HB，伸手动身）、Move out off（触到靶，进入 HOLD）、Red off（回程 go 信号）、Move in on（松靶）、Move in off（回到 HB，同时触发奖励）、End。下排是各 epoch 及时长，两只猴并列标注：装置控制的 FREE、DELAY、HOLD 有固定范围，反应时和运动时给出实测均值，如 RT MOVE OUT 为 S 264±45 ms、F 244±42 ms。这张图是使用数据集前必读的"字典"——HDF5 里存的 event marker 就是这十个标签，epoch 的切法也以此为准。
来源：原文 Fig. 4，PDF 第 6 页。
图 4 · HDF5 文件的层级组织。 把文件结构画成一棵树：/DATA 之下是 unit_01 到 unit_N，每个 unit 下是 condition_01 到 condition_09（对应 9 个靶），每个 condition 下是 trial_01 到 trial_10，最底层是真正装数据的两个 dataset——spike trains 和 event markers。group 相当于文件夹，dataset 相当于文件，attribute 相当于贴在任意一层上的标签。条件数固定为 9，单元数各文件不同，trial 数也可能因单元中途丢失而不足——正文提到有的神经元没撑完整个 run。这张图配合正文列出的 16 类 attribute 清单，基本可以让不看代码的人明白怎么索引数据。
来源：原文 Fig. 5，PDF 第 8 页。
图 5 · 事件标记相对伸手起始的时间分布。 每一行一个事件，蓝色轮廓是该事件相对对齐点 Move out on（灰竖线，t=0）的高斯概率分布，左列 MonkeyS、右列 MonkeyF。离对齐点越远的事件分布越宽——Start 在 S 是 −3559±291 ms，这是把所有 trial 对齐到运动起始的自然代价；靠近对齐点的事件则很紧，如 go 信号 Green to red 在 S 是 −264±45 ms，运动终点 Move out off 在 +388±86 ms。两列分布总体相近（两猴各事件均值只差几毫秒到十几毫秒，如 Fix on 为 −2394±254 对 −2375±249 ms），提供跨动物时程的描述性比较。用这份数据做时间对齐分析时，这张图就是选择对齐事件的参考手册。
来源：原文 Fig. 6，PDF 第 9 页。
图 6 · 示例 raster 图（H5_raster 脚本输出）。 MonkeyS V6A 的 unit_11：横轴时间（约 −5000 到 +3000 ms）、纵轴 trial 序号，每条黑线一个动作电位，彩色圆圈是各事件标记，对应图 3 的十个事件。正文没有描述这个单元的响应特性——它的角色是示范而非展示结果：改动脚本开头的 'str'，就能画整个 unit 的全部 trial、某个条件（靶位）下的 trial，或某个具体 trial。新用户把这一个脚本跑通，基本就掌握了整个文件结构的读法。
来源：原文 Fig. 7，PDF 第 9 页。
图 7 · 示例 ISI 直方图（H5_ISI 脚本输出）。 MonkeyS V6A 全群体的 ISI 分布，横轴相邻锋电位间隔（图示 0–100 ms），纵轴计数，图上直接标出低于 1 ms 阈值的 violation 比例：0.21997%。把这个数字放回正文的全集统计——六个数据集合计只有 1.9% 的 ISI 短于 1 ms、最小 ISI 0.38 ms——可作为分离质量的参考，但不能据此排除重复单元、混合或噪声。这张图同时也是模板：使用者应对自己关心的任何子集做同样的检验。
讨论与解读边界
九个目标结合方向与深度，适合重分析任务期和目标相关响应。使用前需明确单元是否同步记录、实际试次数、脑区与动物层级，以及各时间戳的对齐定义；不能把全部682个单元直接当作同一次记录的群体。本文也没有提供连续手部三维轨迹的完整闭环数据。
事件时长核对和ISI分布是质量检查的一部分。较少的短ISI不能排除单元混合、重复分离或其他噪声；两个动物事件均值相近，也不是严格的跨动物等效性检验。图6、图7是读取与验证示例，不能被当作新的神经机制结果。论文的原始表述和参数应与下载数据分别核验，本次仅审查文献及其图。
思考问题
- 若用这些单元构造伪群体，怎样避免把跨试次平均结构误当成同步群体动力学？
- 在目标方向、深度和注视位置共变的设计中，哪些编码假说可以区分，哪些仍不可辨识？