view m-toolbox/test/test_ao_spcorr.m @ 15:ce3fbb7ebe71
database-connection-manager
Remove broken functions from utils.jmysql
author
Daniele Nicolodi <nicolodi@science.unitn.it>
date
Mon, 05 Dec 2011 16:20:06 +0100 (2011-12-05)
parents
f0afece42f48
children
line source
+ − %
+ − % 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