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);
|