% script to run sigratio

run c:\mstuff\stat\stat185data\francesdata\dodata11 
% leaves linear in data11(:,1) , log10 in logdata11

% note tpvwin is just called here for plots, it is also called inside sigratio
% to identify and remove outliers. See sigratio parameter description below.

tpvwin(diff(data11),30,1,'data11');		

tpvwin(diff(logdata11),30,1,'logdata11');

pv1=sigratio(diff(data11),30,30,.01,1,'data11');

pv2=sigratio(diff(logdata11),30,30,.01,1,'logdata11');

%		function pv=sigratio(x,win,outwin,outalpha,plfg,plstr)
%
%  x = data vector
%  win = half of sigma comparison window
%  outwin = half width of outlier rejection window
%  outalpha = outlier false alarm prob
%  plfg plotflag, 1 to plot
%  plstr = data title string
%  pv is length(x) x 3  pv(:,1) is for sigsq2 > sigsq1, pv(:,2) is sigsq2 < sigsq1
%  and pv(:,3) is for sigsq2 ~= sigsq1, two sided