For each organ, pairwise differential expression analysis between different ages with DESeq2. Fold change represents X over Y (so that X>Y if fold change is positive), if the internal naming scheme is X_vs_Y.
Examples of how to load and view results:
# load results
load("MACA_Bulk_deseq2_results.bin")
# list of tissues
ls(results_list_MACA_bulk)
# list of comparisons for Bone
ls(results_list_MACA_bulk$Bone)
# All results: comparison 3 vs 18 for Bone
results_list_MACA_bulk$Bone$`18_vs_3`$resall
# q<0.05 only: comparison 3 vs 18 for Bone
results_list_MACA_bulk$Bone$`18_vs_3`$ressig
