function times=statin_summt(X,ifshow) % times=statin_summt(X,ifshow) summarizes the event times in a STAToolkit input % % times{category}(firstlast,site)=first and last event times (NaN if none) % ifshow=1 to show category-by-category summary (default) % ifshow=2 to show a min-max summary % ifshow=3 for both % % See also: SWATCH2STATIN, TLS2SWATCH, STATIN_SUMM. % if (nargin<=1) ifshow=1; end % if (ifshow>0) disp(sprintf('number of sites: %3.0f',X.N)); disp(sprintf('number of categories: %3.0f',X.M)); disp(' categ trials [min max] times on each channel') dstring=repmat('[%8.3f %8.3f] ',1,X.N); end trials=zeros(1,X.M); times=cell(0); %make times_sites for each m times_sites=repmat([Inf -Inf]',[1 X.N X.M]); for m=1:X.M trials(m)=size(X.categories(m,1).trials,1); for n=1:X.N for p=1:trials(m) tvals=X.categories(m,1).trials(p,n).list; if length(tvals)>0 times_sites(1,n,m)=min(times_sites(1,n,m),min(tvals)); times_sites(2,n,m)=max(times_sites(2,n,m),max(tvals)); end end end times{m}=times_sites(:,:,m); if (mod(ifshow,2)==1) disp(sprintf(cat(2,'%7.0f %7.0f ',dstring),m,trials(m),times_sites(:,:,m))); end end if (ifshow>=2) disp(sprintf(cat(2,' all %7.0f ',dstring),sum(trials),... [min(times_sites(1,:,:),[],3);max(times_sites(2,:,:),[],3)])); end return