For my textbook, Introduction to Probability by Blitzstein and Hwang, I have the problem where I have the random variables $X$, $Y$, and $Z$ such that $X \sim N(0, 1)$. I am also told that, conditional on $X = x$, $Y$ and $Z$ are independent and identically distribution $N(x, 1)$.
I found the joint PDF of $X$, $Y$, and $Z$ to be
$$P[Y = y, Z = z, X = x] = P[Y = y, Z = z | X = x] P[X = x] = \dfrac{1}{(\sqrt{2 \pi})^3} e^{-1/2 [y^2 + z^2 + 3x^2 - 2x(y + z)]}$$
I now want to find the expected values and variances for each of $X$, $Y$, and $Z$. However, I am unsure of how to do this. I was thinking that, since I have the joint distribution $P[Y = y, Z = z, X = x]$, would I somehow marginalize out each of the other 2 random variables, which would then some how lead to the results I am looking for?
I would greatly appreciate it if someone could please take the time to show and explain how this is done.