Confidence intervals of DEB parameters

The method for obtaining parameter confidence intervals is a two-step procedure. Steps can be performed in any order, but both steps are required for a final result. In the first step, the profile of the loss function for a parameter is obtained. In the second, which is the calibration step, the level of the loss function that corresponds to uncertainty is computed. Finally, the confidence interval is computed and the profile of the loss function and survival function are plotted.

In your working directory you need to have the 3 user-defined files of your species mydata_my_pet.m, predict_my_pet.m, and pars_init_my_pet.m (with the best estimates) + run_CI located in debtool_M/lib/pet_ci/.

Function run_ci performs the following steps:

  1. Calibration step for the interval estimation
    • calls get_lf_distribution
    • saves output to calibrate.mat
  2. Calculates the profile of the loss function for a parameter
    • calls get_profile
    • saves output to profile_nameOFparameter (e.g., profile_v)
  3. Plots the profile of the loss function for a parameter
    • Plots the distribution function of the loss function of Monte Carlo simulations
    • Gives the marginal confidence interval of the parameter
In run.ci the user must specify: species name (line 37); name of the parameter for which the profile will be calculated, the lower and upper limits of the range of simulations (lines 40-42); the confidence level (line 50).

The user can choose the step(s) to simulate (line58) by setting stepCI the value 0, 1, 2, or 3:

  • 0: runs the whole procedure;
  • 1: runs only the calibration step
  • 2: runs only the profile for a parameter step
  • 3: computes the CI and graphs the profile (calibrate.mat and profile_*.mat are required)