Skip to content

Changes

Page history
Update dynesty authored by Gregory Ashton's avatar Gregory Ashton
Show whitespace changes
Inline Side-by-side
pp/dynesty.md
View page @ 2342152e
...@@ -2,141 +2,56 @@ ...@@ -2,141 +2,56 @@
All runs can be found [on cit](https://ldas-jobs.ligo.caltech.edu/~gregory.ashton/bilby_tests/dynesty/). All runs can be found [on cit](https://ldas-jobs.ligo.caltech.edu/~gregory.ashton/bilby_tests/dynesty/).
Runs are generated by drawing simulated signals in Gaussian (aLIGO noise curve) data for H1 and L1. ### runC
* bilby_pipe v0.0.4 (ddf3a61)
* bilby v0.4.4 (57ce7f43)
Then, for the analysis the same prior is used **except** that the `geocent_time` prior is set to * dynesty v0.9.5.3
```
geocent_time = Uniform(minimum=-0.05, maximum=0.05, name='geocent_time', latex_label='$t_c$', unit='$s$')
```
in all cases we use the `bilby.gw.likelihood.GravitationalWaveTransientLikelihood` with no marginalisation.
## Run 1: `nlive=2000` and `walks=50` with no spins
For the first run, we do a no-spin run, i.e. the prior is given by
```python
chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100, unit='$M_{\\odot}$')
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
a_1 = 0
a_2 = 0
tilt_1 = 0
tilt_2 = 0
phi_12 = 0
phi_jl = 0
dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi)
iota = Sine(name='iota')
psi = Uniform(name='psi', minimum=0, maximum=np.pi)
phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi)
luminosity_distance = PowerLaw(alpha=2, minimum=100, maximum=5000, name='luminosity_distance')
geocent_time = 0
```
Here is the full set of sampler kwargs:
```python ```python
{'bound': 'multi', 'sample': 'rwalk', 'verbose': True, 'check_point_delta_t': 600, 'nlive': 2000, 'first_update': None, 'npdim': None, 'rstate': None, 'queue_size': None, 'pool': None, 'use_pool': None, 'live_points': None, 'logl_args': None, 'logl_kwargs': None, 'ptform_args': None, 'ptform_kwargs': None, 'enlarge': None, 'bootstrap': None, 'vol_dec': 0.5, 'vol_check': 2.0, 'facc': 0.5, 'slices': 5, 'walks': 50, 'update_interval': 1200, 'print_func': <bound method Dynesty._print_func of <bilby.core.sampler.dynesty.Dynesty object at 0x2b3c0f61c438>>, 'dlogz': 0.1, 'maxiter': None, 'maxcall': None, 'logl_max': inf, 'add_live': True, 'print_progress': True, 'save_bounds': True} # config.ini
``` label = dynesty_pp_test
outdir = .
### PP plot accounting = ligo.dev.o3.cbc.pe.lalinference
![image](uploads/354c3d9e27257a476f2f8cee838bb1da/image.png) detectors = [H1, L1]
duration = 4
### Sampling time trigger-time = 0
deltaT = 0.1
prior-file = chirp_mass_mass_ratio_spins.prior
![image](uploads/11b36aeaa61d598be5ff8b122862a36c/image.png) sampler = dynesty
sampler-kwargs = {nlive: 1000, walks: 100, n_check_point: 5000}
### SNR distribition create-plots = False
injection = True
n-injection = 500
![image](uploads/966793dc31433597902bd919570192cc/image.png) time-marginalization = True
distance-marginalization = True
## Run 2: `nlive=2000` and `walks=50` with spins phase-marginalization = True
```
```python ```python
# chirp_mass_mass_ratio_spins.prior
chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100, unit='$M_{\\odot}$') chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100, unit='$M_{\\odot}$')
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1) mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
a_1 = Uniform(name='a_1', minimum=0, maximum=0.8) a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
a_2 = Uniform(name='a_2', minimum=0, maximum=0.8) a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
tilt_1 = Sine(name='tilt_1') tilt_1 = Sine(name='tilt_1')
tilt_2 = Sine(name='tilt_2') tilt_2 = Sine(name='tilt_2')
phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi) phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, periodic_boundary=True)
phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi) phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, periodic_boundary=True)
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, periodic_boundary=True)
dec = Cosine(name='dec') dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi)
iota = Sine(name='iota')
psi = Uniform(name='psi', minimum=0, maximum=np.pi)
phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi)
luminosity_distance = PowerLaw(alpha=2, minimum=100, maximum=5000, name='luminosity_distance') luminosity_distance = PowerLaw(alpha=2, minimum=100, maximum=5000, name='luminosity_distance')
geocent_time = 0 theta_jn = Sine(name='iota')
``` geocent_time = Uniform(name='geocent_time', minimum=-0.05, maximum=0.05)
psi = Uniform(name='psi', minimum=0, maximum=np.pi, periodic_boundary=True)
The set of kwargs are phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, periodic_boundary=True)
```
{'bound': 'multi', 'sample': 'rwalk', 'verbose': True, 'check_point_delta_t': 600, 'nlive': 2000, 'first_update': None, 'npdim': None, 'rstate': None, 'queue_size': None, 'pool': None, 'use_pool': None, 'live_points': None, 'logl_args': None, 'logl_kwargs': None, 'ptform_args': None, 'ptform_kwargs': None, 'enlarge': None, 'bootstrap': None, 'vol_dec': 0.5, 'vol_check': 2.0, 'facc': 0.5, 'slices': 5, 'walks': 50, 'update_interval': 1200, 'print_func': <bound method Dynesty._print_func of <bilby.core.sampler.dynesty.Dynesty object at 0x2b07df2b4358>>, 'dlogz': 0.1, 'maxiter': None, 'maxcall': None, 'logl_max': inf, 'add_live': True, 'print_progress': True, 'save_bounds': True}
```
### PP plot
![image](uploads/3acfb087ead9680b9c69c82ba55d961c/image.png)
### Sampling time
![image](uploads/407be4a9d232a86801a0146dd38e39b3/image.png)
### SNR distribition
![image](uploads/e5278c2e1610a1a3fde624e14ce4e2e9/image.png)
## Run 3: `nlive=3000` and `walks=50` with spins
Same prior as in run2, but with more live points
### PP plot
![image](uploads/029f0642e3740af8922b80f55a000438/image.png)
### Sampling time
![image](uploads/a09c6a15dc608053eb45ac033126b1a4/image.png)
### SNR distribution
![image](uploads/aa368dc0f3979c8f831428dea96c5604/image.png)
## Old runs
The following are some old results, leaving here in case we want to come back to them.
### Run 1
Using this kwargs argument for `dynesty`
``` ```
{'bound': 'multi', 'sample': 'rwalk', 'verbose': True, 'check_point_delta_t': 600,
'nlive': 1024, 'first_update': None, 'npdim': None, 'rstate': None, 'queue_size': None,
'pool': None, 'use_pool': None, 'live_points': None, 'logl_args': None,
'logl_kwargs': None, 'ptform_args': None, 'ptform_kwargs': None, 'enlarge': None,
'bootstrap': None, 'vol_dec': 0.5, 'vol_check': 2.0, 'facc': 0.5, 'slices': 5, 'walks': 60,
'update_interval': 614, 'print_func': <bound method Dynesty._print_func of <bilby.core.sampler.dynesty.Dynesty object at 0x2aec680cfef0>>,
'dlogz': 0.1, 'maxiter': None, 'maxcall': None, 'logl_max': inf, 'add_live': True,
'print_progress': True, 'save_bounds': True}
```
Produced this PP plot:
![image](uploads/0cd001c8e5cabf955a2e94a296814e62/image.png)
### Run 2
Running again with more walks
```
{'bound': 'multi', 'sample': 'rwalk', 'verbose': True, 'check_point_delta_t': 600, 'nlive': 1024, 'first_update': None, 'npdim': None, 'rstate': None, 'queue_size': None, 'pool': None, 'use_pool': None, 'live_points': None, 'logl_args': None, 'logl_kwargs': None, 'ptform_args': None, 'ptform_kwargs': None, 'enlarge': None, 'bootstrap': None, 'vol_dec': 0.5, 'vol_check': 2.0, 'facc': 0.5, 'slices': 5, 'walks': 120, 'update_interval': 614, 'print_func': <bound method Dynesty._print_func of <bilby.core.sampler.dynesty.Dynesty object at 0x2b2c0a83aac8>>, 'dlogz': 0.1, 'maxiter': None, 'maxcall': None, 'logl_max': inf, 'add_live': True, 'print_progress': True, 'save_bounds': True}
```
![image](uploads/923eba9fdfde10393dad11c30b0e7086/image.png)
![image](uploads/c86ad17839756e84b804f98db07c836e/image.png)
![image](uploads/c7ded58e8db76251cce48bae297c192a/image.png)
![image](uploads/4ab487489ea384d4f80fdcdaf85f2f3b/image.png)