We might run Quest+ de novo in several distinct experimental sessions, and then want to use its data analysis routines to make estimates from all of the data combined. One way to do this is to set up a questPlus structure and then loop over all the data, updating each time through. Or, we might want to build up a representation directly. Providing some support routines to make this easier could be helpful. Note that in terms of trialData structure, Quest+ needs the trial index to be the row index. So you have to write (e.g.)
trialData(n,1).stim = stimOnThisNthTrial
not
trialData(n).stim = stimOnThisNthTrial
or preallocate the array, or something that has that effect. Maybe we could be more robust where Quest+ uses the information.
We might run Quest+ de novo in several distinct experimental sessions, and then want to use its data analysis routines to make estimates from all of the data combined. One way to do this is to set up a questPlus structure and then loop over all the data, updating each time through. Or, we might want to build up a representation directly. Providing some support routines to make this easier could be helpful. Note that in terms of trialData structure, Quest+ needs the trial index to be the row index. So you have to write (e.g.)
trialData(n,1).stim = stimOnThisNthTrial
not
trialData(n).stim = stimOnThisNthTrial
or preallocate the array, or something that has that effect. Maybe we could be more robust where Quest+ uses the information.