hidelab / diffcoexp Goto Github PK

Hi!
I'm running diffcoexp on a set of 27000 genes. Although all my genes have gene names and count values, I get several errors during the run, resulting in a DCL data.frame with mostly NA's. The DCGs are identified and the data frame with the results looks ok. When I run it in R locally, I get no warnings, but in a cluster I get the following:
<Warning messages:
1: In stats::cor(x, y, use, method) : the standard deviation is zero
2: In stats::cor(x, y, use, method) : the standard deviation is zero
3: In igraph::graph.data.frame(name.DCL) :
In d' NA' elements were replaced with string "NA"
4: In igraph::graph.data.frame(name.same) :
In d' NA' elements were replaced with string "NA">

I assume this is related with the DCL issue, but I could not find a way around it.

Thanks!

Hello,

I encountered this error today:

Error:

Finished running comparecor.
Finished running coexpr.
62060 gene pairs remain after half thresholding.
44049 DCLs identified.
Error in `[<-`(`*tmp*`, g.DCL.name, 2, value = degree.DCL) : 
  subscript out of bounds
In addition: Warning messages:
1: In graph.data.frame(name.DCL) :
  In `d' `NA' elements were replaced with string "NA"
2: In graph.data.frame(name.same) :
  In `d' `NA' elements were replaced with string "NA"