I'm sorry if the title is confusing. I am working with a model with $(X,Y_1,Y_2,...,Y_k)$ such that $X$ is some random variable (a latent factor) and the $Y$'s are generated according to:
$$ Y_i = \mu_i + \lambda_i X + \varepsilon_i$$
with all $\varepsilon$'s and $X$ mutually independent, all $\mu$'s and $\lambda$'s are constant parameters.
My question is: in this case where the $Y$'s are related to $X$ linearly, is it possible to have
$$\mathbb{E}[X|Y_1=y_1, Y_2 = y_2,...,Y_k = y_k] \neq \beta_0 + \beta_1 y_1 + ... +\beta_k y_k$$
(all $\beta$'s are constants and the lowercase y's are the realizations of the big Y's), and what is an example of this case? Any help would be appreciated!
I've played around with copulae on MATLAB but I can't figure this one out... What am I missing here?