fitKernel/mainScript_assembly.m at main · SysNeuroHub/fitKernel · 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

[saveServer, rootFolder] = getReady();
for aa = 1:2
    switch aa
        case 1
            animal = 'hugo'; yidx=1:3;
        case 2
            animal =  'ollie'; yidx = 3;
    end
    for yy = yidx
        switch yy
            case 1
                year = '2021';
            case 2
                year = '2022';
            case 3
                year = '2023';
        end
        [loadNames, months, dates, channels] = getMonthDateCh(animal, year, rootFolder);
        nData = numel(loadNames);
        thisdata = 1:nData;


        saveFolder = fullfile(saveServer, year,animal);%17/6/23
        % load(fullfile(saveFolder, 'assembly20241212.mat'),'param',...
        %     'id_pop','expval_tgt_pop','corr_tgt_pop','corr_tgt_rel_pop',...
        %     'latency_r_pop','avgAmp_hm_pop','p_hm_pop','spkOk_th_pop',...
        %     'spkOkTrials_pop','spkOkUCueTrials_pop','mFiringRate_pop',...
        %     'PtonsetResp_pop','errorIDs');

        for idata = thisdata

            datech = [months{idata} filesep dates{idata} filesep num2str(channels{idata})];
            thisid = [animal '/' year '/' datech];
            disp(thisid);

            saveSuffix = [animal replace(datech,filesep,'_') '_linear_rReg'];

            saveName = fullfile(saveFolder, [saveSuffix '.mat']);
            if ~exist(saveName,'file'); continue; end;
            load(saveName,'expval_tgt','corr_tgt','latencyStats','avgAmp_hm','p_hm','spkOk_th', ...
                'spkOkTrials','spkOkUCueTrials','mFiringRate','cellclassInfo','kernelInfo','ntargetTrials',...
                'corr_tgt_rel','spkNGRate','CueTrRate','nLatencyTrials_pref_success','prefDir',...
                'auc_hm', 'corr_tgt_avg','corr_tgt_avg_rel');
            if ~exist('corr_tgt_avg','var')
                corr_tgt_avg = nan*ones(7,1);
                corr_tgt_avg_rel = nan*ones(7,1);
            end

            if ~exist('mFiringRate','var') || mFiringRate < 5; continue; end

            thisDate = [months{idata} '_' dates{idata}];

            id_pop{idata} = thisid;

            expval_tgt_pop{idata} = expval_tgt;
            corr_tgt_pop{idata} = corr_tgt;
            corr_tgt_rel_pop{idata} = corr_tgt_rel;
            corr_tgt_avg_pop{idata} = corr_tgt_avg;
            corr_tgt_avg_rel_pop{idata} = corr_tgt_avg_rel;

            latency_r_pop{idata} = latencyStats.latency_r;
            latency_p_pop{idata} = latencyStats.latency_p; %15/2/25
            if isfield(latencyStats,'difflatency')
                difflatency_pop{idata} = latencyStats.difflatency; %26/3/25
            else
                difflatency_pop{idata} = nan;
            end
            nLatencyTrials_pref_success_pop{idata} = nLatencyTrials_pref_success;

            avgAmp_hm_pop{idata} = avgAmp_hm;
            p_hm_pop{idata} = p_hm;
            auc_hm_pop{idata} = auc_hm;
            spkOk_th_pop{idata} = spkOk_th;
            spkOkTrials_pop{idata} = spkOkTrials;
            spkOkUCueTrials_pop{idata} = spkOkUCueTrials;
            mFiringRate_pop{idata} = mFiringRate;
            PtonsetResp_pop{idata} = cellclassInfo.PtonsetResp;
            ntargetTrials_pop{idata} = ntargetTrials;
            errorIDs{idata} = 0;

            spkNGRate_pop{idata} = spkNGRate;%(numel(t_tr)-numel(spkOkTrials))/numel(t_tr)*100;
            CueTrRate_pop{idata} = CueTrRate;%(numel(spkOkTrials)-numel(spkOkUCueTrials))/numel(spkOkTrials)*100;
            kernel_pop(:,idata) = kernelInfo.kernel;
            tlags_pop{idata} = kernelInfo.tlags;
            prefDir_pop{idata} = prefDir;

            clear mFiringRate
        end

        save(fullfile(saveFolder, 'assembly20250207.mat'),...
            'id_pop','expval_tgt_pop','corr_tgt_pop','corr_tgt_rel_pop',...
            'latency_r_pop','latency_p_pop','avgAmp_hm_pop','p_hm_pop','spkOk_th_pop',...
            'spkOkTrials_pop','spkOkUCueTrials_pop','mFiringRate_pop',...
            'PtonsetResp_pop','errorIDs','ntargetTrials_pop','spkNGRate_pop',"CueTrRate_pop",...
            'kernel_pop','tlags_pop','nLatencyTrials_pref_success_pop','prefDir_pop', ...
            "auc_hm_pop",'difflatency_pop','corr_tgt_avg_pop','corr_tgt_avg_rel_pop');

        clear  id_pop expval_tgt_pop corr_tgt_pop corr_tgt_rel_pop ...
            latency_r_pop avgAmp_hm_pop p_hm_pop spkOk_th_pop ...
            spkOkTrials_pop spkOkUCueTrials_pop mFiringRate_pop...
            PtonsetResp_pop errorIDs ntargetTrials_pop spkNGRate_pop CueTrRate_pop...
            kernel_pop tlags_pop nLatencyTrials_pref_success_pop prefDir_pop...
            auc_hm_pop difflatency_pop corr_tgt_avg_pop corr_tgt_avg_rel_pop
    end
end