Modeling bivariate beta distributions

Yes, this is exactly why the model denoted model_simple in my previous reply is not able to catch the correlation. But isn’t the model from the paper designed to actually have parameters for correlation in it?