maskpsi2.m

clear all

useEyeLink = false;
useStaircase = true;
usePsi = false;
backgroundColour = [0.5 0.5 0.5];
subject = 'Ian';
fixx = 0;
fixy = 0;
pixelsPerCm = 44; %32=Lab CRT -- 44=27"monitor or Macbook Pro
nBlocks = 10; %number of repeated blocks?
windowed = [1024 768];

%-----spot stimulus
s = spotStimulus();
s.xPosition = -4;
t.contrast = 1;
s.size = 3;

sM = screenManager('verbose',false,'blend',true,'screen',0,'pixelsPerCm',pixelsPerCm,...
	'bitDepth','8bit','debug',true,'antiAlias',0,'nativeBeamPosition',0, ...
	'srcMode','GL_SRC_ALPHA','dstMode','GL_ONE_MINUS_SRC_ALPHA',...
	'windowed',windowed,'backgroundColour',[backgroundColour 0]); %use a temporary screenManager object
screenVals = open(sM); %open PTB screen
setup(s,sM); %setup our stimulus object

if useStaircase
	
	stopCriterion = 'trials';
	stopRule = 25;
	
	alphas = [0:0.05:1];
	priorB = PAL_pdfNormal(alphas, 0.25, 1);
	priorW = PAL_pdfNormal(alphas, 0.75, 1);
	
	PMB = PAL_AMRF_setupRF('xMin', 0,'xMax', 0.6,'PF',@PAL_Weibull,...
		'priorAlphaRange', alphas, 'prior', priorB,...
		'stopCriterion', stopCriterion, 'stopRule', stopRule);
	
	PMW = PAL_AMRF_setupRF('xMin', 0.4,'xMax', 1,'PF',@PAL_Weibull,...
		'priorAlphaRange', alphas, 'prior', priorW,...
		'stopCriterion', stopCriterion, 'stopRule', stopRule);
	
else
% 	%Set up psi
% 	NumTrials = 50;
% 
% 	grain = 201; %grain of posterior, high numbers make method more precise at the cost of RAM and time to compute.
% 	%Always check posterior after method completes [using e.g., :
% 	%image(PAL_Scale0to1(PM.pdf)*64)] to check whether appropriate
% 	%grain and parameter ranges were used.
% 
% 	PF = @PAL_Quick; %assumed psychometric function
% 
% 	%Stimulus values the method can select from
% 	stimRange = 0:0.01:1;
% 
% 	%Define parameter ranges to be included in posterior
% 	priorAlphaRange = linspace(0.1,0.9,grain);
% 	priorBetaRange =  linspace(0.01,3,grain); %Use log10 transformed values of beta (slope) parameter in PF
% 	priorGammaRange = 0;  %fixed value (using vector here would make it a free parameter)
% 	priorLambdaRange = .01; %ditto
% 	%tip: Free parameters sensibly and responsibly (as opposed to just because you can). See www.palamedestoolbox.org/understandingfitting.html.
% 
% 	%Initialize PM structure
% 	PM = PAL_AMPM_setupPM('priorAlphaRange',priorAlphaRange,...
% 		'priorBetaRange',priorBetaRange,...
% 		'priorGammaRange',priorGammaRange,...
% 		'priorLambdaRange',priorLambdaRange,...
% 		'numtrials',NumTrials,...
% 		'PF' , PF,...
% 		'stimRange',stimRange);
end

breakloop = false;
loop = 1;
figure;
hold on;
if rem(loop,2) == 1
	sM.backgroundColour = [0.3 0.3 0.3];
	bG = 0.3;
	s.colourOut = PMB.xCurrent;
else
	sM.backgroundColour = [0.7 0.7 0.7];
	bG = 0.7;
	s.colourOut = PMW.xCurrent;
end
update(s);

try
	%trial loop
	while breakloop ~= true
		drawBackground(sM);
		drawCross(sM,0.3,[1 0 0],fixx,fixy);
		vbl = Screen('Flip',sM.win); %flip the buffer
		WaitSecs(0.5)
		vbl = Screen('Flip',sM.win); %flip the buffer
		vbls = vbl;
		while GetSecs <= vbls+1
			drawBackground(sM);
			draw(s); %draw stimulus
			drawCross(sM,0.3,[1 0 0],fixx,fixy);
			Screen('DrawingFinished', sM.win); %tell PTB/GPU to draw
			nextvbl = vbl + screenVals.halfisi;
			vbl = Screen('Flip',sM.win, nextvbl); %flip the buffer
		end

		sM.backgroundColour = [0.5 0.5 0.5];
		drawBackground(sM);
		Screen('DrawText',sM.win,'BRIGHTER THAN (Left) or DARKER THAN (Right)',0,0);
		Screen('Flip',sM.win); %flip the buffer

		response = -1;
		ListenChar(2);
		[~, keyCode] = KbWait(-1);
		rchar = KbName(keyCode);
		if iscell(rchar);rchar=rchar{1};end
		switch lower(rchar)
			case {'leftarrow','left'}
				response = 1;
			case {'rightarrow','right'}
				response = 0;
			case {'escape','esc'}
				fprintf('\nQUIT!\n');
				breakloop = true;
		end
		ListenChar(0);

		%update PM based on response
		if response > -1
			if bG == 0.3
				if ~PMB.stop
					amplitude = PMB.xCurrent;
					PMB = PAL_AMRF_updateRF(PMB,amplitude,response);
					s.colourOut = PMW.xCurrent; update(s);
					if response == 1
						plot(length(PMB.x),PMB.x(end),'ro','MarkerFaceColor','r');
					else
						plot(length(PMB.x),PMB.x(end),'ko');
					end
					t = sprintf('C=%.2g | R=%i | T=%i | BG=%.2g\n',amplitude,response,loop,bG);
					fprintf(t);
					Screen('DrawText',sM.win,t,0,0);
				end
				Screen('Flip',sM.win); %flip the buffer
			else
				if ~PMW.stop
					amplitude = PMW.xCurrent;
					PMW = PAL_AMRF_updateRF(PMW,amplitude,response);
					s.colourOut = PMB.xCurrent; update(s); %mext trial black
					if response == 1
						plot(length(PMW.x),PMW.x(end),'bo-','MarkerFaceColor','b');
					else
						plot(length(PMW.x),PMW.x(end),'go-');
					end
					t = sprintf('C=%.2g | R=%i | T=%i | BG=%.2g\n',amplitude,response,loop,bG);
					fprintf(t);
					Screen('DrawText',sM.win,t,0,0);
				end
				Screen('Flip',sM.win); %flip the buffer
			end
		end
		
		if PMW.stop && PMB.stop
			fprintf('\nAdaptive Methods finished!\n');
			breakloop = true;
		end

		WaitSecs('YieldSecs',0.5);
		loop = loop + 1;
		if rem(loop,2) == 1
			sM.backgroundColour = [0.3 0.3 0.3];
			bG = 0.3;
		else
			sM.backgroundColour = [0.7 0.7 0.7];
			bG = 0.7;
		end
	end

	Screen('Flip',sM.win);
	Priority(0); ListenChar(0); ShowCursor;
	close(sM); %close screen
	clear sM s
	
catch
	Priority(0); ListenChar(0); ShowCursor;
	close(sM); %close screen
	clear sM s
	
end