p147685
/
CRE_RB_and_CCRE
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

see previous comment

This commit is contained in:
R.J. Renken 2022-11-24 14:10:27 +01:00
parent cbcee6ac27
commit e6e9e593a4
1 changed files with 34 additions and 27 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

61
TestCREClass.m
View File

@ -2,31 +2,36 @@
% show some default behaviour % show some default behaviour
% and test the error handling. % and test the error handling.
%% set some constants %% set some constants
A=randn(100); % A=randn(1000,1);
% A=rand(1000,1);
A=cat(1,randn(1000,1),randn(200,1)+10);
%% create instance of CRE class %% create instance of CRE class
S1=CREClass; % % % S1=CREClass;
%% Calculate by setting labda equal to each of the (unique) values in the data. % % % %% Calculate by setting labda equal to each of the (unique) values in the data.
S1.UseHistProxyFlag=false; % % % S1.UseHistProxyFlag=false;
S1.EqualSizeBinFlag=false; % % % S1.EqualSizeBinFlag=false;
S1.Data=A;% set gaussian data % % % S1.Data=A;% set gaussian data
S1.Calc; % % % S1.Calc; %get RB based on the CRE
disp(S1) % % % S1.CalcRBShannon; %% get RB based on shannon entropy
%% use histogram aproximation with unequal bin sizes in the histogram % % % disp(S1)
S2=CREClass; % % % %% use histogram aproximation with unequal bin sizes in the histogram
S2.UseHistProxyFlag=true; % % % S2=CREClass;
S2.EqualSizeBinFlag=false; % % % S2.UseHistProxyFlag=true;
S2.Data=A;% set gaussian data % % % S2.EqualSizeBinFlag=false;
S2.nBin=numel(A); % set nbin to npoints; This is fine as we use the cummulative distribution and the formulas as defined in Zografos % % % S2.Data=A;% set gaussian data
S2.Calc % % % S2.nBin=numel(A); % set nbin to npoints; This is fine as we use the cummulative distribution and the formulas as defined in Zografos
disp(S2) % % % S2.Calc
%% use histogram aproximation with equal bin sizes in the histogram % % % S2.CalcRBShannon;
S3=CREClass; % % % disp(S2)
S3.UseHistProxyFlag=true; % % % %% use histogram aproximation with equal bin sizes in the histogram
S3.EqualSizeBinFlag=false; % % % S3=CREClass;
S3.Data=A;% set gaussian data % % % S3.UseHistProxyFlag=true;
S3.nBin=numel(A); % set nbin to npoints; This is fine as we use the cummulative distribution and the formulas as defined in Zografos % % % S3.EqualSizeBinFlag=false;
S3.Calc % % % S3.Data=A;% set gaussian data
disp(S3) % % % S3.nBin=numel(A); % set nbin to npoints; This is fine as we use the cummulative distribution and the formulas as defined in Zografos
% % % S3.Calc
% % % S3.CalcRBShannon;
% % % disp(S3)
%% Check recalculation of edges %% Check recalculation of edges
S4=CREClass; S4=CREClass;
S4.UseHistProxyFlag=true; S4.UseHistProxyFlag=true;
@ -34,13 +39,15 @@ S4.EqualSizeBinFlag=true;
S4.Data=A;% set gaussian data S4.Data=A;% set gaussian data
S4.nBin=numel(A); % set nbin to npoints; This is fine as we use the cummulative distribution and the formulas as defined in Zografos S4.nBin=numel(A); % set nbin to npoints; This is fine as we use the cummulative distribution and the formulas as defined in Zografos
S4.Calc S4.Calc
S4.CalcRBShannon;
disp(S4) disp(S4)
S4.nBin=numel(A)/10; S4.nBin=round(numel(A)/10);
S4.CalcEdges ; % force a recalculation of the edges S4.CalcEdges ; % force a recalculation of the edges
S4.Calc; S4.Calc;
S4.CalcRBShannon;
disp(S4) disp(S4)
%%
return;
%% Show effect of scaling or ofsett the data %% Show effect of scaling or ofsett the data
% Scale % Scale
S.Data=100*A; S.Data=100*A;