I have a set of differential equations that describes a physical system - actually a ball being heated. I also have some data taken from real experiments on the system. I would like to try to fit these parameters to the data using Bayesian inference. Over the past few days I've been getting to grips with MCMC and I now have something that seems to work well.
The graph below shows both the both the data (in pink) and the fitted model with credible intervals.
Now I would like to repeat this experiment with a different physical sample in the test machine. The posterior distribution from the first test is my best guess at how the new sample will behave, but after three or four measurements I'll know a lot more about how the new sample behaves. This is where I am uncertain what to do. I think I should use the first posterior as my new prior, but I'm not sure how to do this. I think my options are:
- Fit an analytical distribution (or ensemble of distributions) to the posterior data
- Carry all of the data from both experiments through to the MCMC algorithm
Is this the correct approach? For the second option I'm worried that carrying all of the data from the first experiment will unduly skew the results for the second experiment - broadly speaking I know the physical sample from the second experiment will be different to the first, so I feel that I should be weighting this data more highly. I could do this by scaling the sum of squares error function to make errors from the second experiment more important than first, but I wondered if there was a more correct way of doing this?