comparison multithreshold 1.46/printReport.m @ 22:45f28c49461e master

removing duplicate changes
author Ray Meddis <rmeddis@essex.ac.uk>
date Mon, 13 Jun 2011 18:21:05 +0100
parents fafe69c43108
children 02aa9826efe0
comparison
equal deleted inserted replaced
21:c489ebada16e 22:45f28c49461e
24 printReportGuide.showTracks=experiment.printTracks; 24 printReportGuide.showTracks=experiment.printTracks;
25 printReportGuide.fileName=[]; 25 printReportGuide.fileName=[];
26 if experiment.saveData 26 if experiment.saveData
27 saveFileName=['savedData/' experiment.name '_' experiment.date '_' experiment.paradigm]; 27 saveFileName=['savedData/' experiment.name '_' experiment.date '_' experiment.paradigm];
28 else 28 else
29 % overwrite existing file just in case 29 % save this data (just in case)
30 saveFileName=['savedData/mostRecentResults']; 30 saveFileName=['savedData/mostRecentResults'];
31 end 31 end
32 experiment.minElapsed=etime(clock, betweenRuns.timeNow)/60; 32 experiment.minElapsed=etime(clock, betweenRuns.timeNow)/60;
33 save(saveFileName, 'experiment', 'stimulusParameters', 'betweenRuns', 'withinRuns', 'statsModel', 'expGUIhandles') 33 save(saveFileName, 'experiment', 'stimulusParameters',...
34 'betweenRuns', 'withinRuns', 'statsModel', 'expGUIhandles')
34 disp(['data saved as: ' saveFileName]); 35 disp(['data saved as: ' saveFileName]);
35 36
36 else 37 else
37 % reprint request 38 % reprint request (i.e print out old data)
38 printReportGuide.fileName=fileName; 39 printReportGuide.fileName=fileName;
39 load(printReportGuide.fileName) 40 load(printReportGuide.fileName)
40 saveFileName=printReportGuide.fileName; 41 saveFileName=printReportGuide.fileName;
41 if nargin>1 42 if nargin>1
42 printReportGuide.showTracks=printTracks; 43 printReportGuide.showTracks=printTracks;
44 printReportGuide.showTracks=experiment.printTracks; 45 printReportGuide.showTracks=experiment.printTracks;
45 end 46 end
46 end 47 end
47 48
48 fprintf('******** multithreshold ') 49 fprintf('******** multithreshold ')
50
49 x=pwd; disp(['version ' x(end-3:end)]) 51 x=pwd; disp(['version ' x(end-3:end)])
50 fprintf('\nName:\t%s ', experiment.name) 52 fprintf('\nName:\t%s ', experiment.name)
51 fprintf('\nparadigm:\t%s ', experiment.paradigm) 53 fprintf('\nparadigm:\t%s ', experiment.paradigm)
52 fprintf('\nEar:\t%s ', experiment.ear) 54 fprintf('\nEar:\t%s ', experiment.ear)
55
53 method=experiment.threshEstMethod; 56 method=experiment.threshEstMethod;
54 if stimulusParameters.includeCue && ~strcmp(method(1:6),'2I2AFC') 57 if stimulusParameters.includeCue && ~strcmp(method(1:6),'2I2AFC')
55 method=[method '/ withCue']; 58 method=[method '/ withCue'];
56 end 59 end
57 fprintf('\nmethod:\t%s ', method) 60 fprintf('\nmethod:\t%s ', method)
58 fprintf('\ndate:\t%s ', experiment.date) 61 fprintf('\ndate:\t%s ', experiment.date)
59
60 fprintf('\n\n') 62 fprintf('\n\n')
61 63
62 if isempty(betweenRuns.thresholds) 64 if isempty(betweenRuns.thresholds)
63 disp('no thresholds found') 65 disp('no thresholds found')
64 end 66 end
65 67
66 % prepare results as matrices ready to print 68 % prepare results as matrices ready to print
69 % first sort the actual sequence into a more readable sequence
67 [idx1, idx2, var1values, var2values]=... 70 [idx1, idx2, var1values, var2values]=...
68 sortVariables(betweenRuns.variableList1, betweenRuns.variableList2, betweenRuns.var1Sequence, betweenRuns.var2Sequence); 71 sortVariables(betweenRuns.variableList1, betweenRuns.variableList2,...
69 72 betweenRuns.var1Sequence, betweenRuns.var2Sequence);
70 % if strcmp(betweenRuns.variableName2, 'none') 73
71 % betweenRuns.variableName2=' ';
72 % end
73 header1=betweenRuns.variableName1; 74 header1=betweenRuns.variableName1;
74 header2=betweenRuns.variableName2; 75 header2=betweenRuns.variableName2;
75 header1 = strrep(header1, 'none', ' '); % none is not a useful header 76 header1 = strrep(header1, 'none', ' '); % none is not a useful header
76 header2 = strrep(header2, 'none', ' '); % none is not a useful header 77 header2 = strrep(header2, 'none', ' '); % none is not a useful header
77 headers=strvcat([header1 '/'], header2); 78 headers=strvcat([header1 '/'], header2);
78 79
79 disp('thresholds') 80 disp('thresholds')
80 msg=printTabTable(sortTablesForPrinting(idx1,idx2,var1values,var2values, betweenRuns.thresholds), headers); 81 msg=printTabTable(sortTablesForPrinting(idx1,idx2,...
82 var1values,var2values, betweenRuns.thresholds), headers);
81 addToMsg(msg,0) 83 addToMsg(msg,0)
82 fprintf('\n') 84 fprintf('\n')
83 85
84 % sort tracks into the same order 86 % sort tracks into the same order
85 betweenRuns.levelTracks=betweenRuns.levelTracks(idx1); 87 betweenRuns.levelTracks=betweenRuns.levelTracks(idx1);
86 betweenRuns.responseTracks=betweenRuns.responseTracks(idx1); 88 betweenRuns.responseTracks=betweenRuns.responseTracks(idx1);
87 betweenRuns.bestThresholdTracks=betweenRuns.bestThresholdTracks(idx1); 89 betweenRuns.bestThresholdTracks=betweenRuns.bestThresholdTracks(idx1);
88
89 betweenRuns.levelTracks=betweenRuns.levelTracks(idx2); 90 betweenRuns.levelTracks=betweenRuns.levelTracks(idx2);
90 betweenRuns.responseTracks=betweenRuns.responseTracks(idx2); 91 betweenRuns.responseTracks=betweenRuns.responseTracks(idx2);
91 betweenRuns.bestThresholdTracks=betweenRuns.bestThresholdTracks(idx2); 92 betweenRuns.bestThresholdTracks=betweenRuns.bestThresholdTracks(idx2);
92 93
93 if printReportGuide.structures 94 if printReportGuide.structures
128 header=strvcat(header, 'resp'); 129 header=strvcat(header, 'resp');
129 end 130 end
130 131
131 disp(' '); disp('threshold tracks starting from the first reversal') 132 disp(' '); disp('threshold tracks starting from the first reversal')
132 for i=1:length(betweenRuns.bestThresholdTracks) 133 for i=1:length(betweenRuns.bestThresholdTracks)
133 if printReportGuide.HorizontalTracks
134 end
135 printTabTable(betweenRuns.bestThresholdTracks{i});
136 header=strvcat(header, 'mean');
137 end
138
139 disp(' '); disp('threshold (mean) tracks starting from the first reversal')
140 for i=1:length(betweenRuns.bestThresholdTracks)
141 if printReportGuide.HorizontalTracks
142 end
143 printTabTable(betweenRuns.bestThresholdMeanTracks{i});
144 header=strvcat(header, 'mean');
145 end
146 disp(' '); disp('threshold tracks (median) starting from the first reversal')
147 for i=1:length(betweenRuns.bestThresholdMedianTracks)
148 if printReportGuide.HorizontalTracks 134 if printReportGuide.HorizontalTracks
149 end 135 end
150 printTabTable(betweenRuns.bestThresholdTracks{i}); 136 printTabTable(betweenRuns.bestThresholdTracks{i});
151 header=strvcat(header, 'mean'); 137 header=strvcat(header, 'mean');
152 end 138 end