Skip to content
Snippets Groups Projects
Commit ad316e65 authored by Xisco Jimenez Forteza's avatar Xisco Jimenez Forteza
Browse files

files

parent 48c032b4
No related branches found
No related tags found
No related merge requests found
RDanalysis.nb 0 → 100644
This diff is collapsed.
RDown.m 0 → 100644
This diff is collapsed.
RDown.nb 0 → 100644
This diff is collapsed.
...@@ -2,3 +2,8 @@ ...@@ -2,3 +2,8 @@
Domain dedicated to the Ringdown stacking project. All code, files and plots should be reasonably documented in the README file. Domain dedicated to the Ringdown stacking project. All code, files and plots should be reasonably documented in the README file.
RDown.m: RDown package file that loads the functions needed for our computation.
RDown.nd: RDown mathematica file to add/modify our RDown functions. Every time you modify and save this file, the RDown.m package is automatically updated.
RDanalysis.nb: Mathematica notebook used to generate our results.
mcmc.m: Mathematica package to run mcmc samplers.
l2: Folder containing the modes and tones data. It is taken from Berti & Cardoso.
Mathematica Markov Chain Monte Carlo
====================================
Mathematica package containing a general-purpose [Markov chain Monte Carlo](http://en.wikipedia.org/wiki/Markov_chain_Monte_Carlo) routine I wrote. Includes various examples and documentation.
### Features:
* Convenience wrapper for fitting models to arbitrary-dimensional data with Gaussian errors
* Handles both real-valued and discrete-valued model parameters
* Uses Metropolis algorithm with decaying exponential proposal distribution
* Progress monitor; support for auto save/resume
### Files:
* `mcmc.m`: Package file
* `mcmc_demonst.nb`: Demonstrations and documentation
For a good Python MCMC implementation, check out [emcee](http://dan.iel.fm/emcee).
2 2 0 1.5251d0 -1.1568d0 0.1292d0 0.7000d0 1.4187d0 -0.4990d0
2 2 1 1.3673d0 -1.0260d0 0.1628d0 0.1000d0 0.5436d0 -0.4731d0
2 2 2 1.3223d0 -1.0257d0 0.1860d0 -0.1000d0 0.4206d0 -0.4256d0
2 1 0 0.6000d0 -0.2339d0 0.4175d0 -0.3000d0 2.3561d0 -0.2277d0
2 1 1 0.5800d0 -0.2416d0 0.4708d0 -0.3300d0 0.9501d0 -0.2072d0
2 1 2 0.5660d0 -0.2740d0 0.4960d0 -0.1000d0 0.4173d0 -0.2774d0
2 0 0 0.4437d0 -0.0739d0 0.3350d0 4.0000d0 -1.9550d0 0.1420d0
2 0 1 0.4185d0 -0.0768d0 0.4355d0 1.2500d0 -0.6359d0 0.1614d0
2 0 2 0.3734d0 -0.0794d0 0.6306d0 0.5600d0 -0.2589d0 0.3034d0
2 -1 0 0.3441d0 0.0293d0 2.0010d0 2.0000d0 0.1078d0 5.0069d0
2 -1 1 0.3165d0 0.0301d0 2.3415d0 0.6100d0 0.0276d0 13.1683d0
2 -1 2 0.2696d0 0.0315d0 2.7755d0 0.2900d0 0.0276d0 6.4715d0
2 -2 0 0.2938d0 0.0782d0 1.3546d0 1.6700d0 0.4192d0 1.4700d0
2 -2 1 0.2528d0 0.0921d0 1.3344d0 0.4550d0 0.1729d0 1.3617d0
2 -2 2 0.1873d0 0.1117d0 1.3322d0 0.1850d0 0.1266d0 1.3661d0
3 3 0 1.8956d0 -1.3043d0 0.1818d0 0.9000d0 2.3430d0 -0.4810d0
3 3 1 1.8566d0 -1.2818d0 0.1934d0 0.2274d0 0.8173d0 -0.4731d0
3 3 2 1.8004d0 -1.2558d0 0.2133d0 0.0400d0 0.5445d0 -0.4539d0
3 2 0 1.1481d0 -0.5552d0 0.3002d0 0.8313d0 2.3773d0 -0.3655d0
3 2 1 1.1226d0 -0.5471d0 0.3264d0 0.2300d0 0.8025d0 -0.3684d0
3 2 2 1.0989d0 -0.5550d0 0.3569d0 0.1000d0 0.4804d0 -0.3784d0
3 1 0 0.8345d0 -0.2405d0 0.4095d0 23.8450d0 -20.7240d0 0.03837d0
3 1 1 0.8105d0 -0.2342d0 0.4660d0 8.8530d0 -7.8506d0 0.03418d0
3 1 2 0.7684d0 -0.2252d0 0.5805d0 2.1800d0 -1.6273d0 0.1136d0
3 0 0 0.6873d0 -0.09282d0 0.3479d0 6.7841d0 -3.6112d0 0.09480d0
3 0 1 0.6687d0 -0.09155d0 0.4021d0 2.0075d0 -0.9930d0 0.12297d0
3 0 2 0.6343d0 -0.08915d0 0.5117d0 0.9000d0 -0.3409d0 0.2679d0
3 -1 0 0.5751d0 0.02508d0 3.1360d0 3.0464d0 0.1162d0 -0.2812d0
3 -1 1 0.5584d0 0.02514d0 3.4154d0 1.2000d0 -0.1928d0 0.1037d0
3 -1 2 0.5271d0 0.02561d0 3.8011d0 1.0000d0 -0.4424d0 0.02467d0
3 -2 0 0.5158d0 0.08195d0 1.4084d0 2.9000d0 0.3356d0 2.3050d0
3 -2 1 0.4951d0 0.08577d0 1.4269d0 0.9000d0 0.1295d0 1.6142d0
3 -2 2 0.4567d0 0.09300d0 1.4469d0 0.4900d0 0.0848d0 1.9737d0
3 -3 0 0.4673d0 0.1296d0 1.3255d0 2.5500d0 0.6576d0 1.3378d0
3 -3 1 0.4413d0 0.1387d0 1.3178d0 0.7900d0 0.2381d0 1.3706d0
3 -3 2 0.3933d0 0.1555d0 1.3037d0 0.4070d0 0.1637d0 1.3819d0
4 4 0 2.3000d0 -1.5056d0 0.2244d0 1.1929d0 3.1191d0 -0.4825d0
4 4 1 2.3000d0 -1.5173d0 0.2271d0 0.3000d0 1.1034d0 -0.4703d0
4 4 2 2.3000d0 -1.5397d0 0.2321d0 0.1100d0 0.6997d0 -0.4607d0
4 3 0 1.6869d0 -0.8862d0 0.2822d0 1.4812d0 2.8096d0 -0.4271d0
4 3 1 1.6722d0 -0.8843d0 0.2923d0 0.4451d0 0.9569d0 -0.4250d0
4 3 2 1.6526d0 -0.8888d0 0.3081d0 0.2200d0 0.5904d0 -0.4236d0
4 2 0 1.2702d0 -0.4685d0 0.3835d0 -3.6000d0 7.7749d0 -0.1491d0
4 2 1 1.2462d0 -0.4580d0 0.4139d0 -1.5000d0 2.8601d0 -0.1392d0
4 2 2 1.2025d0 -0.4401d0 0.4769d0 -1.5000d0 2.2784d0 -0.1124d0
4 1 0 1.0507d0 -0.2478d0 0.4348d0 14.0000d0 -9.8240d0 0.09047d0
4 1 1 1.0337d0 -0.2439d0 0.4695d0 4.2000d0 -2.8399d0 0.1081d0
4 1 2 1.0019d0 -0.2374d0 0.5397d0 2.2000d0 -1.4195d0 0.1372d0
4 0 0 0.9175d0 -0.1144d0 0.3511d0 7.0000d0 -2.7934d0 0.1708d0
4 0 1 0.9028d0 -0.1127d0 0.3843d0 2.2000d0 -0.8308d0 0.2023d0
4 0 2 0.8751d0 -0.1096d0 0.4516d0 1.2000d0 -0.4159d0 0.2687d0
4 -1 0 0.7908d0 0.02024d0 5.4628d0 4.6000d0 -0.4038d0 0.4629d0
4 -1 1 0.7785d0 0.02005d0 5.8547d0 1.6000d0 -0.2323d0 0.2306d0
4 -1 2 0.7549d0 0.01985d0 6.5272d0 1.6000d0 -0.8136d0 0.03163d0
4 -2 0 0.7294d0 0.07842d0 1.5646d0 4.0000d0 0.2777d0 2.0647d0
4 -2 1 0.7154d0 0.07979d0 1.5852d0 1.3200d0 0.08694d0 4.3255d0
4 -2 2 0.6885d0 0.08259d0 1.6136d0 0.7500d0 0.05803d0 3.7971d0
4 -3 0 0.6728d0 0.1338d0 1.3413d0 3.700d0 0.5829d0 1.6681d0
4 -3 1 0.6562d0 0.1377d0 1.3456d0 1.1800d0 0.2111d0 1.4129d0
4 -3 2 0.6244d0 0.1454d0 1.3513d0 0.6600d0 0.1385d0 1.3742d0
4 -4 0 0.6256d0 0.1800d0 1.3218d0 3.4000d0 0.8696d0 1.4074d0
4 -4 1 0.6061d0 0.1869d0 1.3168d0 1.0800d0 0.3095d0 1.3279d0
4 -4 2 0.5686d0 0.2003d0 1.3068d0 0.5980d0 0.2015d0 1.3765d0
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment