Description Usage Arguments Value See Also Examples
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.
1 | write_all_de_results(res_list, lfc_filter = FALSE)
|
res_list |
A list of DESeq result sets. |
lfc_filter |
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.
init_cutoffs
, create_dir_struct
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.