Density/validation.m at master · nbul/Density · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
clear variables
clc
currdir = pwd;
addpath(pwd);
warning('off','stats:kmeans:FailedToConvergeRep');
cells = 10;
MTnumber = 25;
I = 25;
setup;


for Eccentricity = 0.7:0.1:0.9
    setcell;
    for I = 15:10:35
        for distribution = 30:10:40
            bundling = 0;
            counterMT = 0;
            Averages_MT = zeros(1,18);
            for MTnumber = 25:25:250
                counterMT = counterMT +1;
                for k=1:cells
                    MT;
                    threshold;
                    calcdens;
                    calcMTSD;
                end
                cd(currdir);
                %% Von Mises
                vonmises;
                SummaryMT;
            end

            summary_filename = ['MTnumber_','SD', num2str(distribution),'_int',...
                num2str(I),'_Ecc',num2str(Eccentricity),'_summary.csv'];

            cd(result_dir);
            csvwrite_with_headers(summary_filename,Averages_MT,headersMT);
            cd(currdir);
            MTnumber = 100;
            Averages_B = zeros(1,18);
            counterB = 0;
            for bundling = 0:10:40
                counterB = counterB + 1;
                for k=1:cells
                    MT;
                    threshold;
                    calcdens;
                    calcMTSD;
                end
                cd(currdir);
                %% Von Mises
                vonmises;
                SummaryB;
            end

            summary_filename = ['Bundling_','SD', num2str(distribution),'_int',...
                num2str(I),'_Ecc',num2str(Eccentricity),'_summary.csv'];
            cd(result_dir);
            csvwrite_with_headers(summary_filename,Averages_B,headersB);
            cd(currdir);
        end
    end
end
clear variables
close all
clc