Hi,
I am new to scRNAseq and I am very confused about few things when it1s time to integrate two datasets/conditions (control vs treated).
So I "characterized" each dataset (control or treated) and identified the various clusters (which cells are what) tutorial. But I need to integrate the two datasets to see what happens to the cells when I treat my tissue with a drug so to my understanding I should follow this integration but:
- why I have to split the dataset into 2 objects since in each dataset has one only condition? Should I not merge first all?
]2. once I merged, will I have to characterized again the dataset (which I assume contains only cells and genes shared across the two conditions) doing PCA, map, cluster identification? is this correct?
- After the integration (harmony) when I look the the pre vd post umps, there is not much difference (below an example). how do I find then only those few genes that change within each cluster?
I am sorry if these questions are stupid but I find a lot of script with very few explanations why should I use X parameter or X approach then I think I am doing random things without understand them.
Thank you
Camilla
Not answering your question but wanted to say that this eBook - https://bioconductor.org/books/release/OSCA/ is highly recommended by scRNA experts on biostars. Check it out.
Please accept the answer(s) that worked for you in your previous question.