G1 CALCULATION GROUP MZ Correlation matrix A factor DA CALC NG=8 MATRICES A Lo 1 1 Free C Lo 1 1 Free I id 1 1 begin algebra; X = A*A' ; Y = C*C' ; end algebra; COMPUTE I|X+Y_ X+Y|I ; Option RS se End G2 CALCULATION GROUP DZ Correlation matrix A factor DA CALC MATRICES A Lo 1 1 =A1 C Lo 1 1 =C1 H fu 1 1 I id 1 1 begin algebra; X = A*A' ; Y = C*C' ; end algebra ; COMPUTE I|H@X+Y_ H@X+Y|I ; Matrix H .5 Option RS End G3 CALCULATION GROUP MZ Correlation matrix AB factor DA CALC MATRICES A Lo 1 1 Free C Lo 1 1 Free I id 1 1 begin algebra; X = A*A' ; Y = C*C' ; end algebra; COMPUTE I|X+Y_ X+Y|I ; Option RS se End G4 CALCULATION GROUP DZ Correlation matrix AB factor DA CALC MATRICES A Lo 1 1 =A3 C Lo 1 1 =C3 H fu 1 1 I id 1 1 begin algebra; X = A*A' ; Y = C*C' ; end algebra ; COMPUTE I|H@X+Y_ H@X+Y|I ; Matrix H .5 Option RS End G5 CALCULATION GROUP MZ Correlation matrix B factor DA CALC MATRICES A Lo 1 1 Free C Lo 1 1 Free I id 1 1 begin algebra; X = A*A' ; Y = C*C' ; end algebra; COMPUTE I|X+Y_ X+Y|I ; Option RS se End G6 CALCULATION GROUP DZ Correlation matrix B factor DA CALC MATRICES A Lo 1 1 =A5 C Lo 1 1 =C5 H fu 1 1 I id 1 1 begin algebra; X = A*A' ; Y = C*C' ; end algebra ; COMPUTE I|H@X+Y_ H@X+Y|I ; Matrix H .5 Option RS End Three disorders model - attempt parameter recovery Data NI=1 NO=1 Matrices A Full 2 2 =%E1 ! corr between A factors C Full 2 2 =%E3 ! corr between AB factors B Full 2 2 =%E5 ! corr between B factors I Iden 1 1 N Full 1 1 ! The scalar 2.0 Q Full 1 1 ! Sample size R Full 10 1! observed data T Full 1 2 ! Threshold for A U Full 1 2 ! Threshold for B V Full 1 2 ! Threshold for AB W Zero 1 2 ! + + These are to control integral type X Zi 1 2 ! + - Z Unit 1 2 ! - - Begin Algebra; D = \muln((A_T_T_Z)) ; E = \muln((A_T_T_X)) ; F = \muln((A_T_T_W)) ; G = \muln((C_U_U_Z)) ; H = \muln((C_U_U_X)) ; J = \muln((C_U_U_W)) ; K = \muln((B_V_V_Z)) ; L = \muln((B_V_V_X)) ; M = \muln((B_V_V_W)) ; Y = \muln(I_T'_I) ; ! below threshold on A O = \muln(I_V'_I) ; ! below threshold on B S = Q @ (D.G.K _ N. D.G.L _ N. E.G.K _ N. (Y.H.O + E.G.L) _ D.G.M _ N. E.G.L _ N. (Y.H.(I-O) + E.G.M) _ F.G.K _ N. ((I-Y).H.O + F.G.L) _ J + F.G.M + N.(I-Y).H.(I-O)) ; end algebra; compute \sum((R-S).(R-S)%S) ; Matrix N 2 Matrix R File=admz.frq Matrix Q File=admz.n Specify T 7 7 Specify U 8 8 Specify V 9 9 Start .6 all Matrix T 1 1 !1.282 1.282 Matrix U 1 1 !2.326 2.326 Matrix V .5 .5 !1 1 !drop @.0 1 Option func=1.e-7 nd=7 !diff=.001 Option user-defined RS End Three disorders model - DZ data Data NI=1 NO=1 Matrices A Full 2 2 =%E2 ! corr between A factors C Full 2 2 =%E4 ! corr between AB factors B Full 2 2 =%E6 ! corr between B factors I Iden 1 1 N Full 1 1 ! The scalar 2.0 Q Full 1 1 ! Sample size R Full 10 1! observed data T Full 1 2 ! Threshold for A U Full 1 2 ! Threshold for B V Full 1 2 ! Threshold for AB W Zero 1 2 ! + + These are to control integral type X Zi 1 2 ! + - Z Unit 1 2 ! - - Begin Algebra; D = \muln((A_T_T_Z)) ; E = \muln((A_T_T_X)) ; F = \muln((A_T_T_W)) ; G = \muln((C_U_U_Z)) ; H = \muln((C_U_U_X)) ; J = \muln((C_U_U_W)) ; K = \muln((B_V_V_Z)) ; L = \muln((B_V_V_X)) ; M = \muln((B_V_V_W)) ; Y = \muln(I_T'_I) ; ! below threshold on A O = \muln(I_V'_I) ; ! below threshold on B S = Q @ (D.G.K _ N. D.G.L _ N. E.G.K _ N. (Y.H.O + E.G.L) _ D.G.M _ N. E.G.L _ N. (Y.H.(I-O) + E.G.M) _ F.G.K _ N. ((I-Y).H.O + F.G.L) _ J + F.G.M + N.(I-Y).H.(I-O)) ; end algebra; compute \sum((R-S).(R-S)%S) ; Matrix N 2 Matrix R File=addz.frq Matrix Q File=addz.n Specify T 7 7 Specify U 8 8 Specify V 9 9 Start .6 all Matrix T 1 1 !1.282 1.282 Matrix U 1 1 !2.326 2.326 Matrix V .5 .5 !1 1 !drop @.0 1 Bound 0 .95 1 2 3 4 5 6 Bound 0 3 7 8 9 !Fix all Option func=1.e-7 nd=7 Option user-defined RS mu !th=-5 nag=10 db=1 End save triace.mxs drop 3 end drop 4 end get triace.mxs drop 4 end