ShellyCoder / cellcall

inferring cell-cell communication from scRNA-seq of ligand-receptor
73 stars 20 forks source link

A problem with TransCommuProfile() function while implementing vignettes codes #33

Open chunxubioinfor opened 2 years ago

chunxubioinfor commented 2 years ago

Hi! I'm new to CellCall and I like this tool after reading the paper. But when I try implementing the codes in the readme file, I encounter some problems.

#  Infer the cell-cell communication score
mt <- TransCommuProfile(object = mt,
                        pValueCor = 0.05,
                        CorValue = 0.1,
                        topTargetCor=1,
                        p.adjust = 0.05,
                        use.type="median",
                        probs = 0.9,
                        method="weighted",
                        IS_core = TRUE,
                        Org = 'Homo sapiens')

Unfortunately, this action can not be done and report the repeated errors like this:

no term enriched under specific pvalueCutoff...
no term enriched under specific pvalueCutoff...

And I try using warnings() to learn more about the errors:

Warning messages:
1: In cor(x, y, use = use, method = method) : the standard deviation is zero
2: In cor(x, y, use = use, method = method) : the standard deviation is zero
3: In cor(x, y, use = use, method = method) : the standard deviation is zero
4: In cor(x, y, use = use, method = method) : the standard deviation is zero
5: In cor(x, y, use = use, method = method) : the standard deviation is zero
6: In cor(x, y, use = use, method = method) : the standard deviation is zero
7: In cor(x, y, use = use, method = method) : the standard deviation is zero
8: In cor(x, y, use = use, method = method) : the standard deviation is zero

I don't how this happened, and I wonder if there's any solution for this. Thank you in advance.

ShellyCoder commented 1 year ago

Sorry for unconsiderable design of Cellcall. And the warning is ok.

  1. We computate the correlation between TF and Target gene. The warning, such as "Warning messages" was caused by the full zero expression data by one gene of TF or target.
  2. And the message "no term enriched under specific pvalueCutoff" means it did not enrich some TF in specific Cell.

This redundant process of information does come to unnecessary distress and I am very sorry. I will make further improvements.

Jiangmeling commented 1 year ago

I met the same problem. How did you solve it? Thank you