G1 Alternate forms model ! ! If you are above threshold on L, you get disorder A with prob p ! If you are above threshold on L, you get disorder B with prob r ! Note: (q=1-p and s=1-r). ! If you are below threshold on L you are immune to A&B ! ! control group correlation matrix DA CALC NG=4 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 SE End G2 CALCULATION GROUP Clinical group 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 Evaluate model control group Data NO=1 NI=1 Matrices A Full 2 2 = %E1 ! corr between A factors I Iden 1 1 N Full 1 1 ! The scalar 2.0 O Full 10 1 ! observed data P Full 1 1 ! probability of getting A if above thresh R Full 1 1 ! probability of getting B if above thresh T Full 1 2 ! Threshold for A W Zero 1 2 ! + + These are to control integral type X Zi 1 2 ! + - Y Full 1 1 ! 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)) ; Q = I - P; S = I - R; K = D + N.E.Q.S + F.Q.S.Q.S _ N.(E.R.Q + F.Q.R.Q.S) _ N.(E.P.S + F.P.S.Q.S) _ N.(E.P.R + F.P.R.Q.S) _ F.Q.R.Q.R _ N.(F.P.Q.R.S) _ N.(F.P.Q.R.R) _ F.P.S.P.S _ N.(F.P.S.P.R) _ F.P.R.P.R ; M = \sum(K); U = M_M_M_M_M_M_M_M_M_M; L = Y@(K%U) ; end algebra; compute \sum((L-O).(L-O)%L) ; !_\sum(K)+\sum(L); Matrix O 244 21 21 15 2 3 1 0 3 3 Matrix Y 313 Matrix N 2 Matrix P .5 Matrix R .5 Specify T 2 2 Specify P 3 Specify R 4 Option RS User-def End G4 Evaluate model, Clinical group Data NI=1 NO=1 Matrices A Full 2 2 = %E1 ! corr between A factors, clinical I Iden 1 1 O Full 9 1 ! observed data P Full 1 1 =P3! probability of getting A if above thresh R Full 1 1 =R3! probability of getting B if above thresh T Full 1 2 =T3! Threshold for A W Zero 1 2 ! + + These are to control integral type X Zi 1 2 ! + - Y Full 1 1 ! Sample size Z Unit 1 2 ! - - b full 1 1 c full 1 1 free Begin Algebra; D = \muln((A_T_T_Z)) ; E = \muln((A_T_T_X)) ; F = \muln((A_T_T_W)) ; Q = I - P; S = I - R; G = B +C; K = b.(E.R.Q + F.Q.R.Q.S) _ c.(E.P.S + F.P.S.Q.S) _ g.(E.P.R + F.P.R.Q.S) _ b.(F.Q.R.Q.R) _ (b+c).(F.P.Q.R.S) _ (b+g).(F.P.Q.R.R) _ c.(F.P.S.P.S) _ (c+g).(F.P.S.P.R) _ g.(F.P.R.P.R) ; M = \sum(K); U = M_M_M_M_M_M_M_M_M; L = Y@(K%U) ; end algebra; compute \sum((L-O).(L-O)%L) ; !_\sum(K)+\sum(L); Matrix b 1 !Matrix c 1 Matrix O 67 17 79 5 5 20 3 6 10 Matrix Y 212 !Start .9 All !Fix all Bound 0 1 1 3 4 Bound 0 3 2 Bound 0 100 5 Option func=.000001 Option RS User-def IT=1000 Start .5 all Option Th=-30 End