Supplement 1. Matlab source code for the algorithm described in our study.

File List likelihood.pdf -Matlab source code saved in a pdf format. likelihood.htm -Matlab source code saved in a html format. Description This is a Matlab function. Please see Matlab Reference (Matlab 2001) for executing a function. The function evaluates likelihood associated with the data given in the function at parameter values also specified within the function. It should return the value 30.0050 if it is running correctly. To evaluate likelihood at different values of parameters and data, do the following: 1) Comment line 6 and 7 (before "param=" and "data=") with "%". 2) Replace the first line with "function [lik]=likelihood (param, data)". 3) Replace the fourth line with "randn ('state', sum(100*clock))". 4) Replace the fifth line with "rand ('state', sum(100*clock))". 5) In Matlab command window, assign size trajectory data to the variable named "data". 6) Similarly, assign parameter values to the variable named "param". 7) Execute the function using command "[lik]=likelihood (param, data)".