A07 · Expanding the primate body schema in sensorimotor cortex by virtual touches of an avatar
Shokur, S., O'Doherty, J.E., Winans, J.A., Bleuler, H., Lebedev, M.A., & Nicolelis, M.A.L. (2013). Proceedings of the National Academy of Sciences. doi:10.1073/pnas.1308459110
同步的视觉—触觉配对后，猕猴 S1/M1 出现对单独虚拟触碰的反应。结果支持跨模态可塑性，但未直接证明猴子产生了橡胶手错觉或对虚拟手臂的身体所有权。
核对版本：Zotero VR 集合中的本地 PDF（2026-09-18）；Zotero 条目 · 打开原文。图下页码均为该 PDF 的文件页序，从 1 开始，可能不同于期刊印刷页码。
背景
视觉中的虚拟身体与真实触觉同步出现，可能改变感觉运动皮层对虚拟事件的响应。神经响应改变与主观身体拥有感是不同问题，需要分别检验。
研究思路
作者安排同步和不同步的视觉—触觉配对，并在配对前后测试仅有虚拟触碰时的S1/M1活动。分析比较响应强度、潜伏期和不同配对条件，以判断跨感觉关联能否在短时间内形成。
方法
两只成年恒河猴（M、N），右侧半球植入 96 通道不锈钢微丝阵列（每半球两块，分别瞄在上肢和下肢代表区），同时采样 M1 和 S1。记录位点用三重手段确认：术中皮层表面解剖标志、棉签手工探查 S1 的 receptive field（猴 M 的 RF 分布在前臂一整片、与刷擦区相当；猴 N 的 RF 小而集中在食指基部——刷擦位置就相应定在前臂和食指）、以及 intracortical microstimulation（双相刺激，100 Hz、75 μA；原文的“100 ms 脉宽”存在单位疑点，见讨论——能引发运动和 EMG 爆发的是 M1，不能的是 S1）。
装置上，猴子手臂以软绑带固定在与水平面成 23° 的扶手上、藏在屏幕下方，23° 的倾角让手臂与屏幕上的虚拟手臂平行；LCD 屏距眼 40 cm、仅比扶手高 10 cm——这些细节都是为了消除可能削弱 RHI 的姿势不协调。avatar 用 Autodesk 3Dsmax 按猴 M 的形态学建模、贴图混合了两只猴的真实毛色纹理。头部不固定，用基于 TLD/Predator 算法的眼动追踪持续监视注视，看屏幕时每 350±50 ms 随机给一次果汁奖励，注视屏幕不足 90% 时长的 trial 被剔除。16 通道 Delsys 表面 EMG 同步记录双臂的肱二头肌、三角肌、肱三头肌、腕屈伸肌，用于检查明显的肌肉活动混淆。每 trial 虚拟球刷擦 1.5 s，间隔 3.3±0.3 s；每猴在分开的日期各做两个同步 session 和一个异步 session。分析上以刺激前 −1.2 到 −0.2 s 为基线，在 PSTH 峰值处取 200 ms 窗计算 ΔR（同时用原始放电率和 z 分数），反应显著性用 Wilcoxon 检验，Vonly 与 V+P 反应的潜伏期差用 PSTH 互相关的峰值时间估计（避免对较弱的 Vonly 反应高估潜伏期）。
主要结果
真实触觉（V+P sync）在两猴共挑出 130 个显著反应的神经元（122 兴奋对 8 抑制），确认记录与刷擦刺激有效。猴 M（S1 62 个、M1 82 个神经元）中，73% 的 S1 神经元和 79% 的 M1 神经元对 V+P sync 兴奋；猴 N（S1 15 个、M1 34 个）比例低一些（S1 33%、M1 26%），可能因为食指 RF 小、只在刷擦末段扫过时才驱动。
核心发现是 Vonly 反应的"习得"。猴 M 的 E1 期间（尚未配对），Vonly 未引起明显的群体反应；猴 N 没有 E1；经过 E2 中与真实刷擦同步配对后，猴 M 有 50% 的 S1 神经元（31/62，全部兴奋）和 70% 的 M1 神经元（57/82，51 兴奋 6 抑制）开始对纯虚拟触碰显著放电，猴 N 为 40% 的 S1 和 38% 的 M1。量值上，猴 M 的 Vonly 反应约为 V+P 反应的一半：S1 从基线 34.9±2.7 sp/s 在 V+P sync 中升到 44.8±3.5（ΔR 9.9±2.7 sp/s），在 Vonly 中升到 39.0±2.9（ΔR 4.2±2.7）；M1 从 17.7±2.0 分别升到 29.6±2.9（ΔR 12.0±2.0）和 22.6±2.4（ΔR 5.5±1.9）。有趣的是猴 N 的 Vonly 反应幅度几乎追平 V+P（S1 ΔR 6.5 对 7.5；M1 8.9 对 9.2），两猴的反应幅度不同，但不能据此量化错觉强弱。EMG 检查削弱了所记录肌肉运动的解释：偶发的肌电与刺激不锁时（P>0.5），且随 session 推进反而减少，而皮层反应保持。
两个对照把结论夹得很紧。异步范式下同样的视觉和触觉刺激、只是不配准，Vonly 反应几乎不出现（猴 M 的 S1 仅 4/62、M1 仅 4/82；猴 N 各 2 个），Vonly 与 V+P async 的反应幅度之间也没有相关；同步范式下两者的幅度却高度相关（V+P 反应越强的神经元，习得的 Vonly 反应也越强，两猴 S1、M1 均 P<0.001）——反应强的细胞更愿意接纳 avatar，这是"同一群神经元在表征两个身体"的旁证。潜伏期的分析提供了通路的线索：Vonly 反应比 V+P sync 反应晚 67±3.5 ms（猴 M S1）和 55±3.8 ms（猴 M M1），猴 N 为 53±20 和 60±12 ms。考虑到对真实触觉的反应潜伏期只有 30±7 ms（S1）和 25±10 ms（M1），虚拟触碰的信号要走约 90 ms 才到 S1/M1——慢出 50–70 ms 的这段路，作者归因于从视觉皮层经后顶叶/前运动区中转的多突触通路，因为猴的视觉皮层并没有到 S1/M1 的直接投射。
逐图解读
来源：原文 Fig. 1，PDF 第 2 页。
图 1 · 实验范式：给猴子造一个 rubber hand。 (A) 装置全景：猴子手臂藏在屏幕下方，屏幕上是第一人称视角、与真臂大小位置匹配的 3D avatar 手臂，虚拟球在 avatar 上滚动，屏幕下方一台 2 自由度机器人刷同步刷真实手臂——注意刷子的轨迹与屏上球的位置严格配准。(B) 刷擦部位按每只猴的 RF 定制：猴 M 刷前臂、猴 N 刷食指基部。(C) 三种试次组合：V+P sync、Vonly、V+P async——异步对照与同步条件的刺激内容完全一致，只破坏时间配准，这是从人类 RHI 文献继承的关键设计。(D) 试次结构：E1 是 100 个纯 Vonly（配对前的基线），E2 是 900 个 V+P 试次中夹杂 100 个 Vonly（配对中的检验），同步范式在黄底、异步在灰底。装置控制了接触位置、肢体姿势和注视，有助于比较时间同步性对神经响应的影响。该实验并未直接测量动物的主观拥有感。
来源：原文 Fig. 2，PDF 第 2 页。
图 2 · 示例神经元：从无到有的虚拟触觉反应。 (A) 单 trial 分辨的 PSTH 热图，猴 M 的一个 S1、两个 M1 神经元加猴 N 的两个神经元；从上到下三行分别是 Vonly E1、Vonly E2、V+P sync E2。关键读法是纵向对比：同一批神经元，E1 里看球刷 avatar 时未见明显变化，E2 里配对建立后，不仅 V+P sync（第三行）有反应，纯虚拟触碰（第二行）也激起了清晰的放电条纹——中间那个 M1 神经元甚至被抑制。星号是 Wilcoxon 检验的放电偏离基线。(B) 猴 M 与猴 N 的平均 PSTH，粉色横线为基线：Vonly（红）的反应形状与 V+P sync（黑）相似，只是幅度打折。这张图是全文结论的单细胞版本：触觉皮层"看得见"avatar 被摸。
来源：原文 Fig. 3，PDF 第 3 页。
图 3 · 群体证据与潜伏期分析。 (A、B、F、G) 是 z 分数化的群体 PSTH 色图，每行一个神经元、按 V+P 反应从弱到强排序（上弱下强），只纳入对真实触觉有显著兴奋反应的神经元（猴 M 同步 session 51 个 S1、59 个 M1；异步 52 个 S1、64 个 M1；猴 N 异步各 6、8 个）。每个色图组自上而下三块：Vonly E1、Vonly E2、V+P——同步范式（黄底）里可以看到"排名越靠下的神经元（V+P 反应越强），E2 的 Vonly 反应越明显"，而 E1 一块整体灰暗；异步范式（灰底）里 E2 的 Vonly 同样几乎空白。(C、D、H、I) 群体平均 PSTH 及其累积和：同步条件下 Vonly（红）抬升清晰，异步条件下则平躺——累积和曲线是判断反应起点的保守工具。(E、J) 潜伏期差直方图（以原图互相关的时移定义为准，正向延迟在文中表示虚拟触碰反应较晚）：红虚线为高斯拟合，箭头标均值，猴 M S1 为 67±3.5 ms、M1 为 55±3.8 ms。这 50–70 ms 的"迟到"是推断多突态视觉-体感通线的定量依据，也是全文机制层面最有分量的一个数字。
来源：原文 Fig. 4，PDF 第 4 页。
图 4 · 幅度相关与时程量化：同步配对的作用。 (A–D) 散点图，每个点一个神经元，横轴 V+P 反应幅度（spikes/s），纵轴 Vonly 反应幅度，分别对应两猴的 S1 和 M1；同一张图里画了 Vonly E1（不相关，回归线近乎水平，P>0.1）、Vonly E2 同步（强相关，斜率显著，P<0.001）、Vonly E2 异步（不相关）三组。"配对前不相关、配对后相关、异步配对后仍不相关"这个三联模式是全文最紧凑的论证。(E–H) 柱状图汇总 Vonly 反应幅度（z 分数）在 E1 与 E2、同步（黄底）与异步（灰底）之间的对比，两猴 S1、M1 全部呈现同一格局：E1 与异步 E2 贴零，同步 E2 显著抬升（***P<0.001）。这四张子图放在一起等于把"虚拟触碰反应的出现需要时间配准的同步配对"这句话翻译成了统计语言。
讨论与解读边界
同步与异步配对的对照，以及 EMG 记录，使“配对经验改变感觉运动皮层响应”的结论较有说服力。虚拟触碰反应比真实触碰晚约 50–70 ms，与较长的视觉到感觉运动皮层通路相容，但潜伏期不能定位具体中继区域。
身体图式扩展是作者提出的解释，联结学习仍是可行的替代解释。只有猴 M 具有配对前 E1；猴 N 仅做 E2，不能把两只猴都写成完整前后对照。本文没有测量主观身体所有权，也未证明这些响应对错觉或假肢使用必需。方法中刺激脉宽原文写作 100 ms，与 100 Hz 的脉冲频率存在明显疑点；本笔记保留该疑点，不擅自改成另一单位，也不将其作为可复用刺激参数。
思考问题
- 加入虚拟手位置偏移后的定位行为测验，能否把身体图式改变与普通视觉—触觉联结学习区分开？
- 同步配对停止后，虚拟触碰响应能保持多久，是否能迁移到不同外观或位置的虚拟手？