Check if all genes are connected

[llrs]

From the https://github.com/llrs/integration-helper see the correct function (although it would require some tweaking to improve efficiency to the scales of gene sets)

[llrs]

Best path is to diagonalize and then compute the S (by adding all the diagonal matrix ( r-r^n)/(1-r). It might take time but there doesn't seem to be any other shortcut.

[llrs]

I don't know why but it returns NA because Reduce('+', red) returns all NaN