I'm analyzing an RNA-seq dataset using DESeq2, using a linear model that includes interaction terms. Some of these interactions are a bit hard for me to visualize (especially the 3-way interactions). I found this really great tutorial on linear models in R that includes some awesome visuals, for example:
But I don't know how to recapitulate such a graph for my data, and the code used to generate the above image was not provided. Ideally, I would like to tell R to make this type of plot for any gene of interest in my data using the model that DESeq2 generates. I'm a total newcomer to R, so any pointers in this direction would be greatly appreciated!