+
−% 
+
−% A test script for spcorr
+
−% 
+
−% L Ferraioli 30-06-2011
+
−%
+
−% $Id: test_ao_spcorr.m,v 1.1 2011/07/06 15:06:42 luigi Exp $
+
−%
+
−%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
−%% useful vars
+
−nsecs = 1000;
+
−fs = 1;
+
−
+
−%% test behaviour
+
−
+
−alp = 0.05;
+
−plcorr = plist('ALPHA',alp);
+
−
+
−a1 = ao.randn(nsecs,fs);
+
−a2 = ao.randn(nsecs,fs);
+
−a3 = ao.randn(nsecs,fs);
+
−a4 = ao.randn(nsecs,fs);
+
−
+
−b = spcorr(a1,a2,a3,a4,plcorr);
+
−
+
−%% test statistic hypotheis test fairness
+
−
+
−Nmc = 1000;
+
−
+
−alp = 0.05;
+
−plcorr = plist('ALPHA',alp);
+
−
+
−TestRes = zeros(Nmc,1);
+
−% run Monte Carlo
+
−for ii=1:Nmc
+
−  a1 = ao.randn(nsecs,fs);
+
−  a2 = ao.randn(nsecs,fs);
+
−  
+
−  b = spcorr(a1,a2,plcorr);
+
−  % get procinfo
+
−  bp = b.procinfo;
+
−  % get test result
+
−  tr = find(bp,'TestRes');
+
−  TestRes(ii) = double(tr);
+
−  
+
−end
+
−
+
−% the sum of the elements of TestRes should be a number near Nmc*alp
+
−sum(TestRes)
+
−Nmc*alp