Michael Creel

michael.creel @ uab.es

Tel. +34 93 581 1696

Despatx (office) B3-1104

Departament d'Economia i d'Història Econòmica

Contact information, working papers,
are here

Information about published papers is here

Computer code is at https://github.com/mcreel

MPFM Econometrics

- My current work is related to the 2013 working paper Indirect
Likelihood Inference , with D. Kristensen.

- a new paper, with additional authors and using considerably
improved methods (see links below), will appear later in 2015

- most recent code for the auction model: https://github.com/mcreel/ABCAuction

- this code uses adaptive importance sampling, as described
in algorithms 2 and 3 of the CS&DA paper, below. This
gives very good results using many fewer simulations than in
the 2013 version of the Indirect Likelihood Inference paper.

- the code also uses local linear nonparametric fitting, as
suggested by Beaumont et al. 2002 and Gao and Hong, 2015.
This gives results with about 20% smaller RMSE than local
constant nonparametric fitting.

- most recent code for the DSGE model: https://github.com/mcreel/ABCDGSE

- this code uses the selection of statistics from the
CS&DA paper, importance sampling, and local linear
regression. The code is now follows theory completely, it is
pretty fast, and it seems to work very well

- here are results for the DSGE model described in the
paper, n=160, using a single set of statistics selected
using the method of the CS&DA paper. pmean, etc., psdev,
etc., refer to the prior, without a p refers to the SBIL
estimator:

`Estimation results (LL mean): rep 1000`

true mean pmean sdev. psdev bias pbias rmse prmse

alpha 0.33000 0.33285 0.30000 0.01251 0.05774 0.00285 -0.03000 0.01282 0.06506

beta 0.99000 0.99094 0.97500 0.00269 0.01443 0.00094 -0.01500 0.00285 0.02082

delta 0.02500 0.02605 0.05500 0.00358 0.02598 0.00105 0.03000 0.00373 0.03969

gam 2.00000 2.03608 2.50000 0.09615 1.44338 0.03608 0.50000 0.10265 1.52753

rho1 0.90000 0.89766 0.50000 0.01016 0.28868 -0.00234 -0.40000 0.01042 0.49329

sigma1 0.02000 0.01922 0.05000 0.00134 0.02887 -0.00078 0.03000 0.00155 0.04163

rho2 0.70000 0.68762 0.50000 0.04476 0.28868 -0.01238 -0.20000 0.04642 0.35119

sigma2 0.01000 0.01010 0.05000 0.00150 0.02887 0.00010 0.04000 0.00150 0.04933

nss 0.33333 0.33437 0.31250 0.00378 0.03608 0.00104 -0.02083 0.00392 0.04167

- The older (2011) version with additional examples is here
- An even older working paper was the origin of this work. It
might still have some elements of interest. The VAR example is
based on flawed code, though, ignore it: "A
Data Mining Approach to Indirect Inference"

- On Selection of Statistics for
Approximate Bayesian Computing (or the Method of Simulated
Moments) (with Dennis Kristensen) (2015, Computational
Statistics & Data Analysis)

- code using written in the Julia language for replication of results or using the method with for your own problem is here: https://github.com/mcreel/SelectStatisticsABC
- contact me if you would like the data file for the jump
diffusion example

- I highly recommend getting to know the Julia language - it has a friendly Matlab-like syntax, but speed like C.
- A
Note on Julia and MPI, with Code Examples (2015,
Computational Economics). Simple Monte Carlo, and also ABC
estimation of the auction model in the Indirect Likelihood
Inference paper.

- Econometrics notes Graduate level notes for econometrics, with numerous examples embedded as links. If you use the notes, you'll find the examples ready to run using virtualization software on the accompanying live image. The live image has GNU Octave, Julia and R installed. Most examples use Octave, but I plan to add more examples using Julia. R is included mostly to run other people's code. The code includes methods such as MLE, GMM, kernel estimation, nearest neighbors, parallel computing, etc., etc.
- (
~~no longer in development~~~~, but works perfectly well~~). Development has been taken over by Aissam Hidoussi. Thanks! The amazing PelicanHPC bootable CD to create a HPC cluster for parallel computing in minutes!