2022 Agua Salud LMER Code RStudio

Published: 20 May 2024| Version 2 | DOI: 10.17632/5sktkmg8x9.2
Contributor:
Claire Beckstoffer

Description

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).

Files

Categories

Forest Management, Tropical Soil, Soil Respiration, Tropical Forest Ecosystem

Funding

Smithsonian Tropical Research Institute

NSF Luquillo LTER DEB-1831952

Licence