search

nanograv / enterprise

ENTERPRISE (Enhanced Numerical Toolbox Enabling a Robust PulsaR Inference SuitE) is a pulsar timing analysis code, aimed at noise analysis, gravitational-wave searches, and timing model analysis.
https://enterprise.readthedocs.io
MIT License
65 stars 64 forks source link

Towards distinguishing CW sources on GWB #53

Open paulthebaker opened 7 years ago

paulthebaker commented 7 years ago

During the astrophysics busyweek, we discussed a paper idea to study the prospects for distinguishing CW sources from the stochastic background.

For phase 2 of this, I would like to use an RJMcMC to adaptively choose the number of CW sources (by frequency and sky location). This could even get crazy and adaptively fit an unknown number anisotropic multipole moments and CW sources.

I'm opening this issue to make sure we think about this sort of functionality as we move forward. I don't know of any out of the box, widely used RJ samplers. @jellis18, your RJSampler from PAL2 is BayesWave inspired, right? It could handle this sort of thing.

paulthebaker commented 7 years ago

Can we use the operator definitions to do RJ moves like:

model = gwb + cw(1)
...
model += cw(2)
...
model -= cw(1)
...

to add and remove CW instances from the model?