0001 function score=nt_repeatability(x,demean_flag)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012 if nargin<1; error('!'); end
0013 if ndims(x)~=3; error('data must be 3D'); end
0014
0015 if nargin<2 || isempty(demean_flag)
0016 demean_flag=1;
0017 end
0018
0019 if demean_flag==1
0020 x=nt_demean(x);
0021 elseif demean_flag==2
0022 x=nt_demean2(x);
0023 end
0024
0025 score=mean(mean(x,3).^2)./mean(nt_unfold(x.^2));
0026
0027 if nargout == 0 ;
0028 plot(score); xlabel('component'); ylabel('score'); title('repeatability');
0029 score=[];
0030 end
0031