Revision 26:b03ef38fe497 parameterStore/MAPparamsNormal.m
| parameterStore/MAPparamsNormal.m | ||
|---|---|---|
| 1 | 1 |
function method=MAPparamsNormal ... |
| 2 |
(BFlist, sampleRate, showParams) |
|
| 2 |
(BFlist, sampleRate, showParams, paramChanges)
|
|
| 3 | 3 |
% MAPparams<> establishes a complete set of MAP parameters |
| 4 | 4 |
% Parameter file names must be of the form <MAPparams> <name> |
| 5 | 5 |
% |
| ... | ... | |
| 11 | 11 |
% output argument |
| 12 | 12 |
% method passes a miscelleny of values |
| 13 | 13 |
|
| 14 |
global inputStimulusParams OMEParams DRNLParams |
|
| 14 |
global inputStimulusParams OMEParams DRNLParams IHC_cilia_RPParams
|
|
| 15 | 15 |
global IHC_VResp_VivoParams IHCpreSynapseParams AN_IHCsynapseParams |
| 16 | 16 |
global MacGregorParams MacGregorMultiParams filteredSACFParams |
| 17 | 17 |
global experiment % used by calls from multiThreshold only |
| 18 |
global IHC_cilia_RPParams |
|
| 18 |
|
|
| 19 | 19 |
|
| 20 | 20 |
currentFile=mfilename; % i.e. the name of this mfile |
| 21 | 21 |
method.parameterSource=currentFile(10:end); % for the record |
| ... | ... | |
| 58 | 58 |
% i.e. a minimum ratio of 0.056. |
| 59 | 59 |
% 'spikes' model: AR based on brainstem spiking activity (LSR) |
| 60 | 60 |
OMEParams.rateToAttenuationFactor=0.006; % * N(all ICspikes) |
| 61 |
% OMEParams.rateToAttenuationFactor=0; % * N(all ICspikes)
|
|
| 61 |
% OMEParams.rateToAttenuationFactor=0; % * N(all ICspikes) |
|
| 62 | 62 |
|
| 63 | 63 |
% 'probability model': Ar based on AN firing probabilities (LSR) |
| 64 | 64 |
OMEParams.rateToAttenuationFactorProb=0.01;% * N(all ANrates) |
| 65 |
% OMEParams.rateToAttenuationFactorProb=0;% * N(all ANrates)
|
|
| 65 |
% OMEParams.rateToAttenuationFactorProb=0;% * N(all ANrates) |
|
| 66 | 66 |
|
| 67 | 67 |
% asymptote should be around 100-200 ms |
| 68 | 68 |
OMEParams.ARtau=.05; % AR smoothing function |
| ... | ... | |
| 104 | 104 |
DRNLParams.rateToAttenuationFactor = .01; % strength of MOC |
| 105 | 105 |
% DRNLParams.rateToAttenuationFactor = 0; % strength of MOC |
| 106 | 106 |
% 'probability' model: MOC based on AN spiking activity (HSR) |
| 107 |
DRNLParams.rateToAttenuationFactorProb = .005; % strength of MOC |
|
| 107 |
DRNLParams.rateToAttenuationFactorProb = .0055; % strength of MOC
|
|
| 108 | 108 |
% DRNLParams.rateToAttenuationFactorProb = .0; % strength of MOC |
| 109 |
DRNLParams.MOCrateThreshold =70; % spikes/s probability only |
|
| 109 |
DRNLParams.MOCrateThresholdProb =70; % spikes/s probability only
|
|
| 110 | 110 |
|
| 111 | 111 |
DRNLParams.MOCtau =.1; % smoothing for MOC |
| 112 | 112 |
|
| ... | ... | |
| 263 | 263 |
end |
| 264 | 264 |
|
| 265 | 265 |
|
| 266 |
%% now accept last minute parameter changes required by the calling program |
|
| 267 |
% paramChanges |
|
| 268 |
if nargin>3 && ~isempty(paramChanges) |
|
| 269 |
nChanges=length(paramChanges); |
|
| 270 |
for idx=1:nChanges |
|
| 271 |
eval(paramChanges{idx})
|
|
| 272 |
end |
|
| 273 |
end |
|
| 274 |
|
|
| 275 |
|
|
| 266 | 276 |
%% write all parameters to the command window |
| 267 | 277 |
% showParams is currently set at the top of htis function |
| 268 | 278 |
if showParams |
| ... | ... | |
| 276 | 286 |
eval(['UTIL_showStructureSummary(' nm{i} ', ''' nm{i} ''', 10)'])
|
| 277 | 287 |
end |
| 278 | 288 |
end |
| 279 |
end |
|
| 280 | 289 |
|
| 281 |
|
|
| 282 |
|
|
| 283 |
% ********************************************************************** comparison data |
|
| 284 |
% store individual data here for display on the multiThreshold GUI (if used) |
|
| 285 |
% the final value in each vector is an identifier (BF or duration)) |
|
| 286 |
if isstruct(experiment) |
|
| 287 |
switch experiment.paradigm |
|
| 288 |
case {'IFMC','IFMC_8ms'}
|
|
| 289 |
% based on MPa |
|
| 290 |
comparisonData=[ |
|
| 291 |
66 51 49 48 46 45 54 250; |
|
| 292 |
60 54 46 42 39 49 65 500; |
|
| 293 |
64 51 38 32 33 59 75 1000; |
|
| 294 |
59 51 36 30 41 81 93 2000; |
|
| 295 |
71 63 53 44 36 76 95 4000; |
|
| 296 |
70 64 43 35 35 66 88 6000; |
|
| 297 |
110 110 110 110 110 110 110 8000; |
|
| 298 |
]; |
|
| 299 |
if length(BFlist)==1 && ~isempty(comparisonData) |
|
| 300 |
availableFrequencies=comparisonData(:,end)'; |
|
| 301 |
findRow= find(BFlist==availableFrequencies); |
|
| 302 |
if ~isempty (findRow) |
|
| 303 |
experiment.comparisonData=comparisonData(findRow,:); |
|
| 304 |
end |
|
| 305 |
end |
|
| 306 |
|
|
| 307 |
case {'TMC','TMC_8ms'}
|
|
| 308 |
% based on MPa |
|
| 309 |
comparisonData=[ |
|
| 310 |
48 58 63 68 75 80 85 92 99 250; |
|
| 311 |
33 39 40 49 52 61 64 77 79 500; |
|
| 312 |
39 42 50 81 83 92 96 97 110 1000; |
|
| 313 |
24 26 32 37 46 51 59 71 78 2000; |
|
| 314 |
65 68 77 85 91 93 110 110 110 4000; |
|
| 315 |
20 19 26 44 80 95 96 110 110 6000; |
|
| 316 |
]; |
|
| 317 |
if length(BFlist)==1 && ~isempty(comparisonData) |
|
| 318 |
availableFrequencies=comparisonData(:,end)'; |
|
| 319 |
findRow= find(BFlist==availableFrequencies); |
|
| 320 |
if ~isempty (findRow) |
|
| 321 |
experiment.comparisonData=comparisonData(findRow,:); |
|
| 322 |
end |
|
| 323 |
end |
|
| 324 |
|
|
| 325 |
case { 'absThreshold', 'absThreshold_8'}
|
|
| 326 |
% MPa thresholds |
|
| 327 |
experiment.comparisonData=[ |
|
| 328 |
32 26 16 18 22 22 0.008; |
|
| 329 |
16 13 6 9 15 11 0.500 |
|
| 330 |
]; |
|
| 331 |
|
|
| 332 |
|
|
| 333 |
otherwise |
|
| 334 |
experiment.comparisonData=[]; |
|
| 290 |
% highlight parameter changes made locally |
|
| 291 |
if nargin>3 && ~isempty(paramChanges) |
|
| 292 |
fprintf('\n Local parameter changes:\n')
|
|
| 293 |
for i=1:length(paramChanges) |
|
| 294 |
disp(paramChanges{i})
|
|
| 295 |
end |
|
| 335 | 296 |
end |
| 336 | 297 |
end |
| 337 | 298 |
|
| 338 |
|
|
| 299 |
% for backward compatibility |
|
| 300 |
experiment.comparisonData=[]; |
|
Also available in: Unified diff