NeuroKit: Offering help to investigate claim that cvxEDA is broken

I am the developer of the cvxEDA code used in NeuroKit.

I would like to offer my help to investigate what sort of issues cvxEDA is giving and to address them, with reference to the comment below in the code: https://github.com/neuropsychology/NeuroKit/blob/5cd089e9d82e6f7fb98d2028e69840aee08f8ab3/neurokit2/eda/eda_phasic.py#L17

There are two related issues in the cvxEDA repository (https://github.com/lciti/cvxEDA/issues/9 and https://github.com/lciti/cvxEDA/issues/10). The error reported there was raised by the cvxopt library with message “ValueError: Rank(A) < p or Rank([P; A; G]) < n”. I have downloaded data that was kindly provided by the user who reported the error and installed the same version of the libraries but I could not reproduce the error under any conditions.

Because cvxEDA depends on cvxopt, which in turn depends on cholmod, which in turn depends on a BLAS/LAPACK implementation, it’s very hard for me to reproduce the error and find a fix, unless I have the exact same data and libraries. I would be grateful if any person who receives this error could please export the exact list of packages with conda env export > env.yaml and attach it in reply to this issue, together with the data and platform. I am only interested in “reasonable” data and values of the parameters (e.g. the data cannot be just a handful of samples long, the sampling frequency cannot be 0.2 Hz or 500 MHz, because in these cases an error is reasonable).

There is the possibility that the error only happened with a specific platform and version of some of the libraries and that the underlying error in those libraries has then been fixed. In this case, I propose we monitor how cvxEDA is doing now and if it is not giving that error any more, we assume the problem has disappeared.

Thanks!