# Drawing from posterior of a Multivariate Distribution

**URL:** <https://discourse.pymc.io/t/drawing-from-posterior-of-a-multivariate-distribution/2335>\
**Category:** Questions\
**Created:** [December 6, 2018, 1:10pm UTC](https://discourse.pymc.io/t/drawing-from-posterior-of-a-multivariate-distribution/2335 "2018-12-06T13:10:47Z")\
**Posts on this page:** 1\
**Showing post:** 3

<div class="post-metadata">

**Author:** ![junpenglao](https://yyz2.discourse-cdn.com/flex036/user_avatar/discourse.pymc.io/junpenglao/32/8_2.png) [@junpenglao](https://discourse.pymc.io/u/junpenglao)\
**Post date:** [December 6, 2018, 2:34pm UTC](https://discourse.pymc.io/t/drawing-from-posterior-of-a-multivariate-distribution/2335/3 "2018-12-06T14:34:17Z")

</div>

Oh maybe try something like:

> [@Performance speedup for updating posterior with new data](https://discourse.pymc.io/t/performance-speedup-for-updating-posterior-with-new-data/1097/2):
>
> I had a look at your code. Unfortunately, I dont see much more room to improve the speed, as the Interpolation rv always need to reinitialized, and the bottom neck in the evaluation of scipy Spline smooth function is also difficult to improve without rewriting everything in theano. Alternatively, if you are happy with using some other distribution as an approximation, it would simplify the code a lot: new\_data = pd.DataFrame(data=[generate\_data() for \_ in range(1)], columns=generate\_data().key…

---

_[View the full topic](https://discourse.pymc.io/t/drawing-from-posterior-of-a-multivariate-distribution/2335)._