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

scirpy.read_10x_vdj(path, filtered=True) #121

Closed
taopeng1100 opened this issue Apr 22, 2020 · 6 comments
Closed

scirpy.read_10x_vdj(path, filtered=True) #121

taopeng1100 opened this issue Apr 22, 2020 · 6 comments

Comments

@taopeng1100
Copy link

Hi I like to read multiple "filtered_contig_annotations.csv" files from 10x; how can I do this in scirpy.read_10x_vdj?

Thx,

Tao

@grst
Copy link
Collaborator

grst commented Apr 22, 2020

Hi Tao,

you can read them in a loop and concatenate them.
It works the same way as for adata files without TCR information, see the solution here:
scverse/scanpy#267

Hope that helps and let me know if you have further problems!

Best,
Gregor

@taopeng1100
Copy link
Author

taopeng1100 commented Apr 22, 2020 via email

@grst grst mentioned this issue Apr 29, 2020
11 tasks
@JLTrincado
Copy link

JLTrincado commented May 4, 2020

Dear authors,

I'm using the solution proposed here for loading multiple 10x vdj files. For some reason, after concatenating some cells loses information. Here I load 2 files and concatenate them. Prior to concatenation, the shown cell has TRA_1_cdr3 region. After concatenating, the region has disappeared. This only happens with some cells, not all. What am I missing?

image

Thanks for your help.
Juan L.

@grst
Copy link
Collaborator

grst commented May 4, 2020

Hi Juan,

could it be just the order of columns?
The concatenation seems to sort columns alphabetically, s.t. the TRB_1_cdr3 column is hidden in the middle (...) or the pandas data frame.

Best,
Gregor

@JLTrincado
Copy link

Hi again,

thanks for the quick reply. It wasn't because of the order in the columns, but thanks to that I realize it was my bad when calculating the TCR neighbors. I was using receptor_arms="all" instead of "any". That's why this cell didn't show up as clone.

Thanks again and congratulations for your work.
Best.

@grst
Copy link
Collaborator

grst commented May 19, 2020

@taopeng1100, I will close this now. Let me know if you have further issues.

@grst grst closed this as completed May 19, 2020
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

3 participants