%% Testing functions

load example_data;

t=2; w=2;
l=1; k=1;

T1=zeros(1,11); %Rank
T2=zeros(1,11); %Partition
T3=zeros(1,11); %KDE
a=zeros(1,11);

for i=1:11

X1=coupling_data(i+33).X;
Y1=coupling_data(i+33).Y;
a(i)=coupling_data(i+33).a;

T1(i) = transferEntropyRank(X1,Y1,l,k,t,w,10);

[T2(i), nPar, dimPar]=transferEntropyPartition(X1,Y1,t,w);

T3(i) = transferEntropyKDE(X1,Y1,t,w,10,1);

end


figure;plot(a,T1,'.-'); hold on; plot(a,T2,'.-r'); plot(a,T3,'.-g'); 
xlabel('a'); ylabel('Transfer Entropy/Bit'); legend('TE Rank','TE Partition','TE KDE'); title('Transfer Entropy Between Time Series');