Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Results form R are different respect online website #81

Open
michela12345 opened this issue Jan 8, 2024 · 6 comments
Open

Results form R are different respect online website #81

michela12345 opened this issue Jan 8, 2024 · 6 comments

Comments

@michela12345
Copy link

When I perform xCell online throug the website i have results like this:

<style> </style>
  X1_11S_A X1_20_A X1_20_C X1_22_A X1_22_C X1_25_A
Adipocytes 0.0184 0.0464 0.0481 0.0357 0.0346 0.0444
Astrocytes 0.0052 0 0 0 0.0099 0
B-cells 0.0356 0 0 0.0056 0.0139 0.0192
Basophils 0.049 0.1879 0.201 0.1506 0.0749 0.1558
`

but when i do it on R i have this:

Adipocytes | 0.000000e+00 | 0.000000e+00 | 1.135611e-19 | 0.000000e+00

Astrocytes | 0.000000e+00 | 3.940393e-18 | 0.000000e+00

B-cells | 0.000000e+00 | 0.000000e+00

Basophils | 1.655104e-02 | 7.256572e-02 | 9.093160e-02 | 5.613090e-02

Do you think are the same? How I can convert the number on R?

And, when I run it online, i can download a p.value file but i do not understand what it is, can you explain me?

Thank you
Michela

@dviraran
Copy link
Owner

dviraran commented Jan 8, 2024

The online version uses rounding. You can try it your self -

scores = round(scores*1000)/1000

P-values are an attempt to assess whether the scores are significantly different than zero. If the p-values are high (>0.2) it means that the presence of this cell type in the mixture is low. I should note that providing confidence here is tough, so I won't make much of it.

Best,
Dvir

@michela12345
Copy link
Author

I round it as you say but the results are still different. My passages are:

devtools::install_github('dviraran/xCell',force = TRUE)
library(xCell)
xcell_r <- xCellAnalysis(d0)
xcell_r = round(xcell_r*1000)/1000

am i doing something wrong? do you think thant is better the version online or on R?

@michela12345
Copy link
Author

michela12345 commented Jan 8, 2024

I used the same matrix for both the analysis, and on the online website i condiser N=64 and I click on the RNA-seq bottom because my matrix is a bulk RNA-seq. My matrix looks like this
Screenshot 2024-01-08 163522
( the data are on excel so they have a , but when I upload the file on the website I use a txt file)

@dviraran
Copy link
Owner

dviraran commented Jan 9, 2024

Hi Michela,

I can't guess what is the issue, but you sent me the data and I can see from it that there are differences in the names of the samples and the order, so my intuition is that you are not using the same TPM matrix for both analyses. Beyond that I can't tell you what is going on. I tried myself one dataset and there where no differences.

I strongly recommend using the R version, with limited cell types which are relevant for your cell types.

https://link.springer.com/protocol/10.1007/978-1-0716-0327-7_19

Best,
Dvir

@michela12345
Copy link
Author

michela12345 commented Jan 9, 2024 via email

@michela12345
Copy link
Author

michela12345 commented Jan 9, 2024 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants