Exercise

DESeq2 visualizations - MA and volcano plots

NOTE: It may take a bit longer to load this exercise.

To explore the results, visualizations can be helpful to see a global view of the data, as well as, characteristics of the significant genes. Usually, we expect to see significant genes identified across the range of mean values, which we can plot using the MA plot. If we only see significant genes with high mean values, then it could indicate an issue with our data. The volcano plot helps us to get an idea of the range of fold changes needed to identify significance in our data.

Let's explore our results using MA plots and volcano plots.

Instructions

100 XP

  • Create an MA plot using the plotMA() function and using the results object, smoc2_res as input.

  • Create a new column as a logical vector regarding whether padj values are less than 0.05 for the results using the mutate() function.

  • Create a volcano plot of the log2 foldchange values versus the -log10 adjusted p-value using ggplot() and coloring the points for the genes by whether or not they are significant.