tResults);

%% Granger causality test (GCT) and instantaneous GCT
%

gct_signif  = 0.001;  % Granger causality test significance level
igct_signif = 0.001;  % Instantaneous GCT significance level

flgPrintResults = 1;

[Tr_gct, pValue_gct]   =  gct_alg(u,A,pf, gct_signif,flgPrintResults);
[Tr_igct, pValue_igct] = igct_alg(u,A,pf,igct_signif,flgPrintResults);


%% Original PDC estimation
%
% PDC analysis results are saved in *c* structure.
% See asymp_pdc.m.

metric = 'euc';  % euc  = original PDC or DTF;
  