I understand. Thanks for your answer.
I still have two questions:
- How to jointly? Is it possible to embed the model in both the first and second ways (jointly and specify a reasonable maximum)?
- Can I include the depth and velocity dimensions as regularization priors in the coefficient matrix? If so how can I implement this in pymc?