I have the following recursive linear model
$X = \varepsilon_1$
$Z = \beta_2 X + \varepsilon_2$
$Y = \beta_3 Z + \varepsilon_3$
By solving the model its easy to see that $E[Y|X, Z]$ is linear too. But what happens in the nonlinear case? Let's say I assume Z and Y to be Probit (or Logit, Poisson, etc)
$X = \varepsilon_1$
$Pr(Z = 1|X) = \Phi(\beta_2 X)$
$Pr(Y = 1|Z) = \Phi(\beta_3 Z)$
Would in this case $Pr(Y = 1|X, Z)$ be Probit too? In preliminary simulations I've run this seems to work. However, I don't know which kind of theoretical result I could apply here. Would be great if you could point me to some literature too. Thanks for your help!