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)".

## 文件列表
1. `likelihood.pdf`：以PDF格式保存的Matlab源代码文件。
2. `likelihood.htm`：以HTML格式保存的Matlab源代码文件。

## 数据集说明
本数据集包含一个Matlab函数，有关函数的执行方法，请参阅《Matlab 2001参考手册》。

该函数会基于函数内部指定的参数值，对给定数据对应的似然值进行计算。若函数运行正常，其返回值应为30.0050。

若需针对不同的参数与数据值计算似然值，请按以下步骤操作：
1) 使用`%`注释第6行与第7行（即`param=`和`data=`之前的两行）。
2) 将首行替换为`function [lik]=likelihood (param, data)`。
3) 将第4行替换为`randn ('state', sum(100*clock))`。
4) 将第5行替换为`rand ('state', sum(100*clock))`。
5) 在Matlab命令窗口中，将轨迹数据赋值给名为`data`的变量。
6) 同理，将参数值赋值给名为`param`的变量。
7) 通过命令`[lik]=likelihood (param, data)`执行该函数。