#define nvars 4 ! Fitting bivariate models to raw data - rectangular data files ! ! /data2/boulder/class95/choleskywithmissing.mx ! ! Analysis of two-wave Neuroticism data - Cholesky model - with missing data DA CALC NG=3 MATRICES W LO nvars nvars FR ! additive genetic path X LO nvars nvars ! non-additive genetic path Y LO nvars nvars FR ! shared environmental path Z LO nvars nvars FR ! within-family environmental path M FU 1 8 FR ! vector of means END MATRICES; BEGIN ALGEBRA; A=W*W'; D=X*X'; C=Y*Y'; E=Z*Z'; G=\stnd(A); END ALGEBRA; SP M 101 102 103 104 101 102 103 104 st 0.05 all st 0.0 m(1,1) m(1,2) m(1,3) m(1,4) st 0.5 z(1,1) z(2,2) z(3,3) z(4,4) st 0.2 w(1,1) w(2,2) w(3,3) w(4,4) st 0.2 y(1,1) y(2,2) y(3,3) y(4,4) BO 0.01 10 z(1,1) z(2,2) z(3,3) z(4,4) interval g(1,2) g(1,3) g(1,4) g(2,3) g(2,4) g(3,4) END; MZ female twin pairs DA NI=8 rectangular fi=mzfepqn.rec MATRICES = group 1 END MATRICES; CO (A+D+C+E | A+D+C_ A+D+C | A+D+C+E) / ME M/ OPTIONS RS END FEMALE DZ PAIRS DA NI=8 rectangular fi=dzfepqn.rec 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/ MAt g 0.5 MAt h 0.25 OPTIONS RS END