teaching/scientificComputing
Archived
3
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

improved plots in population vector solution

Browse Source
This commit is contained in:
Jan Benda 2018-01-30 11:18:23 +01:00
parent c2e37b516e
commit 2103216caf
1 changed files with 7 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
projects/project_populationvector/solution/populationvector.m
View File

@ -1,6 +1,5 @@
close all close all
datapath = '../'; datapath = '../data/';
% datapath = '../code/';
files = dir(strcat(datapath, 'unit*.mat')); files = dir(strcat(datapath, 'unit*.mat'));
for file = files' for file = files'
a = load(strcat(datapath, file.name)); a = load(strcat(datapath, file.name));
@ -82,10 +81,10 @@ ylabel('firing rate')
hold off; hold off;
subplot(1, 3, 2); subplot(1, 3, 2);
hist(angleestimates1); hist(angleestimates1);
xlabel('stimulus angle') xlabel('population vector angle')
subplot(1, 3, 3); subplot(1, 3, 3);
hist(angleestimates2); hist(angleestimates2);
xlabel('stimulus angle') xlabel('max. rate angle')
angle angle
mean(angleestimates1) mean(angleestimates1)
mean(angleestimates2) mean(angleestimates2)
@ -124,9 +123,9 @@ end
figure(); figure();
subplot(1, 2, 1); subplot(1, 2, 1);
scatter(angles, e1m); scatter(angles, e1m);
xlabel('stimuluis angle') xlabel('stimulus angle')
ylabel('estimated angle') ylabel('estimated angle (population vector)')
subplot(1, 2, 2); subplot(1, 2, 2);
scatter(angles, e2m); scatter(angles, e2m);
xlabel('stimuluis angle') xlabel('stimulus angle')
ylabel('estimated angle') ylabel('estimated angle (maximum firing rate)')