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

Some questions about the visualization by group #10

Closed
YiweiNiu opened this issue Aug 1, 2019 · 3 comments
Closed

Some questions about the visualization by group #10

YiweiNiu opened this issue Aug 1, 2019 · 3 comments
Labels
type:Bug Something isn't working

Comments

@YiweiNiu
Copy link

YiweiNiu commented Aug 1, 2019

Hi,

Thanks for this cool tool!

I have two questions about the visualization of .by group data.

Here is my data and metadata.

> names(tcr$data)
 [1] "0619_LN1" "0619_LN2" "0619_LN3" "0619_LN4" "0619_LN5" "0619_LN6" "0619_N1"  "0619_N2"  "0619_N3"  "0619_N4"  "0619_N5" 
[12] "0619_P1"

> tcr$meta
# A tibble: 12 x 4
   Sample   patient source tissue  
   <chr>    <chr>   <chr>  <chr>   
 1 0619_LN1 S0619   LN1    LN      
 2 0619_LN2 S0619   LN2    LN      
 3 0619_LN3 S0619   LN3    LN      
 4 0619_LN4 S0619   LN4    LN      
 5 0619_LN5 S0619   LN5    LN      
 6 0619_LN6 S0619   LN6    LN      
 7 0619_N1  S0619   N1     normal  
 8 0619_N2  S0619   N2     adjacent
 9 0619_N3  S0619   N3     tumor   
10 0619_N4  S0619   N4     tumor   
11 0619_N5  S0619   N5     tumor   
12 0619_P1  S0619   P1     PBMC  

First, I visualized the clonality by proportion.

tcr.imm_pr = repClonality(tcr$data, .method = "clonal.prop")
vis(tcr.imm_pr)

图片

vis(tcr.imm_pr, .by = 'tissue', .meta = tcr$meta, .test = F)

图片

From the individual samples. "0619_N3", "0619_N4" and "0619_N5" (from tumor) all have higher values than that of "0619_P1" (from PBMC). Why does PBMC have higher value than that of tumor after using .by = 'tissue'?

Second, visualized the clonal space homeostasis.

tcr.imm_hom = repClonality(tcr$data, .method = "homeo", .clone.types = c(Small = .0001, Medium = .001, Large = .01, Hyperexpanded = 1))
vis(tcr.imm_hom)

图片

vis(tcr.imm_hom, .by = c('tissue'), .meta = tcr$meta, .test = F)

图片

In the figure of individual samples, "0619_N2" does not have "Small" clones while all others have. But in the group view, "adjacent" (from "0619_N2") has over 20% "Small" clones. Also, nearly all samples have high proportion of "Medium" clones, but only tissue "LN" shows in the group view. So, how to understand this?

I am looking forward to hearing you.

Bests,
Yiwei Niu

@vadimnazarov
Copy link
Contributor

HI, thank you, will look into it and get back with updates or questions as quick as possible!

@vadimnazarov vadimnazarov added the type:Bug Something isn't working label Aug 6, 2019
@vadimnazarov
Copy link
Contributor

Hi, please update to the latest package version and check if everything is OK. The version should be 0.3.3.9006

@YiweiNiu
Copy link
Author

Thank you for your quick reply!

The problem is solved after using the newer version.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
type:Bug Something isn't working
Projects
None yet
Development

No branches or pull requests

2 participants