Hoohm / CITE-seq-Count

A tool that allows to get UMI counts from a single cell protein assay
https://hoohm.github.io/CITE-seq-Count/
MIT License
79 stars 44 forks source link

method keeps looking for whitelist when not provided. #176

Closed naila53 closed 1 year ago

naila53 commented 1 year ago

Hi, I'm trying to run Cite-seq-count without a whitelist argument, however, it keeps being stuck at the looking for whitelist step after finishing the mapping for hours! image

Hoohm commented 1 year ago

Hello,

How many cells are you looking for? Might be that the parameter is set too high

Sent from my phone

On Thu, Feb 23, 2023, 15:47 naila53 @.***> wrote:

Hi, I'm trying to run Cite-seq-count without a whitelist argument, however, it keeps being stuck at the looking for whitelist step after finishing the mapping for hours! [image: image] https://user-images.githubusercontent.com/56326205/221026317-738a5e7a-c901-476d-89ee-535eaa962ef1.png

— Reply to this email directly, view it on GitHub https://github.com/Hoohm/CITE-seq-Count/issues/176, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAJVO2FR5KBLNMRCKM3H7BLWY7EEXANCNFSM6AAAAAAVGDNF74 . You are receiving this because you are subscribed to this thread.Message ID: @.***>

naila53 commented 1 year ago

Hi, I set the cells parameter to 1million, just thinking of retrieving empty droplets. That's too high right?

Hoohm commented 1 year ago

You should set it to the expected cells.

I think I add 20% to be sure to capture all possible cells.

Sent from my phone

On Fri, Feb 24, 2023, 16:13 naila53 @.***> wrote:

Hi, I set the cells parameter to 1million, just thinking of retrieving empty droplets. That's too high right?

— Reply to this email directly, view it on GitHub https://github.com/Hoohm/CITE-seq-Count/issues/176#issuecomment-1444588415, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAJVO2AX35JHD6L7JA2KEOTWZEXAVANCNFSM6AAAAAAVGDNF74 . You are receiving this because you commented.Message ID: @.***>

naila53 commented 1 year ago

hi, I need to retrieve the empty droplets for the sake of normalization, to model noise component. I have run your method about a year ago and you can see I posted a different issue here about totalseqB https://github.com/Hoohm/CITE-seq-Count/issues/158 but the point is I used high number of cells as you can see in my run info to retrieve the empty droplets. I tried again with inputting the raw unfiltered cellranger barcodes list as a whitelist for my data and kept the -cells argument at 1million cells , it still keeps getting stuck at the final step for hours! and it looks like it's trying to process more than 1million! I think something has changed in the new version, because it ran fine on all my previous totalseqA data.. image

Hoohm commented 1 year ago

Interesting.

I can't get into it atm as I'm on holiday. I suggest you try the older version then.

https://stackoverflow.com/questions/5226311/installing-specific-package-version-with-pip#5226504

The old version will give you the same counts as the new one.

Sent from my phone

On Sat, Feb 25, 2023, 00:53 naila53 @.***> wrote:

hi, I need to retrieve the empty droplets for the sake of normalization, to model noise component. I have run your method about a year ago and you can see I posted a different issue here about totalseqB #158 https://github.com/Hoohm/CITE-seq-Count/issues/158 but the point is I used high number of cells as you can see in my run info to retrieve the empty droplets. I tried again with inputting the raw unfiltered cellranger barcodes list as a whitelist for my data and kept the -cells argument at 1million cells , it still keeps getting stuck at the final step for hours! and it looks like it's trying to process more than 1million! I think something has changed in the new version, because it ran fine on all my previous totalseqA data.. [image: image] https://user-images.githubusercontent.com/56326205/221341004-475a0061-210e-41c5-af06-2c88a47cb257.png

— Reply to this email directly, view it on GitHub https://github.com/Hoohm/CITE-seq-Count/issues/176#issuecomment-1445004509, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAJVO2A62D6AZ4UTI422OA3WZGM47ANCNFSM6AAAAAAVGDNF74 . You are receiving this because you commented.Message ID: @.***>

naila53 commented 1 year ago

thank you so much for your prompt response! I tried reinstalling the older version( 1.4.3) and I seem to have resolved this. However, I tried to reproduce my previous analysis, I seem to have a new issue with both new and old version where it fails at the very last step! any ideas? :( image

Hoohm commented 1 year ago

This is related to the python version. I think it has to be under 3.6 to fix that

Sent from my phone

On Sat, Feb 25, 2023, 12:34 naila53 @.***> wrote:

thank you so much for your prompt response! I tried reinstalling the older version and I seem to have resolved this. However, I tried to reproduce my previous analysis, I seem to have a new issue with both new and old version where it fails at the very last step! [image: image] https://user-images.githubusercontent.com/56326205/221371131-8e1dd403-2a08-407c-bd84-e0b6efdefd52.png

— Reply to this email directly, view it on GitHub https://github.com/Hoohm/CITE-seq-Count/issues/176#issuecomment-1445165398, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAJVO2BODYGACSXLTPMALETWZI7AZANCNFSM6AAAAAAVGDNF74 . You are receiving this because you commented.Message ID: @.***>

naila53 commented 1 year ago

yes, it is! installed python 3.7 and worked! thank you so much!