Skip to content
Snippets Groups Projects
Commit 2b67a284 authored by Sylvia Biscoveanu's avatar Sylvia Biscoveanu Committed by Moritz
Browse files

Change example to reflect new structure of CorrelatedUniform

parent 62f575c6
No related branches found
No related tags found
1 merge request!332Resolve "Introduce conditional prior sets"
......@@ -24,7 +24,7 @@ import bilby.gw.prior
def correlation_func_a(mu, a=0):
return mu + a**2 + 2*a + 3
return mu + a**2 + 2 * a + 3
def correlation_func_b(mu, a=0, b=0):
......@@ -44,14 +44,20 @@ corner.corner(np.array(samples))
plt.show()
def correlation_func_min_max(extrema_dict, a, b):
maximum = extrema_dict['maximum'] + a**b
minimum = np.log(b)
return minimum, maximum
a = bilby.core.prior.Uniform(minimum=0, maximum=1)
b = bilby.core.prior.CorrelatedUniform(minimum=0, maximum=1, correlation_func=correlation_func_a)
c = bilby.core.prior.CorrelatedUniform(minimum=0, maximum=1, correlation_func=correlation_func_b)
b = bilby.core.prior.Uniform(minimum=1e-6, maximum=1e-1)
c = bilby.core.prior.CorrelatedUniform(minimum=0, maximum=1, correlation_func=correlation_func_min_max)
correlated_uniform = bilby.core.prior.CorrelatedPriorDict(dictionary=dict(a=a, b=b, c=c))
samples = correlated_uniform.sample(1000000)
pdb.set_trace()
samples = np.array([samples['a'], samples['b'], samples['c']]).T
corner.corner(np.array(samples))
plt.show()
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment