Hi,

I have a oncoprint generated using oncoPrint() from ComplexHeatmap R package for variants for 20 genes from WGS data for 4 groups/phenotypes. Then, I ran differential gene expression for 2 groups and used plotCounts() from DESeq2 R package inside for loop to get 20 boxplots showing expression between 2 groups.

Now, I want to add those boxplots to each row in the oncoprint similar to what's shown here and I believe it is possible via following but I am struggling on how/which output from plotCounts for loop I should input to anno_boxplot():

top_annotation = HeatmapAnnotation(
                          boxplot = anno_boxplot(...output from plotcounts for loop....)

oncoPrint(mat, top_annotation = top_annotation)

Please advise.

Thanks.