Gene set variation analysis

GSVA calculates gene set enrichment scores (GSVA scores) for each sample using the z-statistic matrix.
Differential expression analysis is then conducted on these GSVA scores using limma, aiming to identify differentially regulated pathways across sample groups. Finally, the differential pathway activities across sample groups are visualized using a heatmap.