Teichlab / cellphonedb

MIT License
342 stars 105 forks source link

P-Value Thresholding for Dot-Plot #296

Closed aliciahillturner closed 3 years ago

aliciahillturner commented 3 years ago

I am trying to only show interactions with large circle diameters (i.e. significant interactions), therefore I want my dotplot to filter out any pvalues greater than 0.001 (for example). However, when I run the below code I still generate a dotplot that includes the non significant interactions (see dot plot below).

cellphonedb method statistical_analysis patient21meta.txt T21counts.txt --counts-data=gene_name --iterations=10 --threads=2 –-pvalue=0.001

image

Any help would be greatly appreciated as I am unsure what I am doing wrong with respect to p-value thresholding! :)

prete commented 3 years ago

Hi @aliciahillturner I may be getting your question wrong, but I think your figure looks as expected. If you want something else fileterd, you'll probably have to edit means.txt and pvalues.txt to get a different plot. Is there a row there you think it should not be included in the plot?

Pedramto89 commented 3 years ago

Hello @prete I want to do the same. I mean to insert a threshold for pvalue. For example greater and/or equal to 0.05. What's the command: Is sth like this:

ellphonedb method statistical_analysis meta_subtype2new.txt counts_subtype2new.txt --counts-data=hgnc_symbol --pvalue>0.05

Even I tried this, but did not work:Again I got the same error. ellphonedb method statistical_analysis meta_subtype2new.txt counts_subtype2new.txt --counts-data=hgnc_symbol --pvalue>=0.05

When I use the above command, I get the following error: Error: --pvalue option requires an argument