Constraining Output of Multinomial Distribution

[mathDR]

linking a different discussion that has an analytical integral solution to the Multivariate Wallenius Noncentral Hypergeometric distribution

Biased urn problem & Wallenius' noncentral hypergeometric distribution

I played around and got the integrate_1d to work for this distribution. No doubt that this is way slower than @mathDR approximation. This is my first time using the integrate functionality so that was cool! I’m sure there’s a more numerically stable way to code it (@andrjohns or @bbbales2 may know a good way to use the xc parameter here). functions { real wallenius_integral(real t, // Function argument real xc, // Complement of function argument …