Title: Refined parameter and uncertainty estimation when both variables are subject to error. Case study: estimation of Si consumption and regeneration rates in a marine environment
Authors: de Brauwere, A ×
De Ridder, F
Elskens, M
Schoukens, Joannes
Pintelon, R
Baeyens, W #
Issue Date: Apr-2005
Publisher: Elsevier science bv
Series Title: Journal of marine systems vol:55 issue:3-4 pages:205-221
Abstract: The problem of estimating parameters and their uncertainty from experimental measurements in marine ecosystems is a common task and often necessitates solving nonlinear equations. If the measurements are subject to individually varying errors (i.e., heteroscedastic data), the parameters are often estimated using a Weighted Least Squares (WLS) method. For estimating the parameter uncertainties, a linearized expression for the covariance matrix exists. Yet, both methods assume that the errors on the independent variable, also called "input", is negligible, which is often not true. For instance, in order to determine uptake and regeneration rates of silicic acid by phytoplankton, concentration and isotopic abundance measurements are performed at the beginning (input) and at the end (output) of an incubation experiment. Here, the so-called input and output are measurements of the same quantities, i.e., determined in exactly the same way, only differing by the time at which the measurements were performed. Clearly, there is no reason to assume that the input measurements are subject to less error than the output measurements. We propose a refinement of the two abovementioned estimation methods which enlarges their applicability to cases where input noise is not negligible. The refined methods are evaluated on the uptake and regeneration processes of silicic acid and compared to the original procedures using Monte-Carlo simulations. The results reveal a smaller bias for the refined WLS estimator compared with the original one. An additional advantage of using the refined WLS cost function is that its residual value can be interpreted as a sample from a chi(2) distribution. This property is especially useful because it enables an internal quality control of the results. In addition, the parameter uncertainty estimation is significantly improved. By neglecting the effect of the input noise, a (potentially) important origin of the parameter variation is simply ignored. Therefore, without the refinement, the parameter uncertainties are systematically underestimated. Using the refined method, this systematic error disappears and on the whole, the parameter standard deviations are accurately estimated. (c) 2004 Elsevier B.V. All rights reserved.
ISSN: 0924-7963
Publication status: published
KU Leuven publication type: IT
Appears in Collections:ESAT - STADIUS, Stadius Centre for Dynamical Systems, Signal Processing and Data Analytics
× corresponding author
# (joint) last author

Files in This Item:

There are no files associated with this item.

Request a copy


All items in Lirias are protected by copyright, with all rights reserved.

© Web of science