multivariable distributions - T and P separate or together

[jgcri]

multivariable distributions - T and P separate or together

[abigailsnyder]

per @claudiatebaldi would expect jointly estimated and jointly sampled.

@abigailsnyder will

• update the code in data_raw to jointly estimate - 24 parameters to get a joint multivariable beta distribution of T and P fractions. And add more documentation to those functions.
• refit the models.
• open a PR

Then go back into the monthly_downscaling code and update sampling to be joint, as well as adding options outlined in #16

[abigailsnyder]

In terms of updating the sampling to be joint, it looks like the separate sampling for each variable is happening in the cassandra components code:
https://github.com/JGCRI/cassandra/blob/master/cassandra/components.py Lines 964-977

Which explains why it's harder to tell from the R monthly_downscaling sampling code that the variables are being treated separately than in the data_raw/... training code.

So the R code will have to be updated for the sampling but then the cassandra code will also have to be updated, FYI @crvernon