Formule de probabilité pour une distribution bernoulli multivariée

J'ai besoin d'une formule pour la probabilité d'un événement dans une distribution de Bernoulli à n variables $X\in\{0,1\}^n$ avec des probabilités $P(X_i=1)=p_i$ pour un seul élément et pour des paires d'éléments $P(X_i=1 \wedge X_j=1)=p_{ij}$ . Je pourrais vous donner de manière équivalente moyenne et covariance de $X$ .

J'ai déjà appris qu'il existe de nombreuses distributions $\{0,1\}^n$ ayant les propriétés tout comme il existe de nombreuses distributions ayant une moyenne et une covariance données. Je recherche une canonique sur $\{0,1\}^n$ , tout comme la gaussienne est une distribution canonique pour $R^n$ et une moyenne et une covariance données.

La variable aléatoire prenant des valeurs dans est une variable aléatoire discrète. Sa distribution est entièrement décrite par les probabilités p i = P ( X = i ) avec i ∈ { 0 , 1 } n . Les probabilités p i et p i j que vous donnez sont des sommes de p i pour certains indices i .$\{0,1\}^n$$p_{\mathbf{i}}=P(X=\mathbf{i})$$\mathbf{i}\in\{0,1\}^n$$p_{i}$$p_{ij}$$p_{\mathbf{i}}$$\mathbf{i}$

Il semble maintenant que vous souhaitiez décrire en utilisant uniquement p i et p i j . Ce n'est pas possible sans supposer certaines propriétés sur p i . Pour voir que essayer de tirer fonction caractéristique de X . Si nous prenons n = 3, nous obtenons$p_{\mathbf{i}}$$p_i$$p_{ij}$$p_{\mathbf{i}}$$X$$n=3$

Il n'est pas possible de réorganiser cette expression pour quep

$$\begin{align} Ee^{i(t_1X_1+t_2X_2+t_3X_3)}&=p_{000}+p_{100}e^{it_1}+p_{010}e^{it_2}+p_{001}e^{it_3}\\\\ &+p_{110}e^{i(t_1+t_2)}+p_{101}e^{i(t_1+t_3)}+p_{011}e^{i(t_2+t_3)}+p_{111}e^{i(t_1+t_2+t_3)} \end{align}$$ $p_{\mathbf{i}}$ dissapear. For the gaussian random variable the characteristic function depends only on mean and covariance parameters. Characteristic functions uniquely define distributions, so this is why Gaussian can be described uniquely by using only mean and covariance. As we see for random variable $X$ this is not the case.

 

See the following paper:

J. L. Teugels, Some representations of the multivariate Bernoulli and binomial distributions, Journal of Multivariate Analysis, vol. 32, no. 2, Feb. 1990, 256–268.

Here is the abstract:

Des versions multivariées mais vectorisées pour les distributions de Bernoulli et binomiales sont établies en utilisant le concept de produit Kronecker du calcul matriciel. La distribution de Bernoulli multivariée implique un modèle paramétré, qui fournit une alternative au modèle log-linéaire traditionnel pour les variables binaires.

I don't know what the resulting distribution is called, or if it even has a name, but it strikes me the obvious way to set this up is to think of the model you'd use to model a 2×2×2×…×2 table using a log-linear (Poisson regression) model. As you know the 1st-order interactions only, it's then natural to assume that all higher-order interactions are zero.

Using the questioner's notation, this gives the model:

$$P(X_1=x_1, X_2=x_2,\ldots,X_n=x_n) = \prod_i \left[ p_i^{x_i}(1-p_i)^{1-x_i} \prod_{j<i} \left(\frac{p_{ij}}{p_i p_j}\right)^{x_i x_j} \right]$$