Add an interface for joint priors.

Only available for the new estimation syntax. Something like:

[alpha, beta].prior(shape=gaussian, mean=Vector, variance=Matrix, ...)

This interface is needed for Dirichlet priors over probabilities.