fof-visually-guided/elutils/+stats/qbetween.m

function [y,xi]=qbetween(x, start, finish,refs)
% y=between(x, start, finish)
% returns the values of x that are between start and finish
% if start and finish are vectors y is a cell array
% works for sorted 1D vectors.
% using find you get o(n) , by assuming that the vector is sorted you get
% 2*o(log(n)).   which is WAY better.

if nargin==3
	refs=zeros(size(start));
end



if numel(start)>1
	for sx=1:numel(start)
		[y{sx} xi{sx}]=stats.qbetween(x,start(sx), finish(sx), refs(sx));
	end
else
	if isempty(x) || (x(end)<start) || x(1)>finish || isnan(start) || isnan(finish)
		y=[]; xi=[];
		return
	end
	i=stats.qfind(x, [start finish]);
	
	%qfind returns -1 if the target is less than min(x), since we are
	%getting 'between', we just take the first relevant x
	
	if i(1)>i(2)
		y=[]; xi=[];
		return
	elseif i(1)==0
		i(1)=1;
	end
	
	
	%this code deals with the fact that if there is no exact match, qfind
	%will return the index that is one lower than the target.  since we
	%want between, we just double check the end points.  every other point
	%will be valid.
	
	y=x(i(1):i(2));
	xi=i(1):i(2);
	if y(1)<start
		y=y(2:end);
		xi=xi(2:end);
    end
    
    if isempty(y)
        y = [];
        xi = [];
        return;
    end
	
	if y(end)>finish
		y=y(1:end-1);
		xi=xi(1:end-1);
	end
	
	y=y-refs;
	
	
end