130 likes | 145 Views
Introduction to parameter optimisation. Sabine Beulke, CSL, York, UK FOCUS Work Group on Degradation Kinetics Estimating Persistence and Degradation Kinetics from Environmental Fate Studies in EU Registration Brussels, 26-27 January 2005. Curve fitting. Optimisation. Least squares method:
E N D
Introduction toparameter optimisation Sabine Beulke, CSL, York, UK FOCUS Work Group on Degradation Kinetics Estimating Persistence and Degradation Kinetics from Environmental Fate Studies in EU Registration Brussels, 26-27 January 2005
Optimisation Least squares method: Minimises the sum of squared residuals (RSS) Measured datapoint Calculated line Residual = deviation between calculated and measured data
Optimisation Initial guess (starting value) Calculate curve Calculate RSS Modify parameter
Automatic optimisation Stops when: • Convergence criteria are met Comparison between RSS for actual and previous runs. Convergence reached if difference is smaller than user-specified difference • Termination criteria are met For example, when maximum number of runs has been carried out (user-specified) Good fit not guaranteed!
Non-uniqueness Parameter correlation Parameters strongly related Effects on RSS of changes in one parameter can be compensated by changes in another parameter Inadequate model For example, selection of bi-phasic model not warranted if data follow SFO
Global versus local minimum RSS as a function of changes in 2 parameters The optimisation may find a local “valley” in the RSS surface, but not the absolute, global minimum. Different parameter combinations may be returned for different starting values. Good fit not guaranteed! From: http://www.ssg-surfer.com/
FOCUS recommendations • Always evaluate the visual fit • Avoid over-parameterisation • Aim at finding reasonable starting values • Always use different starting values • Constrain parameter ranges if appropriate • Plausibility checks for parameters and endpoints • Stepwise fitting where necessary • Be aware of differences between software packages
Goodness of fit - statistical criteria • 2 test where C = calculated value O = observed value = mean of all observed values err = measurement error percentage If calculated 2 > tabulated 2 then the model is not appropriate at the chosen level of significance Error percentage unknown Calculate error level at which 2 test is passed (e.g. with Excel spreadsheet provided by FOCUS)
Goodness of fit - statistical criteria • Confidence in parameter estimates Calculate e.g. from ModelMaker output A parameter is significantly different from zero if p (t) < alpha • Others (e.g. model efficiency, F-test)
FOCUS optimisation procedure Enter measured data Select kinetic model & parameters Initial guess (starting values) Change model, fix parameters? Eliminate outliers, weighting? Change starting values Evaluate: Visual fit Statistics Parameters Endpoints Optimise