Many thanks for providing scvi-tools, it do helps a lot when I integrate multiple datasets.
I processed scATAC data with snapatac2, but snapatac2 output scATAC in bin counts (not peak counts). I tired many tools to integrate scRNA and scATAC with bin couns (5k) or gene activity matrix generated from bin counts (5k). I could not get proper integration results.
Could you suggest some about this? Do you think bin count could be used for peakVI or multiVI or could you suggest some preprocess pipeline used to merge multi scATAC data?
Hi, multiVI is intended to integrate e.g. multiome with RNA or ATAC data but not to integrate RNA with ATAC data. Most other models also require some cells with joint measurements. I’m not aware of other tools that I would recommend for this task.
Many thanks for your quick reply. I am integrating scRNA scATAC with multiomic data, so I first try integration with multiVI.
I went through the tutorial of multiVI, but in the tutorial, the test data (scRNA, scATAC and multiomic) are both subset of one multiomic dataset which means it will have the same var names. If I run multiVI with my data, should I align each sample with same var names or some overlapping between scATAC and multiomic is sufficient?