This function accepts a list of DE result sets and writes their outputs individually to file. Relies on init_cutoffs() significance thresholds. Filenames will be generated automatically based on the contrasts performed during the results() function. Output will be written to the /DE/data/ folder.
A list of DESeq result sets.
Also impose a filter on the DE set by a minimum absolute log2foldChange as determined by init_cutoffs(). Default=FALSE.
Output for each result set will be written to file This function does not return a value.
1 2 3 4 5 6 7 8 9 10 11 12 13 14
## 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) #Write differentally expressed gene data for each contrast to file. #Include a minimum fold change filter. This will output 3 files. write_all_de_results(res_list=myResList, lfc_filter=TRUE) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.