Code covered by the BSD License  

Highlights from
ProcessNetwork Version 1.0 Software

image thumbnail

ProcessNetwork Version 1.0 Software

by

 

Functions for the delineation of Dynamical Process Networks using Information Theory

[InormByDist,TnormByDist,SigThreshInormByDist,SigThreshTnormByDist,Ic,Tc,TvsIzero,SigThreshTvsIzero,RelEnt,RelT]=NormTheStats(nBinMat,I,T,SigThreshI,SigThreshT,meanShuffI,sigmaShuffI,meanShuffT,sigmaShuffT,HXt,HYw,HYf)
function [InormByDist,TnormByDist,SigThreshInormByDist,SigThreshTnormByDist,Ic,Tc,TvsIzero,SigThreshTvsIzero,RelEnt,RelT]=NormTheStats(nBinMat,I,T,SigThreshI,SigThreshT,meanShuffI,sigmaShuffI,meanShuffT,sigmaShuffT,HXt,HYw,HYf)

[nSignals,~,nLags]=size(I);

InormByDist=NaN(nSignals,nSignals,nLags);
TnormByDist=NaN(nSignals,nSignals,nLags);
SigThreshInormByDist=NaN(nSignals,nSignals);
SigThreshTnormByDist=NaN(nSignals,nSignals);
Ic=NaN(nSignals,nSignals,nLags);
Tc=NaN(nSignals,nSignals,nLags);
TvsIzero=NaN(nSignals,nSignals,nLags);
SigThreshTvsIzero=NaN(nSignals,nSignals);
RelEnt=NaN(nSignals,nSignals,nLags);
RelT=NaN(nSignals,nSignals,nLags);

for i=1:nSignals
    for j=1:nSignals
        for t=1:nLags
            n=min(nBinMat(i),nBinMat(j));
            InormByDist(i,j,t)=I(i,j,t)/log2(n);
            TnormByDist(i,j,t)=T(i,j,t)/log2(n);
            Ic(i,j,t)=0.5*(1+erf((I(i,j,t)-meanShuffI(i,j))/(sqrt(2)*sigmaShuffI(i,j))));
            Tc(i,j,t)=0.5*(1+erf((T(i,j,t)-meanShuffT(i,j))/(sqrt(2)*sigmaShuffT(i,j))));
            SigThreshInormByDist(i,j)=SigThreshI(i,j)/log2(n);
            SigThreshTnormByDist(i,j)=SigThreshT(i,j)/log2(n);
            TvsIzero(i,j,t)=T(i,j,t)/I(i,j,1);
            SigThreshTvsIzero(i,j)=SigThreshT(i,j)/I(i,j,1);
            RelEnt(i,j,t)=I(i,j,t)/HYf(i,j,t);
            RelT(i,j,t)=T(i,j,t)/HYf(i,j,t);
        end
    end
end

%-----------------------------------------------------------

Contact us