G1 CALCULATION GROUP control Correlation matrix A factor DA CALC NG=10 MATRICES A Lo 1 1 Fixed C Lo 1 1 Free 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 A 0 Matrix H .5 Option RS End G2 CALCULATION GROUP clinical 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 control Correlation matrix AB factor DA CALC MATRICES A Lo 1 1 Fixed C Lo 1 1 Free 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 A 0 Matrix H .5 Option RS End G4 CALCULATION GROUP clinical 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 control Correlation matrix B factor DA CALC MATRICES A Lo 1 1 Fixed C Lo 1 1 Free 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 A 0 Matrix H .5 Option RS End G6 CALCULATION GROUP clinical 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 CALC 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 T Full 1 2 ! Threshold for A U Full 1 2 ! Threshold for AB V Full 1 2 ! Threshold for B 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 = 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 S; Matrix N 2 Specify T 4 4 Specify U 5 5 Specify V 6 6 Start .6 all Matrix T 2 2 !1.282 1.282 Matrix U 2 2 !2.326 2.326 Matrix V 2 2 !1 1 End Three disorders model - clinical data Data CALC 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 p Full 1 1 free 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 = n. D.G.L _ p. E.G.K _ (n+p). (Y.H.O + E.G.L) _ n.(D.G.M) _ (n+p). E.G.L _ ((n+p)+n). (Y.H.(I-O) + E.G.M) _ p.F.G.K _ (p+(n+p)). ((I-Y).H.O + F.G.L) _ (n+p).(J + F.G.M + N.(I-Y).H.(I-O)) ; end algebra; Compute S; Matrix N 1 Specify T 4 4 Specify U 5 5 Specify V 6 6 Start .6 all bound 0 100 7 Matrix T 2 2 !1.282 1.282 Matrix U 2 2 !2.326 2.326 Matrix V 2 2 !1 1 End G9 Fit model to control group data Data No=1 NI=1 Matrices R Full 10 1 ! observed data Q Full 1 1 ! sample size K Full 10 1 =%E7 ! predicted control group proportions begin algebra; M = \sum(K); U = M_M_M_M_M_M_M_M_M_M; S = Q@(K%U); end algebra; compute \sum((R-S).(R-S)%S) ; Matrix R 244 21 21 15 2 3 1 0 3 3 Matrix Q 313 Option func=1.e-7 nd=7 Option user-defined RS mu !th=-5 nag=10 db=1 end G10 Fit model to clinical group data Data No=1 NI=1 Matrices R Full 9 1 ! observed data Q Full 1 1 ! sample size K Full 9 1 =%E8 ! predicted clinical group proportions begin algebra; M = \sum(K); U = M_M_M_M_M_M_M_M_M; S = Q@(K%U); end algebra; compute \sum((R-S).(R-S)%S) ; Matrix R 67 17 79 5 5 20 3 6 10 Matrix Q 212 Bound 0 .95 1 2 3 Bound 0 3 4 5 6 Option func=1.e-7 nd=7 Option user-defined RS mu !th=-5 nag=10 db=1 Option Th=-30 Start .5 all End