Open HMU-lyz opened 4 months ago
Hi @HMU-lyz,
I'm not sure I understand the issue 100%, but I assume you wonder why the receptor is set to FALSE given that the expression proportion for ERBB2 is high?
In this case, likely this complex is getting filtered because while EGFR is found and sufficiently expressed, while ERBB2 likely isn't present in the dataset and is hence set to false.
See, I use two steps in which, first I check if the genes from a complex are sufficiently expressed (I just consider whatever is the minimum expression), and I also check if all of the genes are actually found in the matrix.
I hope this answers the question, if you think that this is not the case, I'd be happy to explore the issue further.
Thanks for the kind words and for using liana :wink:
Thank you for your answer, my problem is solved
Your work is wonderful I encountered some problems when using it, which may not be technical reasons, but misunderstandings. When I opened the running results of liana_pipe, I found that there was only one corresponding gene for the complex subunit. Looking further at the source code, I found that in the recomplexify function, a flag was marked for each complex subunit. I don’t understand the principle of executing this step. In the subsequent steps, the FLASE row in the flag column will be removed, and then my problem arises. Here is the sample code I ran and it still exhibits this behavior:
This picture shows the problem when I run my own seurat: When I use the sample data, I still have similar problems, so I suspect it is a deviation in my understanding: