SCA-IRCM / SingleCellSignalR

R package for Bioconductor submission
29 stars 12 forks source link

cell_signaling() doesn't write into files #22

Open pcr08 opened 1 year ago

pcr08 commented 1 year ago

Good afternoon, I am still working on this package, and I have found a problem that I didn't see before. When applying cell_signaling() function, I get interactions between the different clusters, but then it doens't write them into the files (with write = TRUE parameter). The command line I used was: signal <- cell_signaling(data = data, c.names = c("microglial cells", "macrophages", "glioblastoma AC-like", "glioblastoma MES-like", "glioblastoma NPC-like", "glioblastoma OPC-like"), genes = rownames(data), cluster = clust, write=TRUE)

I asked the program to return the output from this function and also the warnings() that it gives, and this is what I got:

image

I get this for all the different cluster interactions (which is very weird if in the previous image it appears as existing interactions):

image

And the warnings() are all the same caused by the NA of the previous image:

image

Thank you, Best regards,

Paula

SCA-IRCM commented 1 year ago

Hi,

I respond on this thread but I think both of your problems are related because the visualize function cannot work with NAs and the cell_signaling function returns NA filled matrices. You need to check the maximum integer your system is capable of handling by typing .Machine$integer.max in your console and then verify that nrow(data) * ncol(data) isn't higher than this number. If it is the case I suggest you reduce the size of your data matrix by removing some lines. If you can't, you might think of processing your data 2 clusters by 2 clusters by making sub matrices.

Hope this resolves your issue.

SCA