Dear Dr. Muthen, I want't to simulate data for the following model: Y = b1*X + b2*Z + b3*X² + b4*XZ + b5*X²Z
Is it possible to use the xwith command? Or would this lead to biased data due to fact that even when X and Z are normally distributed, X² will not. Thus, b5 will be overestimated because of the sensitivity of the LMS approach to nonnormal data. Or isn't this a problem, when I constrain the b-weights?
MODEL MONTECARLO: [x1-x3*3 z1-z3*3 y1-y3*3]; x1-x3*.3; z1-z3*.3; y1-y3*.3; C BY x1*.7 x2*.7 x3*.7; Z BY z1*.7 z2*.7 z3*.7; Y BY y1*.7 y2*.7 y3*.7; C@1; Y@1; Z@1; CZ | C XWITH Z; Y ON C*.4 Z*.4 CZ*.2; Z WITH C*.3;
When I execute the above code, I get an error: *** ERROR in MODEL POPULATION command Variances for the categorical latent variables are not currently defined. Variance given for: C
However if I change all occurrences of "C" in the code, for example to D, the error goes off. What is happening here? Is "C" reserved as a categorical variable or something like that?