function plotPsycho(results,intensityName,p,functionName)
% plotPsycho(results,intensityName,p,functionName)

if ~exist('intensityName','var')
    intensityName= 'Intensity';
end


intensities = unique(results.intensity);

% Then we'll loop through these intensities calculating the proportion of
% times that 'response' is equal to 1:

nCorrect = zeros(1,length(intensities));
nTrials = zeros(1,length(intensities));

for i=1:length(intensities)
    id = results.intensity == intensities(i) & isreal(results.response);
    nTrials(i) = sum(id);
    nCorrect(i) = sum(results.response(id));
end

pCorrect = nCorrect./nTrials;

clf
hold on



 hold on
 
 sd = pCorrect.*(1-pCorrect)./sqrt(nTrials);  %pq/sqrt(n)
 errorbar(log(intensities),100*pCorrect,100*sd,'bo','MarkerFaceColor','b');
%  %loop through each intensity so each data point can have it's own size.
% for i=1:length(intensities);
%     sz = nTrials(i)+2;
%     plot(log(intensities(i)),100*pCorrect(i),'ko','MarkerFaceColor','b','MarkerSize',sz);
% end

if exist('p','var')
    %plot the parametric psychometric function 
    x = exp(linspace(log(min(results.intensity)),log(max(results.intensity)),101));
    evalStr = sprintf('y=%s(p,x);',functionName);
    eval(evalStr)
end

plot(log(x),100*y,'r-','LineWidth',2);

ylim  = get(gca,'YLim');
xlim = get(gca,'XLim');

pThresh = 100*(1/2)^(1/3);  
plot([xlim(1),log(p.t),log(p.t)],[pThresh,pThresh,ylim(1)],'k-');

set(gca,'XTick',log(intensities));
logx2raw
%set(gca,'YLim',[0,100]);
xlabel('Coherence');
ylabel('Percent Correct');

title(sprintf('Threshold: %5.2g',p.t));