function report=serialcorr(randomseq,maxshift)
%  report=serialcorr(randomseq,maxshift)
%
%  serial correlation function in a random sequence at a range of separations
%
%    author: Qianqian Fang (q.fang <at> neu.edu)
%
%    input:
%        randomseq: a random sequence (1D array)
%        maxshift:  the maximum separation to test with
%
%    output:
%        report: the corr. coeff for the sequence between randomseq 
%            and randomseq(i:end) where i<=maxshift
%
%    this file is part of Monte Carlo eXtreme (MCX)
%    License: GPLv3, see http://mcx.sf.net for details
%    see Boas2002, Heskell1996

report=zeros(maxshift,2);
for i=1:maxshift
    [r,p]=corrcoef(randomseq(1:end-i),randomseq(i+1:end));
    report(i,:)=[r(1,2),p(1,2)];
    fprintf(1,'test shift %d (r=%20.16f p=%20.16f)\n',i,r(1,2),p(1,2));
end