Sampling from two-dimensional distribution

Thanks!

I have some data, and I want to use these data to infer the parameters in this joint distribution.


Actually I don’t know how to use pymc to infer the parameters from a two-dimensional joint distribution.

So I want to start with a simple case that how to infer the two-dimensional normal distribution using generating data firstly.