Can I compare two sets of metadata information at the same time in function irGSEA.integrate()?

[Forsyl]

Many thanks to the authors for providing such a convenient tool for the integration and comparison of multiple gene set enrichment analysis methods!

For example, if there is still a set of metadata named “group” in the pbmc3k.final dataset which is used to mark whether a cell belongs to the control group or the experimental group, I would like to compare the difference between the control group and the experimental group of each cell type and display it by drawing a heatmap with the irGSEA.heatmap() function.

However, it seems that I can only draw a heatmap with only four columns of squares and symmetrical left-right for each cell type by filtering cells of different cell types separately, and then setting the “group.by” parameter of irGSEA.integrate() to “group”.

Can I use irGSEA.integrate() function to compare both “seurat_annotations” and “group” metadata, and finally show the difference between control and experimental groups for each cell type in one heatmap?