Description
Fit a Collection of Curves to Single Cohort Decomposition Data.
Description
There is a long tradition of studying the flux of carbon from the biosphere to the atmosphere by following a particular cohort of litter (wood, leaves, roots, or other organic material) through time. The resulting data are mass remaining and time. A variety of functional forms may be used to fit the resulting data. Some work better empirically. Some are better connected to a process-based understanding. Some have a small number of free parameters; others have more. This package matches decomposition data to a family of these curves using likelihood--based fitting. This package is based on published research by Cornwell & Weedon (2013) <doi:10.1111/2041-210X.12138>.
README.md
litterfitter
R package for fitting and testing alternative models for single cohort litter decomposition data
Installation
#install.packages("remotes")
remotes::install_github("cornwell-lab-unsw/litterfitter")
library(litterfitter)
Getting started
At the moment there is one key function which is fit_litter
which can fit 6 different types of decomposition trajectories. Note that the fitted object is a litfit
object
fit <- fit_litter(time=c(0,1,2,3,4,5,6),
mass.remaining =c(1,0.9,1.01,0.4,0.6,0.2,0.01),
model="weibull",
iters=500)
class(fit)
You can visually compare the fits of different non-linear equations with the plot_multiple_fits
function:
plot_multiple_fits(time=c(0,1,2,3,4,5,6),
mass.remaining=c(1,0.9,1.01,0.4,0.6,0.2,0.01),
model=c("neg.exp","weibull"),
iters=500)
Calling plot
on a litfit
object will show you the data, the curve fit, and even the equation, with the estimated coefficients:
plot(fit)
The summary of a litfit
object will show you some of the summary statistics for the fit.
#> Summary of litFit object
#> Model type: weibull
#> Number of observations: 7
#> Parameter fits: 4.19
#> Parameter fits: 2.47
#> Time to 50% mass loss: 3.61
#> Implied steady state litter mass: 3.71 in units of yearly input
#> AIC: -3.8883
#> AICc: -0.8883
#> BIC: -3.9965
From the litfit
object you can then see the uncertainty in the parameter estimate by bootstrapping.