Mercurial > hg > map
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 |