Hello all,

I have a structure model as below

Code:
gsem ($x <- shock $xfirm_control i.year) ///
($y1 <- $x $xfirm_control_rd i.Industry i.Country i.year, family(gaussian,lcensored(0))) ///
($y2 <- $y1 $x $xfirm_control_rd i.Industry i.Country i.year, family(gaussian,lcensored(0))), nocapslatent
after getting the GSEM result, I want to get the marginal effect of each equation by following the code (by using equation 3 as an example)
Code:
eststo:estpost margins, dydx($y1 $x $xfirm_control_rd) predict(ystar(0 .) eq(#3))
however, the result shows

Code:
option ystar() not allowed
r(198);
could you please give me some advice about this? is there any method to get the marginal effect of these models?

thanks in advance.