annotate multithreshold 1.46/printReport.m @ 38:c2204b18f4a2 tip

End nov big change
author Ray Meddis <rmeddis@essex.ac.uk>
date Mon, 28 Nov 2011 13:34:28 +0000
parents 161913b595ae
children
rev   line source
rmeddis@0 1 function printReport(fileName, printTracks)
rmeddis@0 2
rmeddis@0 3 % End of run report (no args)
rmeddis@0 4 % *or*
rmeddis@0 5 % reprint previous report from file
rmeddis@0 6
rmeddis@32 7 global experiment stimulusParameters betweenRuns withinRuns statsModel
rmeddis@0 8 global LevittControl expGUIhandles
rmeddis@29 9 global paramChanges
rmeddis@0 10
rmeddis@32 11 % NB all globals are used even though this is not obvious
rmeddis@0 12 global inputStimulusParams OMEParams DRNLParams
rmeddis@0 13 global IHC_VResp_VivoParams IHCpreSynapseParams AN_IHCsynapseParams
rmeddis@0 14 global MacGregorParams MacGregorMultiParams filteredSACFParams
rmeddis@0 15 global experiment % used by calls from multiThreshold only
rmeddis@0 16 global IHC_cilia_RPParams
rmeddis@0 17
rmeddis@0 18 printReportGuide.structures=1;
rmeddis@0 19 printReportGuide.showPsychometric=0;
rmeddis@0 20 printReportGuide.HorizontalTracks=1;
rmeddis@0 21
rmeddis@38 22
rmeddis@0 23 if nargin==0
rmeddis@0 24 % print new report
rmeddis@0 25 printReportGuide.showTracks=experiment.printTracks;
rmeddis@0 26 printReportGuide.fileName=[];
rmeddis@33 27 % save this data (just in case)
rmeddis@33 28 saveFileName=['savedData/mostRecentResults'];
rmeddis@0 29 experiment.minElapsed=etime(clock, betweenRuns.timeNow)/60;
rmeddis@20 30 save(saveFileName, 'experiment', 'stimulusParameters',...
rmeddis@20 31 'betweenRuns', 'withinRuns', 'statsModel', 'expGUIhandles')
rmeddis@0 32 disp(['data saved as: ' saveFileName]);
rmeddis@32 33
rmeddis@0 34 else
rmeddis@20 35 % reprint request (i.e print out old data)
rmeddis@0 36 printReportGuide.fileName=fileName;
rmeddis@0 37 load(printReportGuide.fileName)
rmeddis@0 38 saveFileName=printReportGuide.fileName;
rmeddis@0 39 if nargin>1
rmeddis@0 40 printReportGuide.showTracks=printTracks;
rmeddis@0 41 else
rmeddis@0 42 printReportGuide.showTracks=experiment.printTracks;
rmeddis@0 43 end
rmeddis@0 44 end
rmeddis@0 45
rmeddis@0 46 fprintf('******** multithreshold ')
rmeddis@20 47
rmeddis@0 48 x=pwd; disp(['version ' x(end-3:end)])
rmeddis@0 49 fprintf('\nName:\t%s ', experiment.name)
rmeddis@0 50 fprintf('\nparadigm:\t%s ', experiment.paradigm)
rmeddis@0 51 fprintf('\nEar:\t%s ', experiment.ear)
rmeddis@20 52
rmeddis@0 53 method=experiment.threshEstMethod;
rmeddis@0 54 if stimulusParameters.includeCue && ~strcmp(method(1:6),'2I2AFC')
rmeddis@0 55 method=[method '/ withCue'];
rmeddis@0 56 end
rmeddis@0 57 fprintf('\nmethod:\t%s ', method)
rmeddis@0 58 fprintf('\ndate:\t%s ', experiment.date)
rmeddis@0 59 fprintf('\n\n')
rmeddis@0 60
rmeddis@0 61 if isempty(betweenRuns.thresholds)
rmeddis@0 62 disp('no thresholds found')
rmeddis@0 63 end
rmeddis@0 64
rmeddis@0 65 % prepare results as matrices ready to print
rmeddis@20 66 % first sort the actual sequence into a more readable sequence
rmeddis@0 67 [idx1, idx2, var1values, var2values]=...
rmeddis@20 68 sortVariables(betweenRuns.variableList1, betweenRuns.variableList2,...
rmeddis@20 69 betweenRuns.var1Sequence, betweenRuns.var2Sequence);
rmeddis@0 70
rmeddis@0 71 header1=betweenRuns.variableName1;
rmeddis@0 72 header2=betweenRuns.variableName2;
rmeddis@0 73 header1 = strrep(header1, 'none', ' '); % none is not a useful header
rmeddis@0 74 header2 = strrep(header2, 'none', ' '); % none is not a useful header
rmeddis@0 75 headers=strvcat([header1 '/'], header2);
rmeddis@0 76
rmeddis@0 77 disp('thresholds')
rmeddis@28 78 global resultsTable
rmeddis@28 79 resultsTable= sortTablesForPrinting(idx1,idx2,...
rmeddis@28 80 var1values,var2values, betweenRuns.thresholds);
rmeddis@28 81 msg=printTabTable(resultsTable, headers);
rmeddis@0 82 addToMsg(msg,0)
rmeddis@32 83 if ~isempty(paramChanges)
rmeddis@32 84 fprintf('\n')
rmeddis@33 85 fprintf('%s\n', char(paramChanges))
rmeddis@32 86 end
rmeddis@0 87
rmeddis@0 88 % sort tracks into the same order
rmeddis@0 89 betweenRuns.levelTracks=betweenRuns.levelTracks(idx1);
rmeddis@0 90 betweenRuns.responseTracks=betweenRuns.responseTracks(idx1);
rmeddis@0 91 betweenRuns.bestThresholdTracks=betweenRuns.bestThresholdTracks(idx1);
rmeddis@0 92 betweenRuns.levelTracks=betweenRuns.levelTracks(idx2);
rmeddis@0 93 betweenRuns.responseTracks=betweenRuns.responseTracks(idx2);
rmeddis@0 94 betweenRuns.bestThresholdTracks=betweenRuns.bestThresholdTracks(idx2);
rmeddis@0 95
rmeddis@0 96 if printReportGuide.structures
rmeddis@0 97 maxNoArrayValues=30;
rmeddis@32 98 showStructureSummary(stimulusParameters, 'stimulusParameters', ...
rmeddis@32 99 maxNoArrayValues)
rmeddis@0 100 showStructureSummary(experiment, 'experiment',maxNoArrayValues)
rmeddis@0 101 showStructureSummary(betweenRuns, 'betweenRuns',maxNoArrayValues)
rmeddis@0 102 showStructureSummary(withinRuns, 'withinRuns')
rmeddis@32 103
rmeddis@0 104 switch experiment.threshEstMethod
rmeddis@0 105 case {'2I2AFC++', '2I2AFC+++'}
rmeddis@32 106 showStructureSummary(LevittControl, 'LevittControl', ...
rmeddis@32 107 maxNoArrayValues)
rmeddis@0 108 end
rmeddis@32 109
rmeddis@0 110 switch experiment.ear
rmeddis@0 111 case {'statsModelLogistic','statsModelRareEvent'}
rmeddis@32 112 showStructureSummary(statsModel, 'statsModel', ...
rmeddis@32 113 maxNoArrayValues)
rmeddis@0 114 end
rmeddis@0 115 end
rmeddis@0 116
rmeddis@0 117 if printReportGuide.showTracks
rmeddis@0 118 % NB this procedure can only be used if all the tracks are present and
rmeddis@0 119 % of equal length
rmeddis@0 120 bigTable=[]; header=[];
rmeddis@32 121 disp(' ');
rmeddis@32 122 disp('Leveltracks starting from 1 response before the first reversal')
rmeddis@0 123 for i=1:length(betweenRuns.levelTracks)
rmeddis@0 124 if printReportGuide.HorizontalTracks
rmeddis@0 125 printTabTable(betweenRuns.levelTracks{i});
rmeddis@0 126 end
rmeddis@0 127 header=strvcat(header, 'level');
rmeddis@0 128 end
rmeddis@32 129
rmeddis@32 130 disp(' ');
rmeddis@32 131 disp('Response tracks from 1 response before the first reversal')
rmeddis@0 132 for i=1:length(betweenRuns.responseTracks)
rmeddis@0 133 if printReportGuide.HorizontalTracks
rmeddis@0 134 printTabTable(betweenRuns.responseTracks{i});
rmeddis@0 135 end
rmeddis@0 136 header=strvcat(header, 'resp');
rmeddis@0 137 end
rmeddis@32 138
rmeddis@0 139 disp(' '); disp('threshold tracks starting from the first reversal')
rmeddis@0 140 for i=1:length(betweenRuns.bestThresholdTracks)
rmeddis@0 141 if printReportGuide.HorizontalTracks
rmeddis@0 142 end
rmeddis@0 143 printTabTable(betweenRuns.bestThresholdTracks{i});
rmeddis@0 144 header=strvcat(header, 'mean');
rmeddis@0 145 end
rmeddis@32 146
rmeddis@0 147 end
rmeddis@0 148
rmeddis@0 149 switch experiment.ear
rmeddis@0 150 case {'MAPmodelMultiCh', 'MAPmodelSingleCh'}
rmeddis@0 151 % show all parameters but do not compute the model
rmeddis@0 152 nm=UTIL_paramsList(whos);
rmeddis@0 153 for i=1:length(nm)
rmeddis@0 154 try
rmeddis@0 155 eval(['UTIL_showStruct(' nm{i} ', ''' nm{i} ''')'])
rmeddis@0 156 catch
rmeddis@0 157 end
rmeddis@0 158 end
rmeddis@0 159 end
rmeddis@0 160
rmeddis@0 161 if experiment.saveData
rmeddis@0 162 fprintf('\n')
rmeddis@0 163 disp('To reprint this report with tracks use:')
rmeddis@0 164 disp([ 'printReport(''' saveFileName ''',1)'])
rmeddis@0 165 end
rmeddis@0 166
rmeddis@0 167 % print final summary (repeat of above)
rmeddis@0 168 fprintf('\n')
rmeddis@0 169 fprintf('\n')
rmeddis@0 170 disp('thresholds')
rmeddis@32 171 msg=printTabTable(sortTablesForPrinting...
rmeddis@32 172 (idx1,idx2,var1values,var2values, betweenRuns.thresholds), headers);
rmeddis@0 173 addToMsg(msg,0)
rmeddis@0 174 fprintf('\n')
rmeddis@0 175
rmeddis@32 176 if length(var1values)==1 && length(var2values)==1 ...
rmeddis@32 177 && experiment.maxTrials>49
rmeddis@32 178 [psy, levelsBinVector, binFrequencies]= ...
rmeddis@32 179 psychometricFunction(withinRuns.levelsPhaseTwo,...
rmeddis@32 180 withinRuns.responsesPhaseTwo, experiment.psyBinWidth);
rmeddis@0 181 disp('Psychometric function')
rmeddis@0 182 fprintf(' level \tfreq\tprob\n')
rmeddis@32 183 fprintf('%6.0f\t%6.2f\t%6.0f\n',[levelsBinVector; binFrequencies; psy])
rmeddis@38 184
rmeddis@38 185 switch experiment.threshEstMethod
rmeddis@38 186 % only one value required for level change
rmeddis@38 187 case {'MaxLikelihood', 'oneIntervalUpDown'};
rmeddis@38 188 % fprintf('\n')
rmeddis@38 189 % fprintf('k \t %6.2f\n',logistic.bestK)
rmeddis@38 190 % fprintf('g \t%7.5f\n',rareEvent.bestGain)
rmeddis@38 191 end
rmeddis@0 192 fprintf('\n')
rmeddis@32 193
rmeddis@0 194 end
rmeddis@9 195
rmeddis@9 196 fprintf('\nparadigm:\t%s\n ', experiment.paradigm)
rmeddis@32 197 if ~isempty(paramChanges)
rmeddis@33 198 fprintf('%s\n', char(paramChanges))
rmeddis@32 199 end
rmeddis@0 200
rmeddis@32 201 % ------------------------------------------------- sortTablesForPrinting
rmeddis@0 202 function table= sortTablesForPrinting(idx1,idx2, var1values,var2values, x)
rmeddis@0 203 % table converts a vector to a table
rmeddis@0 204 % after sorting according to idx1 and idx2
rmeddis@0 205 % the table is completed by adding row and column values
rmeddis@0 206 x=x(idx1);
rmeddis@0 207 x=x(idx2);
rmeddis@0 208 xMatrix=reshape(x,length(var1values),length(var2values));
rmeddis@0 209
rmeddis@0 210 table=[[-1000 var2values]; [var1values' xMatrix]];
rmeddis@0 211
rmeddis@32 212 % --------------------------------------------------- showStructureSummary
rmeddis@0 213 function showStructureSummary(structure, name, maxNoArrayValues)
rmeddis@0 214 % showStructureSummary prints out the values of a single structure
rmeddis@0 215 % The header is the structure name and each row is a field
rmeddis@0 216 % e.g. showStructureSummary(params,'params')
rmeddis@0 217 % This not the same as 'UTIL_showstruct'
rmeddis@0 218
rmeddis@0 219
rmeddis@0 220 if nargin<3
rmeddis@0 221 maxNoArrayValues=20;
rmeddis@0 222 end
rmeddis@0 223
rmeddis@0 224 fprintf('\n%s:', name)
rmeddis@0 225
rmeddis@0 226 fields=fieldnames(eval('structure'));
rmeddis@0 227 % for each field in the structure
rmeddis@0 228 for i=1:length(fields)
rmeddis@0 229 y=eval([ 'structure.' fields{i}]);
rmeddis@0 230 if isstr(y),
rmeddis@0 231 % strings
rmeddis@0 232 fprintf('\n%s=\t''%s''', fields{i},y)
rmeddis@0 233 elseif isnumeric(y)
rmeddis@0 234 % arrays
rmeddis@0 235 if length(y)>1
rmeddis@0 236 % vectors
rmeddis@0 237 [r c]=size(y);
rmeddis@0 238 if r>c, y=y'; end
rmeddis@32 239
rmeddis@0 240 [r c]=size(y);
rmeddis@0 241 if r>1
rmeddis@32 242 % fprintf('\n%s.%s=\t%g x %g matrix',name, fields{i}, r, c)
rmeddis@0 243 fprintf('\n%s=\t%g x %g matrix',fields{i}, r, c)
rmeddis@32 244
rmeddis@0 245 elseif c<maxNoArrayValues
rmeddis@32 246 % fprintf('\n%s=\t[%s]', fields{i},num2str(y))
rmeddis@0 247 fprintf('\n%s=', fields{i})
rmeddis@0 248 fprintf('\t%g',y)
rmeddis@32 249
rmeddis@0 250 else
rmeddis@32 251 fprintf('\n%s=\t %g... [%g element array]', ...
rmeddis@32 252 fields{i}, y(1),c)
rmeddis@0 253 end
rmeddis@0 254 else
rmeddis@0 255 % single valued arrays
rmeddis@32 256 % fprintf('\n%s.%s=\t%s;', name, fields{i},num2str(y))
rmeddis@0 257 fprintf('\n%s=\t%s', fields{i},num2str(y))
rmeddis@0 258 end
rmeddis@0 259 elseif iscell(y)
rmeddis@0 260 fprintf('\n%s=\t cell array', fields{i})
rmeddis@32 261
rmeddis@0 262 elseif isstruct(y)
rmeddis@0 263 fprintf('\n%s=\t structure', fields{i})
rmeddis@0 264 end,
rmeddis@32 265
rmeddis@0 266 end,
rmeddis@0 267 fprintf('\n')
rmeddis@0 268
rmeddis@0 269
rmeddis@0 270 % ------------------------------------------------------- printTabTable
rmeddis@0 271 function strings= printTabTable(M, headers)
rmeddis@0 272 % printTabTable prints a matrix as a table with tabs
rmeddis@0 273 %headers are optional
rmeddis@0 274 %headers=strvcat('firstname', 'secondname')
rmeddis@0 275 % printTabTable([1 2; 3 4],strvcat('a1','a2'));
rmeddis@0 276 stringCount=1; strings{stringCount}=[];
rmeddis@0 277
rmeddis@0 278 if nargin>1
rmeddis@0 279 [r c]=size(headers);
rmeddis@0 280 for no=1:r
rmeddis@0 281 % print all headers in a row
rmeddis@0 282 fprintf('%s\t',headers(no,:))
rmeddis@32 283 strings{stringCount}=...
rmeddis@32 284 sprintf('%s\t',headers(no,:)); stringCount=stringCount+1;
rmeddis@0 285 end
rmeddis@0 286 fprintf('\n')
rmeddis@0 287 end
rmeddis@0 288
rmeddis@0 289 [r c]=size(M);
rmeddis@0 290
rmeddis@0 291 for row=1:r
rmeddis@0 292 string=[];
rmeddis@0 293 for col=1:c
rmeddis@32 294 if row==1 && col==1 && M(1,1)==-1000
rmeddis@0 295 % Print nothing (tab follows below)
rmeddis@0 296 else
rmeddis@0 297 fprintf('%s',num2str(M(row,col)))
rmeddis@0 298 string=[string ' ' sprintf('%s',num2str(M(row,col)))];
rmeddis@0 299 end
rmeddis@0 300 if col<c
rmeddis@0 301 fprintf('\t')
rmeddis@32 302 %strings{stringCount}=sprintf('\t'); stringCount=stringCount+1;
rmeddis@0 303 end
rmeddis@0 304 end % col
rmeddis@0 305 strings{stringCount}=string; stringCount=stringCount+1;
rmeddis@0 306 fprintf('\n')
rmeddis@0 307 end % row
rmeddis@0 308
rmeddis@0 309 % ------------------------------------------------------- sortVariables
rmeddis@0 310 function [idx1, idx2, var1values, var2values]= ...
rmeddis@0 311 sortVariables(var1values, var2values, var1Sequence, var2Sequence)
rmeddis@0 312
rmeddis@0 313 [x idx1]= sort(var1Sequence);
rmeddis@0 314 var1Sequence= x;
rmeddis@0 315 var2Sequence= var2Sequence(idx1);
rmeddis@0 316 depVarName= 'th';
rmeddis@0 317
rmeddis@0 318 [x idx2]=sort(var2Sequence);
rmeddis@0 319 var2Sequence=x;
rmeddis@0 320 var1Sequence=var1Sequence(idx2);
rmeddis@0 321
rmeddis@0 322 var1values=sort(var1values);
rmeddis@0 323 var2values=sort(var2values);