Covariance of random variable with product of random variables I have the following covariance term:
$Cov(x,yz)$
with $x$, $y$, and $z$ being random variables with a mean and variance.
I found a paper by Bohrnstedt and Goldberger from 1969, On The Exact Covariance Of Products Of Random Variables, who show that
$Cov(x,yz) = E[y]Cov(x,z) + E[z]Cov(x,y) + E[(x-E[x])(y-E[y])(z-E[z])]$
Now, I know that $y$ and $z$ are independent ($x$ instead is dependent on $y$ and $z$). $z$ is i.i.d. standard normal. Is there a way to simplify the last equation in this case?
 A: Sometimes you can work these sorts of things out--and gain some insight--by considering the simplest possible non-trivial cases.  That would mean that $X,Y,$ and $Z$ would have only two possible values.  Let's make the calculations easy by supposing that their means are all zero and they have unit variances.  This easily implies all three variables have equal chances of being $\pm 1$: they are Rademacher variables.
There aren't too many such situations, as we can count by observing


*

*The trivariate distribution of $(X,Y,Z)$ is determined by eight probabilities associated with the eight possible non-negative values $(\pm1,\pm1,\pm1).$

*Those eight values sum to unity (a linear constraint).

*Setting three means to zero adds three more linear constraints.


That still leaves $8-3-1=4$ parameters.  For simplicity, let's suppose $Y$ and $Z$ are independent (one additional restriction) and that $X$ has the same correlation coefficient with each of $Y$ and $Z.$  There's still a one-parameter family of mathematical solutions.  There's no assurance any of them give probability distributions, because the solutions might always have negative components.
But let's push ahead anyway.  If we let $8\rho$ denote the covariance of $X$ with $YZ$ and $4\tau$ the common covariances of $X$ with $Y$ and $Z,$ the mathematical solutions for the eight probabilities are given in this table:
$$\begin{array}{cccc}
X & Y & Z & \text{Probability} \\
-1 & -1 & -1 & 1/8 - \rho + \tau \\
-1 & -1 & 1 & 1/8+ \rho \\
-1 & 1 & -1 & 1/8+ \rho \\
-1 & 1 & 1 & 1/8-\rho-\tau \\
1 & -1 & -1 & 1/8 + \rho - \tau \\
1 & -1 & 1 & 1/8- \rho \\
1 & 1 & -1 & 1/8- \rho \\
1 & 1 & 1 & 1/8+\rho+\tau \\
\end{array} $$
You can easily compute that all the linear constraints hold, no matter what the numbers $\rho$ and $\tau$ might be.  It's also obvious that when $\pm\rho \pm \tau \le 1/8,$ all these probabilities are non-negative.  This gives a range of possible values of $8\rho=\operatorname{Cov}(X, YZ)$ from $-1 + |4\tau|$ to $1 - |4\tau|.$
Just as a check, note that with $\tau=0,$ the two extreme examples $8\rho=\pm 1$ give these two distributions:


*

*$8\rho=-1:$ $\Pr(-1,-1,-1)=\Pr(-1,1,1)=\Pr(1,-1,1)=\Pr(1,1,-1)=1/4.$

*$8\rho=1:$ $\Pr(-1,-1,1)=\Pr(-1,1,-1)=\Pr(1,-1,-1)=\Pr(1,1,1)=1/4.$
In both cases you can easily check that when you select any two of the three variables, all four pairs of their values $(\pm1,\pm1)$ are equally likely, whence in both cases these variables are pairwise independent.  Nevertheless, in the first case $YZ$ always has the opposite sign of $X$ and in the second case it always has the same sign.
This example shows that unless you add some more conditions, you might not be able to say anything at all about $\operatorname{Cov}(X,YZ)$ in terms of the expectations and covariances of the individual variables!
