Tomasovych et al. supplement R code

Source code for model likelihoods and figures written in R language. Six density functions for fitting parameters of models with optimizing likelihoods are followed by age data of four depth assemblages. Fitting of age data with models occurs within the loop, followed by (1) figures showing the relation between the observed medians and IQRs and the values expected under different models, (2) likelihood surfaces for the Weibull and two-phase model, (3) simulations showing the effect of loss rate on the shape of AFD under one-phase and two-phase models, and (4) simulations showing the effects of sample size and fluctuating sinusoidal production on loss-rate parameters.