R code for Statistical Analysis and plots related to the article "Is taurine an aging biomarker?"

This collection contains R code for Statistical Analysis and plots related to our article "Is taurine an aging biomarker?" [ M. E. Fernandez et al., Science 388, eadl2116 (2025). DOI: 10.1126/science.adl2116 ]“Analysis” folderThis folder contains R code for all statistical analysis associated to our article "Is taurine an aging biomarker?"“Main Manuscript Plots” folderThis folder contains the plots provided in the figures of the main manuscript, together with equivalent plots (same data and results) that present different y-axis limits, allowing to better assess variability within each cohort. This is indicated in the name of the plot. E.g., the string “y0-75” in the name “1Ba_Balearic-F(5x11,2.5pt,y0-75)” indicates that the limits of the y-axis were set between 0-75, which corresponds to the range of taurine values in this cohort (in contrast to the main manuscript figure, which presents the y-axis limits between 0-336 for ease of comparison across cohorts). Similarly, the characters in the string “1Ba_” indicate that this plot is equivalent to the plot in Figure 1, panel B, left side, respectively. The full name of the cohort is indicated in the name of the plot. Plots for females and males are indicated with “F” and “M”, respectively.“Supplementary Material Plots” folderThis folder contains the plots provided in the figures of the supplementary material, together with equivalent plots (same data, same results) that present different y-axis limits that allow to better assess variability within each cohort. This is indicated in the name of the plot. E.g., the string “y0-75” in the name “S1Ba_Balearic-F(5x11,2.5pt,y0-75)” indicates that the limits of the y-axis were set between 0-75, , which corresponds to the range of taurine values in this cohort (in contrast to the main manuscript figure, which presents the y-axis limits between 0-336 for ease of comparison across cohorts). Similarly, the characters in the string “1Ba_” indicate that this plot is equivalent to the plot in Figure 1, panel B, left side, respectively. The full name of the cohort is indicated in the name of the plot. Plots for females and males are indicated with “F” and “M”, respectively.“Complementary Plots” folderThis folder contains plots to visualize the results presented in supplementary tables S40-S45 of our article (for details see scripts "Analysis1_BLSA-females-adjDiet.R", "Analysis1_BLSA-males-adjDiet.R", “Analysis1_BLSA-females-excluding-outliers.R”, “Analysis1_BLSA-males-excluding-outliers.R”, “Analysis1_BLSA-females-onlyActive.R”, “Analysis1_BLSA-males-onlyActive.R”, in the folder “Analysis”). None of these plots are provided in the figures of the main manuscript or the supplementary material. The names of the plots follow the same rules described above.