G1 CALCULATION GROUP Control Correlation matrix A factor DA CALC NG=6 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 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 H .5 Matrix A 0 Option RS End G4 CALCULATION GROUP Clinical Correlation matrix B 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 Epiphenomena model - fit model to Control A/D data Data NI=1 NO=1 Matrices A Full 2 2 =%E1 ! corr between A factors B Full 2 2 =%E3 ! corr between B factors I Iden 1 1 N Full 1 1 ! The scalar 2.0 O Full 10 1! observed data P Full 1 1 free ! probability of being comorbid given A R Full 1 1 ! probability of being comorbid given B T Full 1 2 ! Threshold for A U Full 1 2 ! Threshold for B W Zero 1 2 ! + + These are to control integral type X Zi 1 2 ! + - Y Full 1 1 ! for sample size 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((B_U_U_Z)) ; H = \muln((B_U_U_X)) ; J = \muln((B_U_U_W)) ; Q = I - P ; S = I - R ; K = D.G _ N. D.H.S _ N. E.Q.G _ N. (E.(P.G + H) + D.H.R) _ D.J.S.S _ N. E.Q.H.S _ N. (E.S.(P.H + J) + D.J.R.S) _ F.Q.Q.G _ N. (F.(P.Q.G + Q.H) + E.H.R.Q) _ F.(P.P.G + N.P.H + J) + N.E.H.P.R + N.R.E.J + R.R.D.J ; M = \sum(K); V = M_M_M_M_M_M_M_M_M_M; L = Y@(K%V) ; end algebra; compute \sum((L-O).(L-O)%L) ; !_\sum(K)+\sum(L); Specify T 6 6 Specify U 7 7 Matrix T 2 2 Matrix U 2 2 Matrix O 244 21 21 15 2 3 1 0 3 3 Matrix Y 313 Matrix N 2 Matrix P .3 Matrix R 0 Option func=1.e-7 nd=7 !diff=.001 Option user-defined RS End Epiphenomena model - fit model to Clinical A/D data Data NI=1 NO=1 Matrices A Full 2 2 =%E1 ! corr between A factors B Full 2 2 =%E3 ! corr between B factors I Iden 1 1 O Full 9 1! observed data P Full 1 1 =P5 ! probability of being comorbid if you are A R Full 1 1 =R5 ! probability of being comorbid if you are B T Full 1 2 =T5! Threshold for A U Full 1 2 =U5! Threshold for B W Zero 1 2 ! + + These are to control integral type X Zi 1 2 ! + - Y Full 1 1 ! for sample size Z Unit 1 2 ! - - c full 1 1 v full 1 1 free Begin Algebra; D = \muln((A_T_T_Z)) ; E = \muln((A_T_T_X)) ; F = \muln((A_T_T_W)) ; G = \muln((B_U_U_Z)) ; H = \muln((B_U_U_X)) ; J = \muln((B_U_U_W)) ; Q = I - P ; S = I - R ; K = c. D.H.S _ v. E.Q.G _ (c+v). (E.(P.G + H) + D.H.R) _ c.(D.J.S.S) _ (c+v). E.Q.H.S _ ((c+v)+c). (E.S.(P.H + J) + D.J.R.S) _ v.(F.Q.Q.G) _ (v+(c+v)). (F.(P.Q.G + Q.H) + E.H.R.Q) _ (c+v).(F.(P.P.G + (c+c).P.H + J) + (c+c).E.H.P.R + (c+c).R.E.J + R.R.D.J) ; M = \sum(K); N = M_M_M_M_M_M_M_M_M; L = Y@(K%N) ; end algebra; compute \sum((L-O).(L-O)%L) ; !_\sum(K)+\sum(L); Matrix C 1 Matrix O 67 17 79 5 5 20 3 6 10 Matrix P .4 Matrix R 0 Matrix Y 212 Bound .01 .99 3 Bound 0 3 6 7 bound 0 100 8 Option func=1.e-7 nd=7 !diff=.001 Option user-defined RS mu CI 90 option df=18 Option Th=-30 Start .5 all End !Save biepiace.mxs !Drop 6 !End !Get biepiace.mxs !Drop 5 !End