{"action":"created","issue":{"url":"https://api.github.com/repos/ligo-cbc/pycbc/issues/793","repository_url":"https://api.github.com/repos/ligo-cbc/pycbc","labels_url":"https://api.github.com/repos/ligo-cbc/pycbc/issues/793/labels{/name}","comments_url":"https://api.github.com/repos/ligo-cbc/pycbc/issues/793/comments","events_url":"https://api.github.com/repos/ligo-cbc/pycbc/issues/793/events","html_url":"https://github.com/ligo-cbc/pycbc/pull/793","id":151684527,"number":793,"title":"Add inference module","user":{"login":"cdcapano","id":8258446,"avatar_url":"https://avatars.githubusercontent.com/u/8258446?v=3","gravatar_id":"","url":"https://api.github.com/users/cdcapano","html_url":"https://github.com/cdcapano","followers_url":"https://api.github.com/users/cdcapano/followers","following_url":"https://api.github.com/users/cdcapano/following{/other_user}","gists_url":"https://api.github.com/users/cdcapano/gists{/gist_id}","starred_url":"https://api.github.com/users/cdcapano/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/cdcapano/subscriptions","organizations_url":"https://api.github.com/users/cdcapano/orgs","repos_url":"https://api.github.com/users/cdcapano/repos","events_url":"https://api.github.com/users/cdcapano/events{/privacy}","received_events_url":"https://api.github.com/users/cdcapano/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":{"login":"cmbiwer","id":10871764,"avatar_url":"https://avatars.githubusercontent.com/u/10871764?v=3","gravatar_id":"","url":"https://api.github.com/users/cmbiwer","html_url":"https://github.com/cmbiwer","followers_url":"https://api.github.com/users/cmbiwer/followers","following_url":"https://api.github.com/users/cmbiwer/following{/other_user}","gists_url":"https://api.github.com/users/cmbiwer/gists{/gist_id}","starred_url":"https://api.github.com/users/cmbiwer/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/cmbiwer/subscriptions","organizations_url":"https://api.github.com/users/cmbiwer/orgs","repos_url":"https://api.github.com/users/cmbiwer/repos","events_url":"https://api.github.com/users/cmbiwer/events{/privacy}","received_events_url":"https://api.github.com/users/cmbiwer/received_events","type":"User","site_admin":false},"milestone":null,"comments":1,"created_at":"2016-04-28T16:18:54Z","updated_at":"2016-04-28T17:31:49Z","closed_at":null,"pull_request":{"url":"https://api.github.com/repos/ligo-cbc/pycbc/pulls/793","html_url":"https://github.com/ligo-cbc/pycbc/pull/793","diff_url":"https://github.com/ligo-cbc/pycbc/pull/793.diff","patch_url":"https://github.com/ligo-cbc/pycbc/pull/793.patch"},"body":"This adds an \"inference\" module to pycbc, which contains likelihood.py. This file stores classes for estimating likelihood for parameter estimation. Currently, there are two classes, a _BaseLikelihoodEvaluator and a GaussianLikelihood, which inherits from _Base. The base class sets up the general structure, which is to have a data attribute, a waveform_generator attribute, etc. GaussianLikelihood adds a log_likelihood function to this which computes the standard likelihood for Gaussian noise. No marginalization is done, nor is there any PSD fitting. When we want likelihood evaluators that do more fancy things like PSD fitting and/or some internal marginalization (e.g., phase), the intent is these can be added as new classes that inherit from the _Base class and add their own loglikelihood function.\r\n\r\nSince the likelihood evaluator just expects an instance of a waveform generator, it is agnostic to the types of waveforms that are being estimated. It can therefore be used for CBC PE, or ringdown PE, or any other PE for which a generator exists (sine Gaussians, say).\r\n\r\nAn example usage is included in the GaussianLikelihood doc string. The plot created in that example can be found here:\r\n\r\nhttps://www.atlas.aei.uni-hannover.de/~cdcapano/likelihood_eval_test.png\r\n\r\nThe peak in the log likelihood occurs at 3.1s, which is where the example signal was placed."},"comment":{"url":"https://api.github.com/repos/ligo-cbc/pycbc/issues/comments/215503553","html_url":"https://github.com/ligo-cbc/pycbc/pull/793#issuecomment-215503553","issue_url":"https://api.github.com/repos/ligo-cbc/pycbc/issues/793","id":215503553,"user":{"login":"cdcapano","id":8258446,"avatar_url":"https://avatars.githubusercontent.com/u/8258446?v=3","gravatar_id":"","url":"https://api.github.com/users/cdcapano","html_url":"https://github.com/cdcapano","followers_url":"https://api.github.com/users/cdcapano/followers","following_url":"https://api.github.com/users/cdcapano/following{/other_user}","gists_url":"https://api.github.com/users/cdcapano/gists{/gist_id}","starred_url":"https://api.github.com/users/cdcapano/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/cdcapano/subscriptions","organizations_url":"https://api.github.com/users/cdcapano/orgs","repos_url":"https://api.github.com/users/cdcapano/repos","events_url":"https://api.github.com/users/cdcapano/events{/privacy}","received_events_url":"https://api.github.com/users/cdcapano/received_events","type":"User","site_admin":false},"created_at":"2016-04-28T17:31:49Z","updated_at":"2016-04-28T17:31:49Z","body":"Crap... no, we still have to support 2.6. Sorry, forgot about that. I'll fix it in the second bit. I'm going to leave the example string as-is though, since that's just a doc string."}}