Skip to content

Commit

Permalink
Updated all notebooks to use the same run_replicates method.
Browse files Browse the repository at this point in the history
  • Loading branch information
MichaelClerx committed Jul 21, 2021
1 parent 80d181f commit 12986f7
Show file tree
Hide file tree
Showing 8 changed files with 851 additions and 968 deletions.
580 changes: 332 additions & 248 deletions Monomial gamma HMC.ipynb

Large diffs are not rendered by default.

172 changes: 109 additions & 63 deletions differential_evolution.ipynb

Large diffs are not rendered by default.

176 changes: 111 additions & 65 deletions dream.ipynb

Large diffs are not rendered by default.

255 changes: 144 additions & 111 deletions emcee_hammer.ipynb

Large diffs are not rendered by default.

127 changes: 71 additions & 56 deletions haario_ac.ipynb

Large diffs are not rendered by default.

172 changes: 83 additions & 89 deletions haario_bardenet_ac.ipynb

Large diffs are not rendered by default.

287 changes: 0 additions & 287 deletions run-replicates.ipynb

This file was deleted.

50 changes: 1 addition & 49 deletions utils.py
Original file line number Diff line number Diff line change
Expand Up @@ -7,23 +7,7 @@
import pints


def run_replicates(iterations, n_replicates, test):
df = pd.DataFrame(columns=['iterations', 'replicate', 'kld', 'ess'],
index=np.arange(len(iterations) * n_replicates))
k = 0
for it in iterations:
for rep in range(n_replicates):
result = test(it)
df.iloc[k] = {'iterations': it, 'replicate': rep,
'kld': result['kld'], 'ess': result['mean-ess']}
k += 1
df['iterations'] = pd.to_numeric(df['iterations'])
df['kld'] = pd.to_numeric(df['kld'])
df['ess'] = pd.to_numeric(df['ess'])
return df


def run_replicates2(iterations, n_replicates, test, parallel=False):
def run_replicates(iterations, n_replicates, test, parallel=False):
"""
Runs ``test(i)`` for all entries ``i`` in ``iterations``, repeating each
test ``n_replicates`` times.
Expand Down Expand Up @@ -70,35 +54,3 @@ def ecdf_norm_plotter(draws, normal_sd, x=np.linspace(-5, 5, 100)):
plt.ylabel('True cdf')
plt.show()


def run_replicates_distance(iterations, n_replicates, test):
df = pd.DataFrame(columns=['iterations', 'replicate', 'distance', 'ess'],
index=np.arange(len(iterations) * n_replicates))
k = 0
for it in iterations:
for rep in range(n_replicates):
result = test(it)
df.iloc[k] = {'iterations': it, 'replicate': rep,
'distance': result['distance'],
'ess': result['mean-ess']}
k += 1
df['iterations'] = pd.to_numeric(df['iterations'])
df['distance'] = pd.to_numeric(df['distance'])
df['ess'] = pd.to_numeric(df['ess'])
return df


def run_replicates_annulus(iterations, n_replicates, test):
df = pd.DataFrame(columns=['iterations', 'replicate', 'distance', 'ess'],
index=np.arange(len(iterations) * n_replicates))
k = 0
for it in iterations:
for rep in range(n_replicates):
result = test(it)
df.iloc[k] = {'iterations': it, 'replicate': rep, 'distance':
result['distance'], 'ess': result['mean-ess']}
k += 1
df['iterations'] = pd.to_numeric(df['iterations'])
df['distance'] = pd.to_numeric(df['distance'])
df['ess'] = pd.to_numeric(df['ess'])
return df

0 comments on commit 12986f7

Please sign in to comment.