function [binc, mu, lowci, highci, n, ytob]=binbino(x,y, varargin)
%  [binc, mu, se, n]=binbino(x,y,bin_e)
% Takes a vector x and a vector y and returns mean and binomial ci of
% values of y for bins of x.
%
% Input:
% x     1xn vector of x values to bin
% y     1xn vector of y values to average 
%
% Optional Input:
% 
% n_bins Optional # of bins. Default 10
% bin_e Optional bin edges.  Otherwise will return n_bins bins with equal # of samples from
%       min(x) to max(x)
%
% Output:
% binc      1xm bin centers
% mu        1xm The average value of y at that bin 
% se        1xm The standard error of y at that bin
% n         1xm The number of values of y in this bin


bin_e=[];
plot_it=false;
plot_fit='';
plot_fit_x0=[0 1 0 10];
n_bins=10;
alpha=0.01;
clr='k';
func=@nanmean;

utils.overridedefaults(who,varargin);

x=x(:);
y=y(:);
bad = isnan(x) | isnan(y);
x(bad)=[];
y(bad)=[];

ytob=nan+y;
if isempty(bin_e)
    pbins=linspace(0,100,n_bins+1);
    bin_e=prctile(x,pbins);
end
binc=(bin_e(2:end)+bin_e(1:end-1))/2;
mu=nan(size(binc));
se=nan(size(binc));

assert(all(y==0 | y==1 | isnan(y)),'This is for binomial data, use binned instead');

[n,~,yx]=histcounts(x,bin_e);

for nx=1:(numel(n))
    tmp=func(y(yx==nx));
    if isempty(tmp)
        tmp=nan;
    end
    mu(nx)=tmp;
    
    ytob(yx==nx)=nx;
    [~,ci]=binofit(nansum(y(yx==nx)),sum(yx==nx), alpha);         
    lowci(nx)=ci(1);
    highci(nx)=ci(2);
    
    % this is bad for low n.  
end


if plot_it
    
    draw.errorplot(gca,binc, mu, se,'Marker','o','Color',clr);
    if ~isempty(plot_fit)
        beta=nlinfit(x,y,plot_fit,plot_fit_x0);
        xx=linspace(binc(1),binc(end),1000);
        yy=plot_fit(beta,xx);
        hold on;
        plot(xx,yy,'-','Color',clr);
    end
    
    if binomial_data
        ylim([0 1]);
    end
    
    
end