! Fitting univariate models to raw data - rectangular data files ! ! univariateraw_nocovariates.mx ! ! Analysis of Social Non-conformity data - univariate genetic model DA CALC NG=3 BEGIN MATRICES; W LO 1 1 FR ! additive genetic path X LO 1 1 ! non-additive genetic path (fixed in this application) Y LO 1 1 FR ! shared environmental path Z LO 1 1 FR ! within-family environmental path M FU 1 2 FR ! vector of means (intercepts) K FU 6 2 ! vector for definition-variable kontrol variables b fu 1 6 FI ! vector of regression coefficients END MATRICES; SP K -2 0 -3 0 -4 0 0 -2 0 -3 0 -4 BEGIN ALGEBRA; A=W*W'; D=X*X'; C=Y*Y'; E=Z*Z'; N=A+D+C+E; I=N~*A; J=N~*C; L=N~*E; END ALGEBRA; SP M 101 101 !SP B !102 !103 !104 !102 !103 !104 st 0.75 m(1,1) st 0.5 z(1,1) st 0.2 w(1,1) st 0.2 y(1,1) st 0.0 b(1,1) st 0.00 b(1,2) st 0.00 b(1,3) BO 0.01 10 z(1,1) END; MZ female twin pairs DA NI=7 rectangular fi=epqlmzf.rec labels obs epqlt1 epqlt2 age agesq agecb zyg definition_variables obs age agesq agecb zyg / BEGIN MATRICES = group 1 ; END MATRICES; CO (A+D+C+E | A+D+C_ A+D+C | A+D+C+E) / ME M+b*k/ OPTIONS RS END FEMALE DZ PAIRS DA NI=7 rectangular fi=epqldzf.rec labels obs epqlt1 epqlt2 age agesq agecb zyg definition_variables obs age agesq agecb zyg/ begin matrices = group 1 ; g di 1 1 h di 1 1 END MATRICES; CO (A+D+C+E | g*A+h*D+C_ g*A+h*D+C | A+D+C+E) / ME M+b*K/ MAt g 0.5 MAt h 0.25 OPTIONS RS END