Find marginal distribution of $Y$ where $Y\mid X$ is $N(a_1+a_2X,\sigma_1^2)$ and $X$ is $N(\mu,\sigma^2)$?

Let a random variable $X$ be normal $N(\mu,\sigma^2)$ and let the conditional distribution of $Y$ given $X$ be normal $N(a_1+a_2X,\sigma_1^2)$.

a)Find the joint probability density function of $X$ and $Y$.

b)Find the marginal distribution of $Y$ and the correlation coefficient of $X$ and $Y$.

For (a), I just multiplied the conditional density of $Y$ given $X$ and density of $X$; and I think it's ok. For (b), I tried to write their joint density in the form of bivariate normal but couldn't do that. On the other hand, we know that if the random variables $X$ and $Y$ are bivariate normal then, the conditional distribution of $X$ given $Y$ is normal with mean $E[X\mid Y]$ and variance $(1-\rho^2)\sigma_X^2$. But is that true that if $X$ is normal and $Y$ given $X$ is normal, then they are bivariate normal? So that I can do (b) easily, or is there another way to solve this question?

Thanks!

$\newcommand{\var}{\operatorname{var}}\newcommand{\cov}{\operatorname{cov}}$For part (a) you're ok except that I would do some routine simplications to make it clear that the two variables are in symmetrical roles.

You have $Y\mid X\sim N(a_1+a_2X,\sigma_1^2)$ and $X\sim N(\mu,\sigma^2)$.

From $Y\mid X\sim N(a_1+a_2 X,\sigma_1^2)$ we get $Y-(a_1+a_2 X)\mid X\sim \underbrace{\quad N(0,\sigma_1^2)\quad}_{\text{No $X$'' here.}}$.

If the conditional distribution of a random variable given $X$ does not depend on $X$, then we can conclude two things:

• That random variable (in this case $Y-(a_1+x_2 X)$) is independent of $X$; and
• That conditional distribution is also the marginal distribution (in this case the marginal distribution of $Y-(a_1+a_2 X)$). Consequently we have $\var(Y-(a_1+a_2 X))=\sigma_1^2$.

So $$\overbrace{\var(Y) = \var(Y-(a_1 + a_2 X)) + \var(a_1+a_2 X)}^{\text{by independence}} = \sigma_1^2 + a_2^2\var(X) = \sigma_1^2+a_2^2\sigma^2.$$ And $$\operatorname{E}(Y) = \operatorname{E}(\operatorname{E}(Y\mid X)) = \operatorname{E}(a_1+a_2 X) = a_1+a_2 \mu.$$ That gives us the marginal distribution of $Y$.

Next $$0 = \cov(Y-(a_1+a_2 X),X) = \cov(Y,X) - a_2\cov(X,X),$$ so $\cov(Y,X)=a_2\sigma^2$. Then use $\operatorname{corr}(X,Y)=\dfrac{\cov(X,Y)}{\sqrt{\var(X)\var(Y)}}$.

As for bivariate normality, consider linear combinations $cX+dY$. Show that that is a linear combination of the two independent random variables mentioned above.