Hey @y.grabovska
Thank you for posting here. I think your analysis is correct.
- For integration SCVI-Tools need to have the raw counts.
- You can also think about the union of HVG’s per group, not only intersection. how much they differ? But, there are method that also compensate for the group size or batch like Scanpy’s
highly_variable_genes()withbatch_key(Python) or Seurat’sSelectIntegrationFeatures()(R). You can downsample the groups sizes to match each other and repeaet this process several times. - Its ok to manually remove HVG that are known to be in fact, noise.