Dynamic population coding of kinematic structure across executed and observed actions in primate premotor cortex.
The 3 matches · 2 of them tie a paragraph to a whole file, not to given lines: weak matches, whose lines are not tinted
- [1] § RESULTS › Neural population encoding of grasp configurations ↔ Fig_3.m, lines 1–122 · score 0.60 · consecutive bins, neurons selective, selectivity stability, events, Figure 3, Population
- [2] § MATERIALS AND METHODS › Canonical correlation analysis ↔ Fig_4.m, the whole file · a weak match · score 0.57 · canonical correlations, canonical dimension, shuffling, trajectories
- [3] § RESULTS › Encoding of action kinematics ↔ Fig_9.m, the whole file · a weak match · score 0.56 · cross predicted kinematics, Human kinematics predicted, Monkey kinematics predicted, model, Figure 9, contours
Paper
Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
MATLAB · 261 lines · 11 KB · CC-BY-4.0 · 1 match
- clear
- % Load data
- data_dir = 'data\neuronal_selectivity.mat';
- load(data_dir)
- names_to_plot = who('-file', data_dir);
- load('data\event_times.mat');
- graymap = [linspace(0,0.9,100)', linspace(0,0.9,100)', linspace(0,0.9,100)'; % colormap gray
- 1, 1, 1];
- cmap = [linspace(0, 0.3, 26)', linspace(0.2, 0.5, 26)', linspace(0, 0, 26)'; % colormap green-yellow-red
- linspace(0.3, 0.95, 12)' linspace(0.5, 0.9, 12)' linspace(0, 0.1, 12)';
- linspace(0.95, 1, 30)' linspace(0.9, 0.5, 30)' linspace(0.1, 0, 30)';
- linspace(1, 0.9, 35)' linspace(0.5, 0, 35)' linspace(0, 0, 35)'];
- time_points = -550:25:3350;
- time_0 = find(time_points==0);
- % Run for different conditions and populations
- for i = 1:length(names_to_plot)
- % Sort normalized activity
- eval(['curr_result =' names_to_plot{i} ';']);
- isSelective = curr_result.selectivity_p < 0.05;
- [isSelective,order] = keep_only_consecutive_bins(isSelective);
- curr_mat_to_plot = curr_result.selectivity_p .* isSelective;
- curr_mat_to_plot(curr_mat_to_plot==0) = 1;
- percent_selective_units = sum(isSelective, 1)./size(isSelective,1);
- % Plots
- f = figure();
- pos = get(f, 'Position');
- set(f, 'Position', [0, pos(2), 3*pos(3), pos(4)]);
- sgtitle(names_to_plot{i}, 'Interpreter', 'none')
- subplot(1, 3, 1); % per neuron selectivity
- hh=imagesc(curr_mat_to_plot(order,:));
- set(hh, 'AlphaData', ~isnan(curr_mat_to_plot))
- caxis([0 0.05])
- c = colorbar();
- colormap(graymap)
- xlabel('Time (s)')
- ylabel('# Neurons')
- xticks = time_0:20:length(time_points);
- set(gca, 'XTick', xticks);
- set(gca, 'XTickLabel', time_points(xticks)/1000);
- set(gca, 'tickdir', 'out')
- c.Label.String = 'selectivity p-value';
- line([find(time_points==0), find(time_points==0)],get(gca, 'YLim'),'color', [.5 .5 .5],'LineWidth',1, 'LineStyle', '--')
- line([find(time_points==1250), find(time_points==1250)],get(gca, 'YLim'),'color', [.5 .5 .5],'LineWidth',1, 'LineStyle', '--')
- line([find(time_points==800), find(time_points==800)],get(gca, 'YLim'),'color', [.9 .9 .9],'LineWidth',3)
- [~,endMoveIndex] = min(abs(time_points-endMove(i)));
- line([endMoveIndex, endMoveIndex], get(gca, 'YLim'),'color',[.5 .5 .5], 'LineStyle', '--' )
- [~,endHoldIndex] = min(abs(time_points-endHold(i)));
- line([endHoldIndex, endHoldIndex], get(gca, 'YLim'),'color',[.5 .5 .5], 'LineStyle', '--' )
- subplot(1, 3, 2); % selective units per time bin
- plot(time_points/1000, imgaussfilt(percent_selective_units,2), 'color', 'k','LineWidth',2)
- xlim([-0.600 3.400])
- ylim(get(gca, 'YLim'))
- line([0, 0],get(gca, 'YLim'),'color', [.5 .5 .5],'LineWidth',1, 'LineStyle', '--')
- line([1.250, 1.250],get(gca, 'YLim'),'color', [.5 .5 .5],'LineWidth',1, 'LineStyle', '--')
- line([0.800, 0.800],get(gca, 'YLim'),'color', [.9 .9 .9],'LineWidth',3)
- line([endMove(i)/1000, endMove(i)/1000], get(gca, 'YLim'),'color',[.5 .5 .5], 'LineStyle', '--' )
- line([endHold(i)/1000, endHold(i)/1000], get(gca, 'YLim'),'color',[.5 .5 .5], 'LineStyle', '--' )
- xlabel('Time (s)')
- ylabel('selective neurons (%)')
- ax3 = subplot(1, 3, 3); % selectivity stability map
- imagesc(curr_result.stability_map ./ (diag(curr_result.stability_map)+12))
- hold on;
- line([find(time_points==0), find(time_points==0)], get(gca, 'YLim'), 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line(get(gca, 'XLim'), [find(time_points==0), find(time_points==0)], 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line([find(time_points==1250), find(time_points==1250)], get(gca, 'YLim'), 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line(get(gca, 'XLim'), [find(time_points==1250), find(time_points==1250)], 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line([find(time_points==800), find(time_points==800)], get(gca, 'YLim'), 'color', [1 1 1],'LineWidth',3)
- line(get(gca, 'XLim'), [find(time_points==800), find(time_points==800)], 'color', [1 1 1],'LineWidth',3)
- [~,endMoveIndex] = min(abs(time_points-endMove(i)));
- line([endMoveIndex, endMoveIndex], get(gca, 'YLim'), 'color',[1 1 1], 'LineStyle', '--' )
- line(get(gca, 'XLim'), [endMoveIndex, endMoveIndex], 'color',[1 1 1], 'LineStyle', '--' )
- [~,endHoldIndex] = min(abs(time_points-endHold(i)));
- line([endHoldIndex, endHoldIndex], get(gca, 'YLim'), 'color',[1 1 1], 'LineStyle', '--' )
- line(get(gca, 'XLim'), [endHoldIndex, endHoldIndex], 'color',[1 1 1], 'LineStyle', '--' )
- set(gca,'YDir','normal')
- axis image
- ylabel('reference time (s)')
- xlabel('test time (s)')
- ticks = time_0:20:length(time_points);
- set(gca, 'XTick', ticks);
- set(gca, 'XTickLabel', time_points(ticks)/1000);
- set(gca, 'YTick', ticks);
- set(gca, 'YTickLabel', time_points(ticks)/1000);
- set(gca, 'tickdir', 'out')
- box on
- colormap(ax3, cmap)
- caxis([0 0.8])
- cc = colorbar(ax3);
- cc.Label.String = 'selectivity stability index (population average)';
- end
- %% plot congruency map
- data_dir = 'data\selectivity_congruence.mat';
- load(data_dir)
- f = figure();
- pos = get(f, 'Position');
- set(f, 'Position', [0, pos(2), 2*pos(3), pos(4)]);
- subplot(1, 2, 1); % sum of common selectivity periods
- imagesc(common_selectivity)
- hold on;
- line([find(time_points==0), find(time_points==0)], get(gca, 'YLim'), 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line(get(gca, 'XLim'), [find(time_points==0), find(time_points==0)], 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line([find(time_points==1250), find(time_points==1250)], get(gca, 'YLim'), 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line(get(gca, 'XLim'), [find(time_points==1250), find(time_points==1250)], 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line([find(time_points==800), find(time_points==800)], get(gca, 'YLim'), 'color', [1 1 1],'LineWidth',3)
- line(get(gca, 'XLim'), [find(time_points==800), find(time_points==800)], 'color', [1 1 1],'LineWidth',3)
- [~,endMoveIndex] = min(abs(time_points-endMove(1)));
- line([endMoveIndex, endMoveIndex], get(gca, 'YLim'), 'color',[1 1 1], 'LineStyle', '--' )
- [~,endMoveIndex] = min(abs(time_points-endMove(2)));
- line(get(gca, 'XLim'), [endMoveIndex, endMoveIndex], 'color',[1 1 1], 'LineStyle', '--' )
- [~,endHoldIndex] = min(abs(time_points-endHold(1)));
- line([endHoldIndex, endHoldIndex], get(gca, 'YLim'), 'color',[1 1 1], 'LineStyle', '--' )
- [~,endHoldIndex] = min(abs(time_points-endHold(2)));
- line(get(gca, 'XLim'), [endHoldIndex, endHoldIndex], 'color',[1 1 1], 'LineStyle', '--' )
- ylabel('time (s) - GO')
- xlabel('time (s) - GE')
- title('common selectivity')
- set(gca,'YDir','normal')
- axis image
- ticks = time_0:20:length(time_points);
- set(gca, 'XTick', ticks);
- set(gca, 'XTickLabel', time_points(ticks)/1000);
- set(gca, 'YTick', ticks);
- set(gca, 'YTickLabel', time_points(ticks)/1000);
- set(gca, 'tickdir', 'out')
- box on
- subplot(1, 2, 2); % sum of congruent selectivity periods
- imagesc(congruence)
- hold on;
- line([find(time_points==0), find(time_points==0)], get(gca, 'YLim'), 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line(get(gca, 'XLim'), [find(time_points==0), find(time_points==0)], 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line([find(time_points==1250), find(time_points==1250)], get(gca, 'YLim'), 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line(get(gca, 'XLim'), [find(time_points==1250), find(time_points==1250)], 'color', [1 1 1],'LineWidth',1, 'LineStyle', '--')
- line([find(time_points==800), find(time_points==800)], get(gca, 'YLim'), 'color', [1 1 1],'LineWidth',3)
- line(get(gca, 'XLim'), [find(time_points==800), find(time_points==800)], 'color', [1 1 1],'LineWidth',3)
- [~,endMoveIndex] = min(abs(time_points-endMove(1)));
- line([endMoveIndex, endMoveIndex], get(gca, 'YLim'), 'color',[1 1 1], 'LineStyle', '--' )
- [~,endMoveIndex] = min(abs(time_points-endMove(2)));
- line(get(gca, 'XLim'), [endMoveIndex, endMoveIndex], 'color',[1 1 1], 'LineStyle', '--' )
- [~,endHoldIndex] = min(abs(time_points-endHold(1)));
- line([endHoldIndex, endHoldIndex], get(gca, 'YLim'), 'color',[1 1 1], 'LineStyle', '--' )
- [~,endHoldIndex] = min(abs(time_points-endHold(2)));
- line(get(gca, 'XLim'), [endHoldIndex, endHoldIndex], 'color',[1 1 1], 'LineStyle', '--' )
- ylabel('reference time (s) - GO')
- xlabel('test time (s) - GE')
- title('congruent selectivity')
- set(gca,'YDir','normal')
- axis image
- ticks = time_0:20:length(time_points);
- set(gca, 'XTick', ticks);
- set(gca, 'XTickLabel', time_points(ticks)/1000);
- set(gca, 'YTick', ticks);
- set(gca, 'YTickLabel', time_points(ticks)/1000);
- set(gca, 'tickdir', 'out')
- box on
- c = colorbar();
- c.Label.String = 'sum of common/congruent selectivity periods';
- colormap(cmap)
- caxis([0 35])
- %%
- function [isSig, order] = keep_only_consecutive_bins(isSig)
- first_significance = nan(1,size(isSig,1));
- for iUnit= 1:size(isSig,1)
- %%%% keep only consecutive bin significant periods
- bin = 1;
- min_num_of_bins = 7; % choose minimum number of significant bins
- while bin <= size(isSig,2)
- if isSig(iUnit,bin) == 1
- cond = 1;
- k = 0;
- while cond % check the next bins
- k = k+1;
- if bin+k>size(isSig,2) % dont let the bin go above limit
- if k<min_num_of_bins
- isSig(iUnit,bin:bin+k-1) = 0;
- end
- break
- end
- if isSig(iUnit,bin+k) == 0 % find the end of the significant period
- cond = 0;
- if k<min_num_of_bins
- isSig(iUnit,bin:bin+k) = 0; % delete short significant periods
- end
- end
- end
- bin = bin+k+1;
- else
- bin = bin+1;
- end
- end % for all bins
- %%% find first significance after movement preparation
- if max(isSig(iUnit,60:157))==1
- first_index = 59+find(isSig(iUnit,60:157),1);
- % [max_pev, max_index] = max(omegaSq{iUnit}{iCond}.*significant{iUnit}{iCond});
- % first_significance(iUnit,iCond) = max_index;
- ind = 0;
- while isnan(first_significance(iUnit))
- ind = ind+1;
- if isSig(iUnit,first_index-ind)==0 && isSig(iUnit,first_index-ind-1)==0 && isSig(iUnit,first_index-ind-2)==0 && isSig(iUnit,first_index-ind-3)==0
- first_significance(iUnit) = first_index-ind+1;
- elseif first_index-ind-3==1
- first_significance(iUnit) = first_index-ind+1;
- end
- end
- end
- end
- [~,order] = sort(first_significance);
- end
Fig_3.m, under CC-BY-4.0 · at the source
Overview
- Lab of Movement Physiology, Medical School, University of Crete, Heraklion, Greece
- Graduate Program in the Brain and Mind Sciences, University of Crete, Heraklion, Greece
- Institute of Applied and Computational Mathematics, Foundation for Research and Technology-Hellas, Heraklion, Greece
Abstract
Neurons active during both action execution and observation [mirror neurons (MirNs)] are central to theories of action understanding, yet what they represent remains debated: abstract goals, static grips, or movement kinematics? We recorded 433 neurons from macaque premotor cortex during execution and observation of reach-to-grasp actions. Population analyses revealed grasp-specific information in both conditions, broadly distributed across neurons and dynamically reconfigured over time. Generalization across task phases was limited, indicating time-specific and evolving population codes rather than static representations. Execution and observation were linked by a shared, partially overlapping population geometry that supported reliable cross-condition classification, with the strongest alignment emerging during movement and hold. Neural activity was systematically related to multidimensional hand kinematics, and these relationships generalized across neuronal populations and across agents. Together, these findings support a dynamic population-level account of MirN function in which premotor circuits integrate visual and motor signals to represent and anticipate the unfolding structure of others’ actions.
Reproduced under the paper's license (CC BY), from the paper cited above.
Repository
Its files are read in the Code ↔ Paper reader above, with 3 matches between paragraphs and lines of code.
Zenodo 18983284
Availability: 1 check, the latest on 27 September 2026: the link answers (HTTP 200)
- 27 September 2026: the link answers (HTTP 200)
9 files
The paper's code and data availability statement is in the Data section.
Tracing map
Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.
What the map holds:
- 1 repository of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 9 scripts, each with its path and the digest of its content;
- 3 matches between paragraphs of the paper and lines of the code (method lexical-v1);
- neither the text of the paper nor the code itself.
Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.
Data
No dataset and no data link were found in the paper.
Data, code, and materials availability
All data and code needed to evaluate and reproduce the results in the paper are present in the paper and/
Reproduced under the paper's license (CC BY), from the paper cited above.
Versions
The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.
Version 2, 28 September 2026
- Funding: added European Social Fund: MIS 5048512
Version 1, 27 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 5 authors, 8 MeSH terms, 86 references.
Cite
This paper
Chatzimichail, K., Paschalidis, C., Tzamali, E., Papadourakis, V., & Raos, V. (2026). Dynamic population coding of kinematic structure across executed and observed actions in primate premotor cortex. Science advances, 12(25), eaed9309. https://
BibTeX
@article{chatzimichail20
author = {Chatzimichail, Konstantinos and Paschalidis, Christos and Tzamali, Eleftheria and Papadourakis, Vassilis and Raos, Vassilis},
title = {{Dynamic population coding of kinematic structure across executed and observed actions in primate premotor cortex}},
journal = {Science advances},
year = {2026},
month = jun,
volume = {12},
number = {25},
pages = {eaed9309},
publisher = {American Association for the Advancement of Science},
issn = {2375-2548},
doi = {10.1126/
url = {https://
pmid = {42319945},
pmcid = {PMC13281793}
}
RIS
TY - JOUR
AU - Chatzimichail, Konstantinos
AU - Paschalidis, Christos
AU - Tzamali, Eleftheria
AU - Papadourakis, Vassilis
AU - Raos, Vassilis
TI - Dynamic population coding of kinematic structure across executed and observed actions in primate premotor cortex
T2 - Science advances
J2 - Sci Adv
PY - 2026
DA - 2026/
VL - 12
IS - 25
SP - eaed9309
SN - 2375-2548
PB - American Association for the Advancement of Science
DO - 10.1126/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1126/
"type": "article-journal",
"title": "Dynamic population coding of kinematic structure across executed and observed actions in primate premotor cortex",
"container-title": "Science advances",
"author": [
{
"family": "Chatzimichail",
"given": "Konstantinos"
},
{
"family": "Paschalidis",
"given": "Christos"
},
{
"family": "Tzamali",
"given": "Eleftheria"
},
{
"family": "Papadourakis",
"given": "Vassilis"
},
{
"family": "Raos",
"given": "Vassilis"
}
],
"container-title-short":
"volume": "12",
"issue": "25",
"page": "eaed9309",
"DOI": "10.1126/
"PMID": "42319945",
"PMCID": "PMC13281793",
"ISSN": "2375-2548",
"publisher": "American Association for the Advancement of Science",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
6,
19
]
]
}
}
The tracing map gets a citation of its own once an author has validated it and it has a DOI.
Similar papers
The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.
- [1] doi:10.1016/j.isci.2026.116776 [code]
- Selective perturbation of mirror and non-mirror neurons in an in silico model of the action observation network.Journal: iScienceIn common: 12 references
- [2] doi:10.1126/sciadv.adw3876 [code]
- Intracortical brain-computer interface for navigation in virtual reality in macaque monkeys.Journal: Science advancesIn common: non-human primate, systems, 5 references
- [3] doi:10.1126/sciadv.aea1927 [code]
- An evolutionary conserved neural mechanism for interpersonal coordination in primates.Journal: Science advancesIn common: Statistics and Machine Learning Toolbox, non-human primate, systems, 4 references
- [4] doi:10.1038/s41598-026-55225-1 [code]
- Benchmarking criteria to determine latent linear dimensionality in neural data.Journal: Scientific reportsIn common: Statistics and Machine Learning Toolbox, 3 references
- [5] doi:10.1007/s00429-026-03111-x [code]
- Direction selectivity in naturalistic action observation: distributed representations across the action observation network.Journal: Brain structure & functionIn common: Image Processing Toolbox, Statistics and Machine Learning Toolbox, 2 references
- [6] doi:10.1016/j.neuroimage.2026.122004
- Preventing visual feedback of the moving limb during grasping preferentially activates parietofrontal premotor areas.Journal: NeuroImageIn common: systems, 3 references
- [7] doi:10.1113/jp288070 [code]
- History bias and its perturbation of the stimulus representation in the macaque prefrontal cortex.Journal: The Journal of physiologyIn common: non-human primate, 3 references
- [8] doi:10.1162/imag.a.1199 [code]
- Sustained alpha oscillations serve attentional prioritization in working memory, not maintenance.Journal: Imaging neuroscience (Cambridge, Mass.)In common: Image Processing Toolbox, Statistics and Machine Learning Toolbox, 2 references
- [9] doi:10.1038/s41586-026-10653-x [code]
- A mosaic of whole-body representations on the human precentral gyrus.Journal: NatureIn common: Statistics and Machine Learning Toolbox, systems, 2 references
- [10] doi:10.1038/s41467-026-73037-9 [code]
- Sensorimotor transformation of number in the primate parietal cortex.Journal: Nature communicationsIn common: Statistics and Machine Learning Toolbox, non-human primate, 2 references
Contribute
The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.
Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 1 repository of the authors' code, each at its verified commit and with its license, 9 scripts, and 3 matches between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:29175aa9704e485f…
Add the badge to its README
The badge links the code to this page. Copy one of these into the README of the paper's code: only you decide where it goes, and nothing is changed for you.
Markdown
[.
Discussion, reproductions, activity
Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.
Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
