! Multiformity model. ! You can express B&A if you are above second threshold on A or ! above the second threshold on B ! A & B are otherwise independent. ! ! Control correlation matrix A factor Title 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 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 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 Compute predicted control group proportions Data NI=6 ordinal fi=control.dat Labels age1 age2 sex1 sex2 type dummy Definition_variables age1 age2 sex1 sex2 type / Matrices A Full 2 2 =%E1 ! corr between A factors B Full 2 2 =%E3 ! corr between B factors C Full 1 1 Free ! Estimated Lower Threshold for A D Full 1 1 Free ! Age difference in lower threshold for A J Full 1 1 Free ! Sex difference in Lower threshold for A E Full 1 1 ! Estimated Positive Deviation of upper from lower Threshold for A !F Full 1 1 Free ! Age difference in positive deviation of upper from lower threshold for A !L Full 1 1 Free ! Sex difference in positive deviation of upper from lower threshold for A R Full 1 1 Free ! Estimated Lower Threshold for B S Full 1 1 Free ! Age difference in lower threshold for B K Full 1 1 Free ! Sex difference in lower threshold for B T Full 1 1 Free ! Estimated Positive Deviation of upper from lower Threshold for B U Full 1 1 Free ! Age difference in positive deviation of upper from lower threshold for B M Full 1 1 Free ! Sex difference in positive deviation of upper from lower threshold for B V Full 1 1 ! variance of dummy variable O Full 1 1 ! matrix for age definition variable sibling 1 P Full 1 1 ! matrix for age definition variable sibling 2 G Full 1 1 ! matrix for sex definition variable sibling 1 Y Full 1 1 ! matrix for sex definition variable sibling 2 H Full 4 1 ! for extracting right part of I W Full 1 1 ! 0 These are to control integral type X Full 1 1 ! 1 These are to control integral type Z Full 1 1 ! 2 These are to control integral type Begin Algebra; ! C+(O.D)+(G.J) ! Lower threshold for A sibling 1 ! C+(P.D)+(Y.J) ! Lower threshold for A sibling 2 ! E ! positive deviation for threshold A sibling 1 ! C+(O.D)+(G.J)+E ! E ! positive deviation for threshold A sibling 2 ! C+(P.D)+(Y.J)+E ! R+(O.S)+(G.K) ! Lower threshold for B sibling 1 ! R+(P.S)+(Y.K) ! Lower threshold for B sibling 2 ! T+(O.U)+(G.M) ! positive deviation for threshold B sibling 1 ! R+(O.S)+(G.K)+T+(O.U)+(G.M) ! T+(P.U)+(Y.M) ! positive deviation for threhsold B sibling 2 ! R+(P.S)+(Y.K)+T+(P.U)+(Y.M) ! X|X ! lower/lower ! X|Z ! lower/middle ! Z|X ! middle/lower ! X|W ! lower/upper ! W|X ! upper/lower ! Z|Z ! middle/middle ! Z|W ! middle/upper ! W|Z ! upper/middle ! W|W ! upper/upper N = (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|X))); ! LOWER/LOWER STRIPE OF A ! (\muln((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))) ; ! lower/middle stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))) ; ! middle/lower stripe of A ! (\muln((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_X|W))) ; ! lower/upper stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|X))) ; ! upper/lower stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|Z))) ; ! middle/middle stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_Z|W))) ; ! middle/upper stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|Z))) ; ! upper/middle stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_W|W))) ; ! upper/upper stripe of A ! (\muln((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X))) ; ! lower/lower stripe of B ! (\muln((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z))) ; ! lower/middle stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X))) ; ! middle/lower stripe of B ! (\muln((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W))) ; ! lower/upper stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))) ; ! upper/lower stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z))) ; ! middle/middle stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))) ; ! middle/upper stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))) ; ! upper/middle stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_W|W))); ! upper/upper stripe of B I = N.(\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X)))_ ! 1 N.(\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))_ ! 2 N.(\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))_ ! 3 (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X)))_ ! 4 (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))). (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X)))_ ! 5 (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|X))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))))+ N.(\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X)))_ ! 6 (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_X|W))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X)))+ (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))+ (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W))))+ (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))+ (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W))))+ N.(\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W)))_ ! 7 N.(\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))_ ! 8 (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))_ ! 9 (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))). (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))_ ! 10 (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|X))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))))+ N.(\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z)))_ ! 11 (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_X|W))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))))+ (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))))+ N.(\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W)))_ ! 12 (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|Z))). (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X)))_ ! 13 (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|Z))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|Z))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))))+ (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))). (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X)))_ ! 14 (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_Z|W))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X)))+ (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))+ (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|Z))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))+ (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W)))_ ! 15 (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_W|W)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_W|W))) - (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_W|W))). (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_W|W)))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_Z|W))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))))+ (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_X|W))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|Z))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))+ (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|Z))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))))+ (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))). (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z)))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|X))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))); ! 16 End Algebra; Matrix H 1 1 1 1 Matrix V .159154943 Specify C 100 Specify D 101 Specify J 102 !Specify E 103 !Specify F 104 !Specify L 105 Specify R 106 Specify S 107 Specify K 108 Specify T 109 Specify U 110 Specify M 111 Matrix C 1.5 Matrix E 3.0 Matrix R 1.5 Matrix T .5 Matrix W 0 Matrix X 1 Matrix Z 2 Specify O age1; Specify P age2; Specify G sex1; Specify Y sex2; Specify H -5 0 -5 0; Means W; Covariance V; Weight (\part(I,H)) / Bound 0 3 103 109 Option RS End G6 Compute predicted clinical group proportions Data NI=6 ordinal fi=clinical.dat Labels age1 age2 sex1 sex2 type dummy Definition_variables age1 age2 sex1 sex2 type / Matrices A Full 2 2 =%E1 ! corr between A factors B Full 2 2 =%E3 ! corr between B factors Q Full 1 1 free C Full 1 1 =C5 ! Estimated Lower Threshold for A D Full 1 1 =D5 ! Age difference in lower threshold for A J Full 1 1 =J5 ! Sex difference in Lower threshold for A E Full 1 1 =E5 ! Estimated Positive Deviation of upper from lower Threshold for A !F Full 1 1 =F5 ! Age difference in positive deviation of upper from lower threshold for A !L Full 1 1 =L5 ! Sex difference in positive deviation of upper from lower threshold for A R Full 1 1 =R5 ! Estimated Lower Threshold for B S Full 1 1 =S5 ! Age difference in lower threshold for B K Full 1 1 =K5 ! Sex difference in lower threshold for B T Full 1 1 =T5 ! Estimated Positive Deviation of upper from lower Threshold for B U Full 1 1 =U5 ! Age difference in positive deviation of upper from lower threshold for B M Full 1 1 =M5 ! Sex difference in positive deviation of upper from lower threshold for B V Full 1 1 ! variance of dummy variable O Full 1 1 ! matrix for age definition variable sibling 1 P Full 1 1 ! matrix for age definition variable sibling 2 G Full 1 1 ! matrix for sex definition variable sibling 1 Y Full 1 1 ! matrix for sex definition variable sibling 2 H Full 4 1 ! for extracting right part of I W Full 1 1 ! 0 These are to control integral type X Full 1 1 ! 1 These are to control integral type Z Full 1 1 ! 2 These are to control integral type Q Full 1 1 free ! ascertainment parameter Begin Algebra; N = (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|X))); ! LOWER/LOWER STRIPE OF A ! (\muln((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))) ; ! lower/middle stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))) ; ! middle/lower stripe of A ! (\muln((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_X|W))) ; ! lower/upper stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|X))) ; ! upper/lower stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|Z))) ; ! middle/middle stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_Z|W))) ; ! middle/upper stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|Z))) ; ! upper/middle stripe of A ! (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_W|W))) ; ! upper/upper stripe of A ! (\muln((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X))) ; ! lower/lower stripe of B ! (\muln((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z))) ; ! lower/middle stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X))) ; ! middle/lower stripe of B ! (\muln((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W))) ; ! lower/upper stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))) ; ! upper/lower stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z))) ; ! middle/middle stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))) ; ! middle/upper stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))) ; ! upper/middle stripe of B ! (\muln((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_W|W))); ! upper/upper stripe of B I = (x).(N.(\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X))))_ ! 2 (q).((\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X))))_ ! 4 (x+q).((\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|X))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))))+ N.(\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))))_ ! 6 (x).(N.(\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z))))_ ! 8 (q).((\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z))))_ ! 9 (x).((\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))). (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X))))_ ! 10 (x+q).((\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|X))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))))+ N.(\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))))_ ! 11 (x).((\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_X|W))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))))+ (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))))+ N.(\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))))_ ! 12 (q).((\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|Z))). (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X))))_ ! 13 (x+q).((\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|Z))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|Z))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))))+ (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))). (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X))))_ ! 14 (q).((\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_Z|W))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|X)))+ (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))+ (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|Z))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))+ (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W))))_ ! 15 (x+q).((\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_W|W)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_W|W))) - (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_W|W))). (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_W|W)))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_Z|W))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))))+ (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_X|W))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|X)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|Z))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_X|Z)))+ (\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))))+ (\muln((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|Z))). ((\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)_Z|Z)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z))))+ (\MULN((A_C+(O.D)+(G.J)|C+(P.D)+(Y.J)+E_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_X|Z))). (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)_W|Z)))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_W|X))). ((\MULN((B_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_X|W)))+ (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W))))+ (\MULN((A_C+(O.D)+(G.J)+E|C+(P.D)+(Y.J)_C+(O.D)+(G.J)|C+(P.D)+(Y.J)_Z|X))). (\MULN((B_R+(O.S)+(G.K)+T+(O.U)+(G.M)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_R+(O.S)+(G.K)|R+(P.S)+(Y.K)+T+(P.U)+(Y.M)_Z|W)))); ! 16 end algebra; Matrix H 1 1 1 1 Matrix V .159154943 Matrix W 0 Matrix X 1 Matrix Z 2 Matrix Q .5 Specify O age1; Specify P age2; Specify G sex1; Specify Y sex2; Specify H -5 0 -5 0; Means W; Covariance V; Weight (\part((I%( \sum(I)_\sum(I)_\sum(I)_\sum(I)_\sum(I)_\sum(I)_\sum(I)_\sum(I)_\sum(I)_\sum(I)_\sum(I)_\sum(I))),H)) / Option func=1.e-7 Option RS Option Th=-3 !Start .9 all End