I am trying to figure out what I need to assume in a project work I am writing.
Suppose I have $N$ $0-1$ Bernoulli trials $\{B_i\}$ with success probabilities $p_i$. I know that they are positively correlated since I know their conditional distribution given another Bernoulli variable $X$ (and given my particular parameters this indeed happens). We have $B_i|X=^d Be(q_i)$ where $q_i$ is computable/known. I know that in general, the joint distribution of an $N$-dimensional Bernoulli vector is determined by the "moments" $E(B_iB_j...B_k)$ (see Multivariate normal and multivariate Bernoulli, Joint distribution of dependent Bernoulli Random variables). However, one should think that knowing the conditional distribution is quite a bit stronger than just knowing the second moment. Is there any way I can extract more information from conditional distribution above about the joint distribution? If not, what sort of "information" would I require to compute these mixed moments?
Any help or reference would be much appreciated!
EDIT: Would it for instance suffice to know the joint distribution of $\{B_i\}$ conditional on $X$?