Hello all
I would like to discuss something about the function "filtdir.m" for detection of coupling direction, downloaded from http://tocsy.pik-potsdam.de/.
I tried to reproduce the results for the problem of "three oscillators in a ring", from the publication ' Identification of coupling direction: Application to cardiorespiratory interaction, Phys. Rev. E, 65, 2002, 041909'.
But for every set of data i am only getting zeros for the couplings C1 & C2 and Nan for directionality index.
Should something be done in the mldivide step (i am finding zeros from that step...). could someone help me with this problem.
bindu