disp('Running MATLAB script file SiCur1fig1.m') ;
%
%    Makes Figure 1 for the SIgnificance of CURvature paper
%
%    Copied from Others\Old\rmtimes1.m
%
%    This analyzes Richard Smith's Raleigh Marathon Data



pstr = 'SiCur1fig1' ;



figh = figure(1) ;
clf ;


%  Read in data
%
infstr = '\Research\SiZer\Data\rmtimes.txt' ;
fid = fopen(infstr,'rt') ;
          %  'rt' is for "read only" and "text"
fstr = '%f %f' ;
[inmat,cnt] = fscanf(fid,fstr) ;
          %  can add a "size" parameter when don't want to read all.
          %  cnt tells how many reads were done (2 * 3 fields = 6)
fclose(fid) ;

n = cnt / 2 ;
inmat = reshape(inmat,2,n)' ;


%inmat
%gppause
          %  used this to test input

vdat = inmat(:,1) ;
titstr = 'Half Marathon Times' ;





%  do SiZer analysis
%
paramstruct = struct('iout',2,...
                     'imovie',0,...
                     'hhighlight',-1,...
                     'ibigdot',1,...
                     'famoltitle',[],...
                     'sizertitle',[],...
                     'curvsizertitle',[],...
                     'xlabelstr','time (min)',...
                     'ylabelstr',[],...
                     'labelfontsize',12,...
                     'ibdryadj',0,...
                     'iscreenwrite',1) ;

sizerSM(vdat,paramstruct) ;





  orient portrait ;
  set(gcf,'PaperSize',[5.0, 9.0]) ; 
  set(gcf,'PaperPosition',[0.25, 0.25, 4.5, 8.5]) ; 
print('-dpsc', [pstr '.ps']) ;


