VisualTraceCheck/initialAnalyses.m at main · NERD-CO/VisualTraceCheck · 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
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
% Main Directory

% Raw isolated recordings

% Load the data

mainDir = 'Z:\Crane_Summer2020\CaseData\SpikeProcessing\SpikeDates2Cluster';
allEphysDir = 'Z:\Crane_Summer2020\CaseData\SpikeProcessing\';
summaryDir = 'Z:\Crane_Summer2020\SummaryData\';
cd(mainDir);
matFiles = dir('*.mat');
matFiles2 = {matFiles.name};


%% FIRST ANALYSIS - Fraction of active recordings

% May further follow up with specific analysis of implanted track

allFracAct = zeros(length(matFiles2),1);

for mi = 1:length(matFiles2)

    tmpMfile = matFiles2{mi};

    load(tmpMfile)

    % Get number of recording depths
    tmpDepthIND = DepthIND;
    % Get number of electrodes
    numElecs = numel(unique(fedTable.ElNum));

    totalNumRec = length(tmpDepthIND) * numElecs;

    activeNum = height(fedTable);

    fracActive = activeNum/totalNumRec;

    allFracAct(mi) = fracActive;

end
%%
cd(summaryDir)
save('ActiveRecFraction.mat','allFracAct')

%% Second ANALYSIS - Analysis the normalized root mean square of the
%%% voltage

for mi2 = 1:length(matFiles2)

    % Change directory back to surgery directory
    cd(mainDir);
    % Get a temporay file name
    tmpMfile = matFiles2{mi2};

    % Get file information
    tmpStrgs = strsplit(tmpMfile,'_');

    dateName = tmpStrgs{1};
    setNUM = str2double(tmpStrgs{3});
    surgNum = str2double(tmpStrgs{2}(1));

    % Load in surgery file for depth information
    load(tmpMfile)

    dataTable = fedTable;

    % Load depth variable
    tmpDepthIND = DepthIND;
    tmpDepthID = DepthID;

    if setNUM == 0
        addFold = '';
    elseif setNUM == 1
        addFold = 'Set1';
    elseif setNUM == 2
        addFold = 'Set2';
    end

    if surgNum == 0
        surgAdd = '';
    else
        surgAdd = ['_',tmpStrgs{2}(1)];
    end

    newDname = [dateName(1:2), '_', dateName(3:4), '_', dateName(5:end), surgAdd];

    tmpLoc = [allEphysDir  , 'RAW_Ephys_Files' , filesep, newDname , filesep , addFold];

    cd(tmpLoc)

    % Loop through our active depths
    for di = 1:height(fedTable)

        load(fedTable.Depth{di})

        % list of depth file
%         tmpIND = tmpDepthIND{di}; % this throws error

        rVolt = eval(fedTable.ElNum{di});

        oRMs = rms(rVolt);

        dataTable.RMS(di) = oRMs;

    end

    cd([summaryDir,'RMS_summary'])

    saveName = [strjoin(tmpStrgs(1:3),'_'),'.mat'];

    save(saveName,'dataTable')


end