# Combining multiple posterior distributions

I am new to Bayesian statistics, and thus have problems to come up with a solution for the following problem:

Using Approximate Bayesian Computation (ABC), I generate a posterior distribution from one set of observed summary statistics, and a set of about a million simulated summary statistics. The prior distribution is the variable value of the model parameter used in the simulations. The posterior distribution is used to obtain point estimates (e.g., mode, mean, etc.) of the underlying model parameter.

The biological problem I am working on is that the true parameter value cannot be estimated from one observation only. Hence, I have N observed data sets, and obtain N posterior distributions. In ABC, I use the same simulated data for each observation data set.

Now, my problem is to combine the N posterior distributions in a way to estimate the true parameter value, thus capturing the information contained in all N densities.

Any ideas?

I think you are getting the problem wrong and this has nothing to do with ABC. If your posterior distribution is such that the parameter $\theta$ is not entirely identifiable from one dataset, $x_1$, but is identifiable from a $N$ sample, $x_1,\ldots,x_N$, what you need to consider is the "single" posterior conditional on all of those datasets: $$\pi(\theta|x_1,\ldots,x_N)$$ not a product of posteriors (which is not defined from a Bayesian viewpoint anyway).
1. Generate theta from the prior

because your data is then the $N$ dataset.