**MX-RS6000 version 1.15** The following MX control lines have been read: ! E5. Independent Pathway Model ! Australian Atopy Data (asthma,hayfever,dust allergy,eczema) of Female Twins G1: GENETIC STRUCTURE DATA CALC NGROUPS=5 MATRICES A FULL 4 1 FREE ! GENETIC COMMON FACTOR G DIAG 4 4 FREE ! GENETIC SPECIFICS COMPUTE A*A' + G*G' / END The following MX control lines have been read: G2: DOMINANCE STRUCTURE DATA CALC MATRICES D FULL 4 1 FREE ! DOMINANCE COMMON FACTOR COMPUTE D*D' / END The following MX control lines have been read: G3: SPECIFIC ENVIRONMENTAL STRUCTURE DATA CALC MATRICES E FULL 4 1 FREE ! SPECIFIC ENVIRONMENTAL COMMON FACTOR F DIAG 4 4 FREE ! SPECIFIC ENVIRONMENTAL SPECIFICS COMPUTE E*E' + F*F'/ END The following MX control lines have been read: G4: FEMALE MZ TWIN PAIRS DATA NINPUT_VARS=8 NOBSERVATIONS=1232 LABELS ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ECZEMA2 PMATRIX FILE=AHDEMZF.COR ACOV FILE=AHDEMZF.ACV MATRICES A SYMM 4 4 = %E1 ! EXPECTED MATRIX OF GROUP 1 D SYMM 4 4 = %E2 ! EXPECTED MATRIX OF GROUP 2 E SYMM 4 4 = %E3 ! .. COVARIANCES A+D+E | A+D _ A+D | A+D+E / OPTION RS END The following MX control lines have been read: G5: FEMALE DZ TWIN PAIRS DATA NINPUT_VARS=8 NOBSERVATIONS=751 LABELS ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ECZEMA2 PMATRIX FILE=AHDEDZF.COR ACOV FILE=AHDEDZF.ACV MATRICES A SYMM 4 4 = %E1 D SYMM 4 4 = %E2 E SYMM 4 4 = %E3 H FULL 1 1 Q FULL 1 1 COVARIANCES A+D+E | H@A+Q@D _ H@A+Q@D | A+D+E / MATRIX H .5 MATRIX Q .25 START .3 ALL OPTIONS DFREEDOM=-12 NDECIMALS=4 ITERATIONS=200 OPTION RS DEFAULT END PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: GENETIC STRUCTURE MATRIX A This is a FULL matrix of order 4 by 1 1 1 1 2 2 3 3 4 4 MATRIX G This is a DIAGONAL matrix of order 4 by 4 1 2 3 4 1 5 2 0 6 3 0 0 7 4 0 0 0 8 GROUP NUMBER: 2 G2: DOMINANCE STRUCTURE MATRIX D This is a FULL matrix of order 4 by 1 1 1 9 2 10 3 11 4 12 GROUP NUMBER: 3 G3: SPECIFIC ENVIRONMENTAL STRUCTURE MATRIX E This is a FULL matrix of order 4 by 1 1 1 13 2 14 3 15 4 16 MATRIX F This is a DIAGONAL matrix of order 4 by 4 1 2 3 4 1 17 2 0 18 3 0 0 19 4 0 0 0 20 GROUP NUMBER: 4 G4: FEMALE MZ TWIN PAIRS MATRIX A This is a constrained FULL matrix of order 4 by 4 It has no free parameters specified MATRIX D This is a constrained FULL matrix of order 4 by 4 It has no free parameters specified MATRIX E This is a constrained FULL matrix of order 4 by 4 It has no free parameters specified GROUP NUMBER: 5 G5: FEMALE DZ TWIN PAIRS MATRIX A This is a constrained FULL matrix of order 4 by 4 It has no free parameters specified MATRIX D This is a constrained FULL matrix of order 4 by 4 It has no free parameters specified MATRIX E This is a constrained FULL matrix of order 4 by 4 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Q This is a FULL matrix of order 1 by 1 It has no free parameters specified WARNING - I am not sure I have found a solution that satisfies Kuhn-Tucker conditions for a minimum. NAG's IFAIL parameter is 4 The number of iterations has exceeded maximum. Try increasing it with IT=n on the OU line. Current value is 60 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: GENETIC STRUCTURE MATRIX A This is a FULL matrix of order 4 by 1 1 1 .4300 2 .7717 3 .5155 4 .2217 MATRIX G This is a DIAGONAL matrix of order 4 by 4 1 2 3 4 1 .4422 2 .0000 .0245 3 .0000 .0000 .2952 4 .0000 .0000 .0000 .7122 GROUP NUMBER: 2 G2: DOMINANCE STRUCTURE MATRIX D This is a FULL matrix of order 4 by 1 1 1 .4663 2 .0958 3 .4319 4 .2603 GROUP NUMBER: 3 G3: SPECIFIC ENVIRONMENTAL STRUCTURE MATRIX E This is a FULL matrix of order 4 by 1 1 1 .3206 2 .4957 3 .6579 4 .0908 MATRIX F This is a DIAGONAL matrix of order 4 by 4 1 2 3 4 1 .5471 2 .0000 .3862 3 .0000 .0000 .1672 4 .0000 .0000 .0000 .6064 GROUP NUMBER: 4 G4: FEMALE MZ TWIN PAIRS MATRIX A This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 .3805 .3319 .2217 .0953 2 .3319 .5961 .3979 .1711 3 .2217 .3979 .3529 .1143 4 .0953 .1711 .1143 .5563 MATRIX D This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 .2175 .0447 .2014 .1214 2 .0447 .0092 .0414 .0249 3 .2014 .0414 .1865 .1124 4 .1214 .0249 .1124 .0677 MATRIX E This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 .4021 .1589 .2109 .0291 2 .1589 .3948 .3261 .0450 3 .2109 .3261 .4608 .0597 4 .0291 .0450 .0597 .3760 OBSERVED COVARIANCE MATRIX ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ASTHMA1 1.0000 HAYFVR1 .5556 1.0000 DUSTAL1 .5731 .7581 1.0000 ECZEMA1 .2728 .2636 .3089 1.0000 ASTHMA2 .5918 .3656 .3980 .2316 1.0000 HAYFVR2 .4107 .5931 .4512 .1452 .5485 1.0000 DUSTAL2 .4336 .4209 .5180 .1916 .6402 .7700 1.0000 ECZEMA2 .0874 .1963 .1933 .5885 .1454 .1218 .2183 ECZEMA2 ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ECZEMA2 1.0000 EXPECTED COVARIANCE MATRIX ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ASTHMA1 1.0000 HAYFVR1 .5354 1.0001 DUSTAL1 .6340 .7654 1.0002 ECZEMA1 .2458 .2410 .2864 1.0001 ASTHMA2 .5979 .3765 .4231 .2167 1.0000 HAYFVR2 .3765 .6053 .4392 .1960 .5354 1.0001 DUSTAL2 .4231 .4392 .5394 .2267 .6340 .7654 1.0002 ECZEMA2 .2167 .1960 .2267 .6241 .2458 .2410 .2864 ECZEMA2 ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ECZEMA2 1.0001 RESIDUAL MATRIX ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ASTHMA1 .0000 HAYFVR1 .0202 -.0001 DUSTAL1 -.0609 -.0072 -.0002 ECZEMA1 .0270 .0226 .0225 -.0001 ASTHMA2 -.0061 -.0109 -.0251 .0149 .0000 HAYFVR2 .0342 -.0123 .0119 -.0509 .0131 -.0001 DUSTAL2 .0105 -.0184 -.0214 -.0351 .0062 .0047 -.0002 ECZEMA2 -.1293 .0003 -.0334 -.0355 -.1004 -.1192 -.0681 ECZEMA2 ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ECZEMA2 -.0001 This group has a function value of 15.8603 Where the fit function is Asymptotic Weighted Least Squares Part of value due to diagonal of S is .0001 GROUP NUMBER: 5 G5: FEMALE DZ TWIN PAIRS MATRIX A This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 .3805 .3319 .2217 .0953 2 .3319 .5961 .3979 .1711 3 .2217 .3979 .3529 .1143 4 .0953 .1711 .1143 .5563 MATRIX D This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 .2175 .0447 .2014 .1214 2 .0447 .0092 .0414 .0249 3 .2014 .0414 .1865 .1124 4 .1214 .0249 .1124 .0677 MATRIX E This is a constrained FULL matrix of order 4 by 4 1 2 3 4 1 .4021 .1589 .2109 .0291 2 .1589 .3948 .3261 .0450 3 .2109 .3261 .4608 .0597 4 .0291 .0450 .0597 .3760 MATRIX H This is a FULL matrix of order 1 by 1 1 1 .5000 MATRIX Q This is a FULL matrix of order 1 by 1 1 1 .2500 OBSERVED COVARIANCE MATRIX ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ASTHMA1 1.0000 HAYFVR1 .5240 1.0000 DUSTAL1 .5882 .7490 1.0000 ECZEMA1 .2911 .3142 .2790 1.0000 ASTHMA2 .2617 .1703 .0410 .1388 1.0000 HAYFVR2 .1286 .3176 .2616 .0932 .3952 1.0000 DUSTAL2 .0788 .1706 .2144 .0189 .6835 .7228 1.0000 ECZEMA2 .2168 .1145 .0875 .3129 .2540 .2176 .2757 ECZEMA2 ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ECZEMA2 1.0000 EXPECTED COVARIANCE MATRIX ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ASTHMA1 1.0000 HAYFVR1 .5354 1.0001 DUSTAL1 .6340 .7654 1.0002 ECZEMA1 .2458 .2410 .2864 1.0001 ASTHMA2 .2446 .1771 .1612 .0780 1.0000 HAYFVR2 .1771 .3004 .2093 .0918 .5354 1.0001 DUSTAL2 .1612 .2093 .2231 .0852 .6340 .7654 1.0002 ECZEMA2 .0780 .0918 .0852 .2951 .2458 .2410 .2864 ECZEMA2 ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ECZEMA2 1.0001 RESIDUAL MATRIX ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ASTHMA1 .0000 HAYFVR1 -.0115 -.0001 DUSTAL1 -.0458 -.0164 -.0002 ECZEMA1 .0453 .0732 -.0074 -.0001 ASTHMA2 .0171 -.0068 -.1202 .0608 .0000 HAYFVR2 -.0485 .0172 .0523 .0014 -.1402 -.0001 DUSTAL2 -.0824 -.0386 -.0087 -.0664 .0495 -.0425 -.0002 ECZEMA2 .1388 .0227 .0022 .0178 .0082 -.0234 -.0107 ECZEMA2 ASTHMA1 HAYFVR1 DUSTAL1 ECZEMA1 ASTHMA2 HAYFVR2 DUSTAL2 ECZEMA2 -.0001 This group has a function value of 22.6288 Where the fit function is Asymptotic Weighted Least Squares Part of value due to diagonal of S is .0001 *WARNING* Minimization may not be successful. See above CODE BLUE - not enough iterations Your model has 20 estimated parameters and 60 Observed statistics Note that Observed statistics includes your DF= adjustment of -12 Chi-squared fit of model = 38.489 Degrees of freedom = 40 Probability = .538 Akaike's Information Criterion = -41.511 This problem used .3% of my workspace and took 4.83 seconds of CPU time ______________________________________________________________________________ ______________________________________________________________________________