Updated Palamedes parameters

TaskPipeline_v3.m

@@ -4,7 +4,7 @@
 %% enter subject ID %%
 %%%%%%%%%%%%%%%%%%%%%%
 
-subj = 'TEST_v1';
+subj = 'JQK_170707_v2';
 
 %% add required paths
 

TaskPipeline_v4.m

+% Pipeline including task practice/threshold estimation, MAT pre, checkerboard and MAT post assessment
+
+%%%%%%%%%%%%%%%%%%%%%%
+%% enter subject ID %%
+%%%%%%%%%%%%%%%%%%%%%%
+
+subj = 'JQK_170707_v2';
+
+%% add required paths
+
+pn.root = '/Volumes/LNDG-1/Julian/Projects/DYNSAT/C_Paradigm/';
+pn.CB = [pn.root, 'checker/'];
+pn.MAT = [pn.root, 'dotsx/'];
+
+addpath(genpath([pn.root, 'functions/']));
+addpath(pn.CB);
+addpath(pn.MAT);
+
+addpath(genpath('/Users/Shared/Psychtoolbox/')); % PTB 3.0.13 (160606)
+
+%% create results directory
+
+subjectsPath = [pn.root, '/data/', subj '_' datestr(now, 'yymmdd_HHMM') '/']; mkdir(subjectsPath);
+
+%% set up PTB
+
+Screen('Preference', 'SkipSyncTests', 1);
+oldLevel = Screen('Preference', 'Verbosity', 4); % output debugging info
+PsychDebugWindowConfiguration(0,0.3)
+
+%% practice MAT task
+
+cd(pn.root);
+
+% staircase parameters
+
+P.PM.numtrials          = 60;                                               % (numtrials/2) needs to be round
+P.PM.PF                 = @PAL_Gumbel;                                      % assumed psychometric function
+P.PM.grain              = 201;                                              % tradeoff between precision & computation time
+P.PM.stimRange          = .5:.02:.7;                                        % only allow values in range from 50 to 80%
+P.PM.priorAlphaRange    = linspace(.5,.7,P.PM.grain);
+P.PM.priorBetaRange     = -1:.05:3; 
+P.PM.priorGammaRange    = 0.5;                                              % fixed chance level
+P.PM.priorLambdaRange   = .05;                                              % fixed lapse rate
+
+dotInfo = [];
+dotInfo = createDotInfo_JQK_170704;
+dotInfo.dataDir = subjectsPath;
+
+DEBUG = 1; opacity = 1;
+cursubj = [subj, '_MAT_Practice'];
+thresholds = MAT_practice_170707(dotInfo, cursubj, DEBUG, opacity, P);
+
+disp(thresholds); disp('Please check whether they are reasonable.')
+
+%% run MAT Task (Pre-CB) 
+
+dotInfo = createDotInfo_JQK_170704(thresholds);
+dotInfo.dataDir = [pn.root, 'data/'];
+
+cursubj = [subj, '_MAT_Pre'];
+DEBUG = 1; opacity = 1;
+MAT_experiment_170704(dotInfo, cursubj, DEBUG, opacity)
+
+%% practice checkerboard task
+
+setup.subj          = subj;
+setup.design        = 1;
+setup.ET            = 0;
+setup.in_scanner    = 0;
+setup.dummyScans    = 0;
+setup.add_display   = 0; % if not in scanner: additional display connected?
+setup.test_trigger  = 0; % test scanner trigger (1) or not (0)
+setup.DEBUG         = 1; % debug mode?
+setup.transparency  = .99; % only relevant if DEBUG == 1
+setup.verbosity     = 4; % magnitude of PTB output
+setup.EPI           = 0;
+setup.BGmode        = 'black'; % black/grey mode
+setup.mode          = 'task';
+
+[checkerInfo_practice, setup] = PureCheckerParameters_Practice_170403(setup);
+setup.catchPractice = 0;
+setup.subjectsPath = subjectsPath;
+checkerBoard_Practice_Threshold_170705(checkerInfo_practice, setup);
+
+%% set up and run checkerboard
+
+setup               = [];
+setup.subj          = subj;
+setup.design        = 1;
+setup.ET            = 0;
+setup.in_scanner    = 0;
+setup.dummyScans    = 0;
+setup.add_display   = 0; % if not in scanner: additional display connected?
+setup.test_trigger  = 0; % test scanner trigger (1) or not (0)
+setup.DEBUG         = 1; % debug mode?
+setup.transparency  = .99; % only relevant if DEBUG == 1
+setup.verbosity     = 4; % magnitude of PTB output
+setup.EPI           = 0;
+setup.BGmode        = 'black'; % black/grey mode
+setup.mode          = 'task';
+
+[checkerInfo, setup] = PureCheckerParameters_170403(setup);
+setup.outputName = 'Task';
+setup.subjectsPath = subjectsPath;
+checkerContrastComp_170706(checkerInfo, setup)
+
+%% run MAT Task (post-CB)
+
+dotInfo         = [];
+dotInfo         = createDotInfo_JQK_170704(thresholds);
+dotInfo.dataDir = subjectsPath;
+
+cursubj = [subj, '_MAT_Post'];
+DEBUG = 1; opacity = 1;
+MAT_experiment_170704(dotInfo, cursubj, DEBUG, opacity)

dotsx/MAT_experiment_170704.m

@@ -13,6 +13,8 @@
 
 % adapted from dotsExperiment_JQK_170327
 % 170704    | cleanup, adapted
+% 170707    | removed saving for each trial; instead save in catch part
+%           | now saves DisplayInfo
 
 function MAT_experiment_170704(dotInfo, subj, DEBUG, opacity)
 
@@ -125,13 +127,12 @@ function MAT_experiment_170704(dotInfo, subj, DEBUG, opacity)
                     Screen('Flip', screenInfo.curWindow);
                 end
             end
-      save(sessionFile, 'ExperimentProtocol', 'dotInfo','ResultMat');
-      end % block
+      end % trial presentation
       
       sessEndTime = GetSecs;
       ExperimentProtocol = [ExperimentProtocol; {'SessionEnd'}, {sessEndTime}, {[]}, {[]}, {[]}, {indTrial}, {[]}];
 
-      save(sessionFile, 'ExperimentProtocol', 'dotInfo', 'ResultMat');
+      save(sessionFile, 'ExperimentProtocol', 'dotInfo', 'ResultMat', 'DisplayInfo');
 
     % inform the subject that the experiment is over
     DrawFormattedText(screenInfo.curWindow, 'task complete!', 'center', 'center');
@@ -159,6 +160,7 @@ function MAT_experiment_170704(dotInfo, subj, DEBUG, opacity)
 
   catch exception
     getReport(exception) % show stack trace
+    save(sessionFile);
   end
 
   closeExperiment; % close drawing canvas

dotsx/MAT_practice_170704.m

@@ -349,10 +349,14 @@ function [thresholds] = MAT_practice_170704(dotInfo, subj, DEBUG, opacity, P)
      %% 4th dimension
     
     % intialize psi parameters
-    PM = PAL_AMPM_setupPM('priorAlphaRange' , P.PM.priorAlphaRange, ...
-        'stimRange', P.PM.stimRange, 'numtrials', P.PM.numtrials, ...
-        'gamma', P.PM.gamma, 'priorLambdaRange', P.PM.priorLambdaRange, ...
-        'marginalize', 'slope', 'marginalize', 'lambda');
+    PM = PAL_AMPM_setupPM(...
+        'priorAlphaRange' , P.PM.priorAlphaRange, ...
+        'stimRange', P.PM.stimRange, ...
+        'numtrials', P.PM.numtrials, ...
+        'gamma', P.PM.gamma, ...
+        'priorLambdaRange', P.PM.priorLambdaRange, ...
+        'marginalize', 'slope', ...
+        'marginalize', 'lambda');
      
     % intiate parameters
     indTrial = 1;

dotsx/MAT_practice_170707.m

+%_______________________________________________________________________
+%
+% Run random dot motion experiment defined in createDotInfo
+%_______________________________________________________________________
+%
+% Input
+%
+% dotInfo | experiment configuration (struct)
+%    subj | identifier of subject (string)
+% condSet | conditions (cell array of strings)
+%_______________________________________________________________________
+%
+% 170707    | updated palamedes to take in more inputs
+
+function [thresholds] = MAT_practice_170707(dotInfo, subj, DEBUG, opacity, P)
+
+  %%
+  %% sanitize function parameters
+  %%
+  if nargin == 0
+    eval('help dotsExperiment')
+    return
+  end
+
+  if not(exist('dotInfo','var'))
+    error('Missing input: dotInfo')
+  elseif isempty(dotInfo) || not(isstruct(dotInfo))
+    error('Invalid input: dotInfo (requires non-empty struct)')
+  end
+
+  if not(exist('subj','var'))
+    error('Missing input: subj')
+  elseif not(isempty(subj)) && ischar(subj)
+    dotInfo.subj = subj;
+  else
+    error('Invalid input: subj (requires non-empty string)')
+  end
+
+  if not(isfield(dotInfo,'deployed'))
+    dotInfo.deployed = true;
+  end
+
+  %%
+  %% for saving behavioural data
+  %%
+
+  sessionFile = [dotInfo.dataDir, subj '_' datestr(now, 'yymmdd_HHMM'), '.mat'];
+    
+  %%
+  %% prepare presentation
+  %%
+  displayList = Screen('screens');
+  displayIdx = displayList(1);
+  KbName('UnifyKeyNames'); % use portable key naming scheme
+  oldVerbosityLevel = Screen('Preference', 'Verbosity', 2); % show errors and warnings
+  arrowKeys = [dotInfo.keyLeft dotInfo.keyRight];
+  baseDotInfo = dotInfo; % keep for saving later
+  
+    if DEBUG == 1
+        Screen('Preference', 'SkipSyncTests', 1);
+        PsychDebugWindowConfiguration(0, opacity);
+    else
+        clear Screen; %PsychDebugWindowConfiguration([], 1);
+    end
+
+  try
+
+    %%
+    %% prepare canvas
+    %%
+    screenInfo = openExperiment(50,50,displayIdx); clc; % open drawing canvas
+    if numel(Screen('screens')) == 1 || dotInfo.deployed
+      HideCursor(screenInfo.curWindow);
+    end
+    Screen('TextFont', screenInfo.curWindow, 'Helvetica');
+    Screen('TextSize', screenInfo.curWindow, 20);
+    Screen('TextStyle', screenInfo.curWindow, 0); % regular
+    Screen('TextColor', screenInfo.curWindow, 255); % white
+    
+    ResultMat = [];
+    
+    ExperimentProtocol = cell(0);
+    ExperimentProtocol = [ExperimentProtocol; {'SessionOnset'}, {GetSecs}, {[]}, {[]}, {[]}, {NaN}, {[]}];
+
+    %% Intro slide
+
+    DrawFormattedText(screenInfo.curWindow, '\nWillkommen bei der Multi-Attribut Aufgabe\n\nIm folgenden werden Sie Punktewolken sehen, die verschiedene Charakteristika aufweisen: \n\n Bewegungsrichtung \n Helligkeit \n Gre \n Farbe \n\n Im Folgenden sehen Sie ein (lngeres) Beispiel', 'center', 'center');
+    Screen('Flip', screenInfo.curWindow);
+
+    % wait for experimenter to resume session
+    while true
+    [exitSession, resumeSession] = checkKeys(dotInfo);
+    if resumeSession
+      break
+    elseif exitSession
+      closeExperiment;
+      return
+    end
+    end
+    
+    %% Show first 5 second example
+    
+    dotInfo.durCue = 0;     % duration of cue
+    dotInfo.durPres = 5;	% duration of presentation
+    dotInfo.durResp = 0;	% duration of question
+    dotInfo.durConf = 0;	% duration of confidence
+    
+    dotInfo.MAT.percAtt1H = .65;
+    dotInfo.MAT.percAtt1L = 1-dotInfo.MAT.percAtt1H;
+    dotInfo.MAT.percAtt2H = .65;
+    dotInfo.MAT.percAtt2L = 1-dotInfo.MAT.percAtt2H;
+    dotInfo.MAT.percAtt3H = .65;
+    dotInfo.MAT.percAtt3L = 1-dotInfo.MAT.percAtt3H;
+    dotInfo.MAT.percAtt4H = .65;
+    dotInfo.MAT.percAtt4L = 1-dotInfo.MAT.percAtt4H;
+    dotInfo.MAT.color = [255 255 255; 255 0 0];
+    dotInfo.MAT.coherence = 1;  % define dot movement (coherence) [useless now!]
+    dotInfo.MAT.direction = [180 0]; % left and right
+    dotInfo.MAT.size = [5 8]; % define dot size
+    dotInfo.MAT.luminance = [.4 .8]; % define dot luminance
+    dotInfo.MAT.attNames = {'color'; 'direction'; 'size'; 'luminance'};
+    
+    TrialInitiation = GetSecs;
+    ExperimentProtocol = [ExperimentProtocol; {'TrialInitiation'}, {TrialInitiation}, {[]}, {[]}, {[]}, {0}, {[]}];
+    targets = makeDotTargets(screenInfo, dotInfo); % initialize fixation targets
+    showTargets(screenInfo, targets, 1);
+    DisplayInfo.CombPositionByTrial = [];
+    DisplayInfo.CombSamples = [];
+    [~, dotInfo, ExperimentProtocol,ResultMat,DisplayInfo] = dotsX_JQK_MAT_170704(screenInfo, dotInfo, targets, 1, ExperimentProtocol,ResultMat,DisplayInfo);              
+
+    %% Instruction
+    
+    DrawFormattedText(screenInfo.curWindow, 'Nun sehen Sie die Attribute einzeln und in der spteren Anzeigedauer. \n Bitte geben Sie an, von welcher der beiden Optionen mehr vorhanden ist.', 'center', 'center');
+    Screen('Flip', screenInfo.curWindow);
+
+    pause(1);
+    
+    % wait for experimenter to resume session
+    while true
+    [exitSession, resumeSession] = checkKeys(dotInfo);
+    if resumeSession
+      break
+    elseif exitSession
+      closeExperiment;
+      return
+    end
+    end
+    
+    %% Instruction
+    
+    DrawFormattedText(screenInfo.curWindow, 'Wir beginnen mit der Farbe. \n\n Falls die weien Punkte dominieren, drcken Sie die linke Taste. \n Falls die roten Punkte dominieren, dcken Sie die rechte Taste. \n\n Im Anschluss an Ihre Antwort erhalten Sie eine Rckmeldung.', 'center', 'center');
+    Screen('Flip', screenInfo.curWindow);
+
+    pause(1);
+    
+    % wait for experimenter to resume session
+    while true
+    [exitSession, resumeSession] = checkKeys(dotInfo);
+    if resumeSession
+      break
+    elseif exitSession
+      closeExperiment;
+      return
+    end
+    end
+    
+    %% 1st dimension
+    
+    % intialize psi parameters
+    PM = PAL_AMPM_setupPM(...
+        'priorAlphaRange' , P.PM.priorAlphaRange, ...
+        'priorBetaRange', P.PM.priorBetaRange,...
+        'gamma', P.PM.priorGammaRange,...
+        'priorLambdaRange', P.PM.priorLambdaRange, ...
+        'stimRange', P.PM.stimRange, ...
+        'numtrials', P.PM.numtrials, ...
+        'PF' , P.PM.PF,...
+        'marginalize', 'slope', 'marginalize', 'lambda');
+     
+    % intiate parameters
+    indTrial = 1;
+    DisplayInfo.CombPositionByTrial = [];
+    DisplayInfo.CombSamples = [];
+    ResultMat = [];
+    % run the staircase algorithm until criterion is reached
+    dotInfo = [];
+    dotInfo = createDotInfo_TEST_170704(PM, 1);
+    dotInfo.deployed = true;
+    dotInfo.feedback = 1;
+    
+    ExperimentProtocol = [ExperimentProtocol; {'TrialInitiation'}, {GetSecs}, {[]}, {[]}, {[]}, {0}, {[]}];
+    while ~PM.stop
+        % Set the difference in amount of evidence for the target feature
+        % to the adaptive value. Note that for the remaining attributes,
+        % this is meaningless, as the different options are equal anyways.
+        disp(num2str(PM.xCurrent))
+        dotInfo.MAT.percAtt1H = PM.xCurrent;
+        dotInfo.MAT.percAtt1L = 1-dotInfo.MAT.percAtt1H;
+        targets = makeDotTargets(screenInfo, dotInfo); % initialize fixation targets
+        showTargets(screenInfo, targets, 1);
+        [~, dotInfo, ExperimentProtocol,ResultMat,DisplayInfo] = dotsX_JQK_MAT_170704(screenInfo, dotInfo, targets, indTrial, ExperimentProtocol,ResultMat,DisplayInfo);              
+        % update accuracy for the staircase code
+        trial.resp = ResultMat(indTrial,4);
+        % update PM structure
+        PM = PAL_AMPM_updatePM(PM, trial.resp);
+        indTrial = indTrial+1;
+    end
+    
+    Results{1,1} = dotInfo;
+    Results{2,1} = ExperimentProtocol;
+    Results{3,1} = ResultMat;
+    Results{4,1} = PM;
+    Results{5,1} = DisplayInfo;
+    
+    % estimate threshold
+    thresholds.color = PM.threshold(end);
+    
+    %% Instruction
+    
+    DrawFormattedText(screenInfo.curWindow, 'Als nchstes die Bewegungsrichtung. \n\n Falls die Bewegung nach links dominiert, drcken Sie die linke Taste. \n Falls die Bewegung nach rechts dominiert, dcken Sie die rechte Taste.', 'center', 'center', [255 255 255]);
+    Screen('Flip', screenInfo.curWindow);
+
+    pause(1);
+    
+    % wait for experimenter to resume session
+    while true
+    [exitSession, resumeSession] = checkKeys(dotInfo);
+    if resumeSession
+      break
+    elseif exitSession
+      closeExperiment;
+      return
+    end
+    end
+    
+    %% 2nd dimension
+    
+    % intialize psi parameters
+    PM = PAL_AMPM_setupPM(...
+        'priorAlphaRange' , P.PM.priorAlphaRange, ...
+        'priorBetaRange', P.PM.priorBetaRange,...
+        'gamma', P.PM.priorGammaRange,...
+        'priorLambdaRange', P.PM.priorLambdaRange, ...
+        'stimRange', P.PM.stimRange, ...
+        'numtrials', P.PM.numtrials, ...
+        'PF' , P.PM.PF,...
+        'marginalize', 'slope', 'marginalize', 'lambda');
+     
+    % intiate parameters
+    indTrial = 1;
+    DisplayInfo.CombPositionByTrial = [];
+    DisplayInfo.CombSamples = [];
+    ResultMat = [];
+    % run the staircase algorithm until criterion is reached
+    dotInfo = [];
+    dotInfo = createDotInfo_TEST_170704(PM, 2);
+    dotInfo.deployed = true;
+    dotInfo.feedback = 1;
+    
+    ExperimentProtocol = [ExperimentProtocol; {'TrialInitiation'}, {GetSecs}, {[]}, {[]}, {[]}, {0}, {[]}];
+    while ~PM.stop
+        disp(num2str(PM.xCurrent))
+        dotInfo.MAT.percAtt2H = PM.xCurrent;
+        dotInfo.MAT.percAtt2L = 1-dotInfo.MAT.percAtt2H;
+        targets = makeDotTargets(screenInfo, dotInfo); % initialize fixation targets
+        showTargets(screenInfo, targets, 1);
+        [~, dotInfo, ExperimentProtocol,ResultMat,DisplayInfo] = dotsX_JQK_MAT_170704(screenInfo, dotInfo, targets, indTrial, ExperimentProtocol,ResultMat,DisplayInfo);
+        % update accuracy for the staircase code
+        trial.resp = ResultMat(indTrial,4);
+        % update PM structure
+        PM = PAL_AMPM_updatePM(PM, trial.resp);
+        indTrial = indTrial+1;
+    end
+    
+    Results{1,2} = dotInfo;
+    Results{2,2} = ExperimentProtocol;
+    Results{3,2} = ResultMat;
+    Results{4,2} = PM;
+    Results{5,2} = DisplayInfo;
+    
+    % estimate threshold
+    thresholds.direction = PM.threshold(end);
+    
+    %% Instruction
+    
+    DrawFormattedText(screenInfo.curWindow, 'Als nchstes die Punktgre. \n\n Falls die kleineren Punkte dominieren, drcken Sie die linke Taste. \n Falls die greren Punkte dominieren, dcken Sie die rechte Taste.', 'center', 'center', [255 255 255]);
+    Screen('Flip', screenInfo.curWindow);
+
+    pause(1);
+    
+    % wait for experimenter to resume session
+    while true
+    [exitSession, resumeSession] = checkKeys(dotInfo);
+    if resumeSession
+      break
+    elseif exitSession
+      closeExperiment;
+      return
+    end
+    end
+    
+    %% 3rd dimension
+    
+    % intialize psi parameters
+    PM = PAL_AMPM_setupPM(...
+        'priorAlphaRange' , P.PM.priorAlphaRange, ...
+        'priorBetaRange', P.PM.priorBetaRange,...
+        'gamma', P.PM.priorGammaRange,...
+        'priorLambdaRange', P.PM.priorLambdaRange, ...
+        'stimRange', P.PM.stimRange, ...
+        'numtrials', P.PM.numtrials, ...
+        'PF' , P.PM.PF,...
+        'marginalize', 'slope', 'marginalize', 'lambda');
+     
+    % intiate parameters
+    indTrial = 1;
+    DisplayInfo.CombPositionByTrial = [];
+    DisplayInfo.CombSamples = [];
+    ResultMat = [];
+    % run the staircase algorithm until criterion is reached
+    dotInfo = [];
+    dotInfo = createDotInfo_TEST_170704(PM, 3);
+    dotInfo.deployed = true;
+    dotInfo.feedback = 1;
+    
+    ExperimentProtocol = [ExperimentProtocol; {'TrialInitiation'}, {GetSecs}, {[]}, {[]}, {[]}, {0}, {[]}];
+    while ~PM.stop
+        disp(num2str(PM.xCurrent))
+        dotInfo.MAT.percAtt3H = PM.xCurrent;
+        dotInfo.MAT.percAtt3L = 1-dotInfo.MAT.percAtt3H;
+        targets = makeDotTargets(screenInfo, dotInfo); % initialize fixation targets
+        showTargets(screenInfo, targets, 1);
+        [~, dotInfo, ExperimentProtocol,ResultMat,DisplayInfo] = dotsX_JQK_MAT_170704(screenInfo, dotInfo, targets, indTrial, ExperimentProtocol,ResultMat,DisplayInfo);
+        % update accuracy for the staircase code
+        trial.resp = ResultMat(indTrial,4);
+        % update PM structure
+        PM = PAL_AMPM_updatePM(PM, trial.resp);
+        indTrial = indTrial+1;
+    end
+    
+    Results{1,3} = dotInfo;
+    Results{2,3} = ExperimentProtocol;
+    Results{3,3} = ResultMat;
+    Results{4,3} = PM;
+    Results{5,3} = DisplayInfo;
+    
+    % estimate threshold
+    thresholds.size = PM.threshold(end);
+    
+    %% Instruction
+    
+    DrawFormattedText(screenInfo.curWindow, 'Als letztes die Punkthelligkeit. \n\n Falls die dunkleren Punkte dominieren, drcken Sie die linke Taste. \n Falls die helleren Punkte dominieren, dcken Sie die rechte Taste.', 'center', 'center', [255 255 255]);
+    Screen('Flip', screenInfo.curWindow);
+
+    pause(1);
+    
+    % wait for experimenter to resume session
+    while true
+    [exitSession, resumeSession] = checkKeys(dotInfo);
+    if resumeSession
+      break
+    elseif exitSession
+      closeExperiment;
+      return
+    end
+    end
+    
+     %% 4th dimension
+    
+    % intialize psi parameters
+    PM = PAL_AMPM_setupPM(...
+        'priorAlphaRange' , P.PM.priorAlphaRange, ...
+        'priorBetaRange', P.PM.priorBetaRange,...
+        'gamma', P.PM.priorGammaRange,...
+        'priorLambdaRange', P.PM.priorLambdaRange, ...
+        'stimRange', P.PM.stimRange, ...
+        'numtrials', P.PM.numtrials, ...
+        'PF' , P.PM.PF,...
+        'marginalize', 'slope', 'marginalize', 'lambda');
+     
+    % intiate parameters
+    indTrial = 1;
+    DisplayInfo.CombPositionByTrial = [];
+    DisplayInfo.CombSamples = [];
+    ResultMat = [];
+    % run the staircase algorithm until criterion is reached
+    dotInfo = [];
+    dotInfo = createDotInfo_TEST_170704(PM, 4);
+    dotInfo.deployed = true;
+    dotInfo.feedback = 1;
+    
+    ExperimentProtocol = [ExperimentProtocol; {'TrialInitiation'}, {GetSecs}, {[]}, {[]}, {[]}, {0}, {[]}];
+    while ~PM.stop
+        disp(num2str(PM.xCurrent))
+        dotInfo.MAT.percAtt4H = PM.xCurrent;
+        dotInfo.MAT.percAtt4L = 1-dotInfo.MAT.percAtt4H;
+        targets = makeDotTargets(screenInfo, dotInfo); % initialize fixation targets
+        showTargets(screenInfo, targets, 1);
+        [~, dotInfo, ExperimentProtocol,ResultMat,DisplayInfo] = dotsX_JQK_MAT_170704(screenInfo, dotInfo, targets, indTrial, ExperimentProtocol,ResultMat,DisplayInfo);        % update accuracy for the staircase code
+        trial.resp = ResultMat(indTrial,4);
+        % update PM structure
+        PM = PAL_AMPM_updatePM(PM, trial.resp);
+        indTrial = indTrial+1;
+    end
+    
+    Results{1,4} = dotInfo;
+    Results{2,4} = ExperimentProtocol;
+    Results{3,4} = ResultMat;
+    Results{4,4} = PM;
+    Results{5,4} = DisplayInfo;
+    
+    % estimate threshold
+    thresholds.luminance = PM.threshold(end);
+    
+    %% Instruction
+    
+    DrawFormattedText(screenInfo.curWindow, 'Nun kennen Sie alle Attribute. \n  Diese werden Ihnen im Folgenden gemeinsam prsentiert. \n Bitte geben Sie wie oben immer an, welche der beiden Optionen eines Attributs den Eindruck dominiert. \n Das relevante Attribut wird Ihnen dabei nach der Prsentation prsentiert.', 'center', 'center', [255 255 255]);
+    Screen('Flip', screenInfo.curWindow);
+
+    pause(1);
+    
+    % wait for experimenter to resume session
+    while true
+    [exitSession, resumeSession] = checkKeys(