Add example of adding responses

examples/other_examples/add_multiple_results.py

+#!/bin/python
+"""
+An example of running two sets of posterior sample estimations and adding them
+"""
+from __future__ import division
+import tupak
+import numpy as np
+
+outdir = 'outdir'
+
+
+def model(time, m, c):
+    return time * m + c
+
+
+injection_parameters = dict(m=0.5, c=0.2)
+sigma = 1
+sampling_frequency = 10
+time_duration = 10
+time = np.arange(0, time_duration, 1/sampling_frequency)
+N = len(time)
+data = model(time, **injection_parameters) + np.random.normal(0, sigma, N)
+
+likelihood = tupak.core.likelihood.GaussianLikelihood(
+    time, data, model, sigma=sigma)
+
+priors = {}
+priors['m'] = tupak.core.prior.Uniform(0, 1, 'm')
+priors['c'] = tupak.core.prior.Uniform(-2, 2, 'c')
+
+resultA = tupak.run_sampler(
+    likelihood=likelihood, priors=priors, sampler='emcee', nsteps=1000,
+    nburn=500, outdir=outdir, label='A')
+
+resultB = tupak.run_sampler(
+    likelihood=likelihood, priors=priors, sampler='emcee', nsteps=1000,
+    nburn=500, outdir=outdir, label='B')
+
+resultA.plot_walkers()
+result = resultA + resultB
+result.plot_corner()
+print(result)
+
+