compops / qnmh-sysid2018

Constructing Metropolis-Hastings proposals using damped BFGS updates

Constructing Metropolis-Hastings proposals using damped BFGS updates

This code was downloaded from https://github.com/compops/qnmh-sysid2018 and contains the code and data used to produce the results in the paper:

J. Dahlin, A. Wills and B. Ninness, Constructing Metropolis-Hastings proposals using damped BFGS updates. Proceedings of the 18th IFAC Symposium on System Identification, Stockholm, Sweden, July 2018.

The paper is available as a preprint from https://arxiv.org/abs/1801.01243.

Python code (python/)

This code is used to set-up and run all the experiments in the paper. The code can possibly also be modified for other models. See the README.md file for more information.

Docker

A simple method to reproduce the results is to make use of the Docker container build from the code in this repository when the paper was published. Docker enables you to recreate the computational environment used to create the results in the paper. Hence, it automatically downloads the correct version of Python and all dependencies.

First, you need to download and installer Docker on your OS. Please see https://docs.docker.com/engine/installation/ for instructions on how to do this. Secondly, you can run the Docker container by running the command

docker run -v <<LOCALPATH>>:/app/results -e EXPERIMENT='<<EXPERIMENT_NUMBER>>' -e FULLRUN='<<FULL_RUN>>' --name qnmh-sysid2018-run compops/qnmh-sysid2018:final

where <<LOCALPATH>> is replaced with a local path to where the results are to be stored, e.g., /home/username/tmp on a Linux computer. <<EXPERIMENT_NUMBER>> is replaced with either 1, 2 or 3 corresponding to the three experiments in the paper. <<FULL_RUN>> can be either 0 or 1. Selecting 0 runs only one experiment (for practical reasons) and 1 runs 25 Monte Carlo simulations just like in the paper.

Note that running with the flag -e EXPERIMENT='2' -e FULLRUN='1', i.e. the second experiments (LGSS model with particle methods) and 25 Monte Carlo runs will probably take around a day to complete.

To reproduce the plots from the paper, move the contents of the results folder in <<LOCALPATH>> into a folder called results in a cloned version of the GitHub repository. Follow the instruction for the R code to create pdf versions of the plots. Note that the plot scripts are written for the case when -e FULLRUN='1' for the first and second example, which means that they might require some changes when running with -e FULLRUN='0'.

R code (r/)

This code is used to generate diagnostic plot as well as plots and table for the paper. See the README.md file for more information.

Binaries and results from simulations

The data generated from each run of the algorithms is about 100 megabytes. Therefore the generated data cannot be easily distributed via GitHub. Please contact the authors if you would like to receive a copy of the out from the simulations. Otherwise, you should be able to reproduce all runs yourself within a few hours by running the Docker container.

License

See the file LICENSE for more information.

###############################################################################
#    Constructing Metropolis-Hastings proposals using damped BFGS updates
#    Copyright (C) 2018  Johan Dahlin < uni (at) johandahlin [dot] com >
#
#    This program is free software: you can redistribute it and/or modify
#    it under the terms of the GNU General Public License as published by
#    the Free Software Foundation, either version 3 of the License, or
#    any later version.
#
#    This program is distributed in the hope that it will be useful,
#    but WITHOUT ANY WARRANTY; without even the implied warranty of
#    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
#    GNU General Public License for more details.
#
#    You should have received a copy of the GNU General Public License
#    along with this program.  If not, see <https://www.gnu.org/licenses/>.
###############################################################################

qnmh-sysid2018's Issues

issue with docker run

hey mate, cheers for the repository.
getting an error on the docker run.

gcc -pthread -shared build/temp.linux-x86_64-3.6/state/particle_methods/cython_lgss_helper.o -L/usr/local/lib -lpython3.6m -o /app/python/state/particle_methods/cython_lgss_helper.cpython-36m-x86_64-linux-gnu.so
building 'state.particle_methods.cython_sv_helper' extension
gcc -pthread -Wno-unused-result -Wsign-compare -DNDEBUG -g -fwrapv -O3 -Wall -Wstrict-prototypes -fPIC -I/usr/local/lib/python3.6/site-packages/numpy/core/include -I/usr/local/include/python3.6m -c state/particle_methods/cython_sv_helper.c -o build/temp.linux-x86_64-3.6/state/particle_methods/cython_sv_helper.o
gcc -pthread -shared build/temp.linux-x86_64-3.6/state/particle_methods/cython_sv_helper.o -L/usr/local/lib -lpython3.6m -o /app/python/state/particle_methods/cython_sv_helper.cpython-36m-x86_64-linux-gnu.so
File "run_script.py", line 10
elif int(sys.argv[1]) == 2:
^
IndentationError: unindent does not match any outer indentation level
File "run_script.py", line 10
elif int(sys.argv[1]) == 2:
^
IndentationError: unindent does not match any outer indentation level
File "run_script.py", line 10
elif int(sys.argv[1]) == 2:
^
IndentationError: unindent does not match any outer indentation level
tar: Removing leading `/' from member names
tar: /app/results/*: Cannot stat: No such file or directory
tar: Exiting with failure status due to previous errors