Say I have the following sample from a continuous variable $X$ and a categorical (dichotomous) variable $Y$:
X Y 0.5 1 2.3 2 2.2 2 1.8 1
Moreover, the covariance matrix between $X$ and $Y$ is also known:
$$ \Sigma_{XY} = \begin{pmatrix} \sigma_{XX} & \sigma_{XY} \\ \sigma_{YX} & \sigma_{YY} \end{pmatrix} = \begin{pmatrix} 0.5 & 1.2 \\ 1.2 & 0.4 \end{pmatrix} $$
Now let's say we have dummy-coded $Y$, so our data looks like this:
X Y1 Y2 0.5 1 0 2.3 0 1 2.2 0 1 1.8 1 0
My question is: how can we calculate the covariance matrix between $X$, $Y_1$ and $Y_2$?
$$ \Sigma_{XY_1Y_2} = \begin{pmatrix} \sigma_{XX} & \sigma_{XY_1} & \sigma_{XY_2} \\ \sigma_{Y_1X} & \sigma_{Y_1Y_1} & \sigma_{Y_1Y_2} \\ \sigma_{Y_2X} & \sigma_{Y_2Y_1} & \sigma_{Y_2Y_2} \end{pmatrix} = \begin{pmatrix} 0.5 & ? & ? \\ ? & ? & ? \\ ? & ? & ? \end{pmatrix} $$
In my actual problem, I can often bootstrap my way into getting a fair estimate for this matrix, but it is very important for me to find the analytical solution as well. For what it's worth, the proportion of each category in the population is also known.