% Extracts and saves a field from each file
% Same as standard auto_despike.m, except we have
% the paramaters optimized for the case when there
% is no downsampling of the model data
fnarr{1} = 'fb0s_onlybase';
fnarr{2} = 'fb1s_onlyPSP';
fnarr{3} = 'fb2s_det_PSPs';
fnarr{4} = 'fb4s_stoch_PSPs';
fnarr{5} = 'fb5s_det_PSPs';
fnarr{6} = 'fb6s_bstoch_PSPs';
fnarr{7} = 'fb7s_bdet_PSPs';
fnarr{8} = 'fb8s_bstochPSP_m';
fnarr{9} = 'fb8s_bstochPSP_m2';
fnarr{10} = 'fb0s_balance_first_try';
fnarr{11} = 'fb3s_fix_GABAB';
fnarr{12} = 'fb19_w2f_all';
fnarr{13} = 'fb0s_onlyion';
fnarr{14} = 'fb1s_stoch';
fnarr{15} = 'fb2s_det';
fnarr{16} = 'fb3sonly_stoch';
fnarr{17} = 'fb4sonly_det';
fnarr{18} = 'fb0s_base_m50m45';
fnarr{19} = 'fb0s_base_m45m45';
fnarr{20} = 'fb0s_base_m40m40';
fnarr{21} = 'fb0s_base_m35m35';
fnarr{22} = 'acsf';
fnarr{21} = 'fb0s_ap62';
fnarr{22} = 'fb0s_ap62_det';
fnarr{23} = 'fb0s_ap62_det_soma';
fnarr{24} = 'fb0s_ap62_soma';
fnarr{25} = 'fb0s_ap92';
fnarr{26} = 'fb0s_ap92_det';
fnarr{27} = 'fb0s_ap92_det_soma';
fnarr{28} = 'fb0s_ap92_soma';
fnarr{29} = 'fb0s_bas12';
fnarr{30} = 'fb0s_bas12_det';
fnarr{31} = 'fb0s_bas12_det_soma';
fnarr{32} = 'fb0s_bas12_soma';
fnarr{91} = 'fb5s_stoch_noAMPA';
fnarr{92} = 'fb5sb_stoch_defaulttest';
datastart = 1; % AFter the PID is finished...
downsamp = 1;
plotting = 1;
% Test it on the first one
i = 5
fn = fnarr{i};
if ~(isempty(fn))
data_orig = load (fn);
data_orig_start = data_orig(1:datastart,:);
s = data_orig(datastart+1:end,:);
s = downsample(s,downsamp);
[jcs spk] = spikinator_dav3(s(:,2),10000,5,0, 0, 0.5, 0.2,2500); % Everything scaled up due to smaller dt
jcs = [ s(1:length(jcs),1) jcs ];
data = [data_orig_start; jcs];
if (plotting)
figure; plotarr([data_orig_start; s]);
hold on; plotarr(data, 'r');
% hold on; plotarr(data_orig,'g');
legend ('original','despiked');
end
channelnames = 'Vm';
save (fn, 'data', 'channelnames');
end
for i = 50:54
i
fn = fnarr{i};
if ~(isempty(fn))
data_orig = load (fn);
data_orig_start = data_orig(1:datastart,:);
s = data_orig(datastart+1:end,:);
s = downsample(s,downsamp);
[jcs spk] = spikinator_dav3(s(:,2),2000,10,0, 0, 0.5, 0.2,500);
jcs = [ s(1:length(jcs),1) jcs ];
data = [data_orig_start; jcs];
if (plotting)
figure; plotarr([data_orig_start; s]);
hold on; plotarr(data, 'r');
hold on; plotarr(data_orig,'g');
legend ('original','despiked', 'original');
end
channelnames = 'Vm';
save (fn, 'data', 'channelnames');
end
end