# NOT RUN {
#Prepare a result list for aggregation.
res.day1 <- results(dds, contrast=c("Condition_Time", "day1_disease", "day1_control"))
res.day2 <- results(dds, contrast=c("Condition_Time", "day2_disease", "day2_control"))
res.day3 <- results(dds, contrast=c("Condition_Time", "day3_disease", "day3_control"))
myResList <- list(res.day1, res.day2, res.day3)
#Plot differentially expressed genes for the aggregate result sets.
de_boxplot(res_list=myResList, filename="DE_condition_time_boxplot.pdf",
theme=2, returnData=FALSE)
# }
Run the code above in your browser using DataLab