2022 Agua Salud LMER Code RStudio

This is an RStudio script for a linear mixed effects model analysis using lmer. This linear mixed effects model was constructed to analyze the relative impacts of soil temperature, soil moisture, forest floor mass and soil organic C on soil CO2 fluxes (see dataset entitled 2022 Agua Salud CO2 Fluxes - Raw Data). First, a full mixed effects model was constructed using soil temperature (T), soil moisture (Moist), forest floor mass (Litter), and soil organic C (SOC) and their two-level interactions as fixed effects with Site as a random effect. Parameter estimates and sums of squares for each term were determined using summary(model) and anova(model) functions for the full model. To investigate the directions of the significant interaction terms identified in the full model, the deltaMethod function from the car package was used to generate transformations of the parameter estimates, resulting in point estimates and standard errors for CO2 flux. Data tables for the interaction plots were generated and plotted in a separate file (see 2022 Agua Salud Interaction Plots Code RStudio).