The following MX script lines were read for group 1
#DEFINE NVAR 2
G1: CALCULATION
DATA CALC NGROUPS=3
BEGIN MATRICES;
X LOWER NVAR NVAR FREE ! GENETIC STRUCTURE
Y LOWER NVAR NVAR FREE ! COMMON ENVIRONMENTAL STRUCTURE
Z LOWER NVAR NVAR FREE ! SPECIFIC ENVIRONMENTAL STRUCTURE
G FULL 1 NVAR FREE ! MEANS
H FULL 1 1
END MATRICES;
MATRIX H .5
ST 1 ALL
START 400 G 1 1
START 70 G 1 2
ST 18 X 1 1 Y 1 1 Z 1 1
ST 4 X 2 2 Y 2 2 Z 2 2
BEGIN ALGEBRA;
A= X*X' ;
C= Y*Y' ;
E= Z*Z' ;
P = A%(A+C+E)| C%(A+C+E) | E%(A+C+E) ; !STANDARDIZED ESTIMATES
K=\STND(A) ; !GENETIC CORRELATION
L=\STND(C); !COMMON ENVIRONMENTAL CORRELATION
M=\STND(E); !SPECIFIC ENVIRONMENTAL CORRELATION
END ALGEBRA;
END
The following MX script lines were read for group 2
G2: MZ TWIN PAIRS AND SIB
DA NINPUT_VARS=14
MISSING=-200.00
NOTE: Missing value *string* set to '-200.00'
RECTANGULAR FILE=MRI4.DAT
Rectangular continuous data read initiated
NOTE: Rectangular file contained 322 records with data
that contained a total of 1835 observations
LABELS
FAMNR ZYG WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 WM3 BBGMT3 BBWMT3 WM4 BBGMT4 BBWMT4
SELECT IF ZYG < 3
NOTE: Select if yields 144 data vectors for analysis
NOTE: Vectors contain a total of 809 observations
SELECT BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 BBGMT4 WM4 ;
MATRICES = G 1
NOTE: Selection yields 144 data vectors for analysis
NOTE: Vectors contain a total of 401 observations
MEANS G|G|G|G; !MODEL FOR MEANS
COVARIANCES ! MODEL FOR COVARIANCES
A+C+E | A+C | H@A+C |H@A+C_
A+C | A+C+E | H@A+C |H@A+C_
H@A+C | H@A+C | A+C+E |H@A+C _
H@A+C | H@A+C | H@A+C | A+C+E;
OPTIONS RSIDUAL
END
The following MX script lines were read for group 3
G3: DZ TWIN PAIRS AND SIB
DA NINPUT_VARS=14
MISSING=-200.00
NOTE: Missing value *string* set to '-200.00'
RECTANGULAR FILE=MRI4.DAT
Rectangular continuous data read initiated
NOTE: Rectangular file contained 322 records with data
that contained a total of 1835 observations
LABELS
FAMNR ZYG WM1 BBGMT1 BBWMT1 WM2 BBGMT2 BBWMT2 WM3 BBGMT3 BBWMT3 WM4 BBGMT4 BBWMT4
SELECT IF ZYG > 2
NOTE: Select if yields 178 data vectors for analysis
NOTE: Vectors contain a total of 1026 observations
SELECT BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3 BBGMT4 WM4 ;
MATRICES = G 1
NOTE: Selection yields 178 data vectors for analysis
NOTE: Vectors contain a total of 532 observations
MEANS G|G|G|G; !MODEL FOR MEANS
COVARIANCES ! MODEL FOR COVARIANCES
A+C+E |H@A+C | H@A+C |H@A+C_
H@A+C | A+C+E | H@A+C |H@A+C_
H@A+C | H@A+C | A+C+E |H@A+C _
H@A+C | H@A+C | H@A+C | A+C+E;
OPTIONS RSIDUAL
END
Summary of VL file data for group 2
BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3
Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000
Number 54.0000 106.0000 54.0000 107.0000 12.0000 55.0000
Mean 705.8689 62.3434 697.0530 61.3925 708.9617 61.8987
Variance 3234.4423 157.6016 3779.2227 201.7276 2227.0761 196.3816
Minimum 570.6100 26.1900 523.6600 27.7400 608.9900 30.1200
Maximum 860.5200 90.4800 836.6400 91.7700 805.5700 97.2700
WM4
Code 8.0000
Number 13.0000
Mean 56.6969
Variance 191.3088
Minimum 30.9500
Maximum 75.3600
Summary of VL file data for group 3
BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3
Code 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000
Number 57.0000 146.0000 58.0000 139.0000 19.0000 90.0000
Mean 699.6267 62.3736 682.9521 61.8701 706.4447 66.1716
Variance 2234.7938 181.6186 2733.4405 127.3058 2679.9778 135.6165
Minimum 562.2500 24.9400 568.7000 28.5700 606.5200 32.5700
Maximum 799.9600 94.3700 775.9000 86.6600 798.1000 95.3000
BBGMT4 WM4
Code 7.0000 8.0000
Number 4.0000 19.0000
Mean 694.2000 61.3105
Variance 3169.6562 161.9674
Minimum 630.2600 32.5000
Maximum 779.0800 80.1200
PARAMETER SPECIFICATIONS
GROUP NUMBER: 1
G1: calculation
MATRIX A
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX C
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX E
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX G
This is a FULL matrix of order 1 by 2
1 2
1 10 11
MATRIX H
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX K
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX L
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX M
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX P
This is a computed FULL matrix of order 2 by 6
It has no free parameters specified
MATRIX X
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 1
2 2 3
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 4
2 5 6
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 7
2 8 9
GROUP NUMBER: 2
G2: MZ twin pairs and sib
MATRIX A
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX C
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX E
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX G
This is a FULL matrix of order 1 by 2
1 2
1 10 11
MATRIX H
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX K
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX L
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX M
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX P
This is a computed FULL matrix of order 2 by 6
It has no free parameters specified
MATRIX X
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 1
2 2 3
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 4
2 5 6
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 7
2 8 9
GROUP NUMBER: 3
G3: DZ twin pairs and sib
MATRIX A
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX C
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX E
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX G
This is a FULL matrix of order 1 by 2
1 2
1 10 11
MATRIX H
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX K
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX L
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX M
This is a computed FULL matrix of order 2 by 2
It has no free parameters specified
MATRIX P
This is a computed FULL matrix of order 2 by 6
It has no free parameters specified
MATRIX X
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 1
2 2 3
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 4
2 5 6
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 7
2 8 9
Mx starting optimization; number of parameters = 11
MX PARAMETER ESTIMATES
GROUP NUMBER: 1
G1: calculation
MATRIX A
This is a computed FULL matrix of order 2 by 2
[=X*X']
1 2
1 2030.0534 209.1538
2 209.1538 114.1243
MATRIX C
This is a computed FULL matrix of order 2 by 2
[=Y*Y']
1 2
1 309.7520 5.2243
2 5.2243 0.0881
MATRIX E
This is a computed FULL matrix of order 2 by 2
[=Z*Z']
1 2
1 529.3391 -28.7396
2 -28.7396 56.4224
MATRIX G
This is a FULL matrix of order 1 by 2
1 2
1 696.5500 62.4723
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX K
This is a computed FULL matrix of order 2 by 2
[=\STND(A)]
1 2
1 1.0000 0.4345
2 0.4345 1.0000
MATRIX L
This is a computed FULL matrix of order 2 by 2
[=\STND(C)]
1 2
1 1.0000 1.0000
2 1.0000 1.0000
MATRIX M
This is a computed FULL matrix of order 2 by 2
[=\STND(E)]
1 2
1 1.0000 -0.1663
2 -0.1663 1.0000
MATRIX P
This is a computed FULL matrix of order 2 by 6
[=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)]
1 2 3 4 5
1 7.0755E-01 1.1267E+00 1.0796E-01 2.8142E-02 1.8449E-01
2 1.1267E+00 6.6882E-01 2.8142E-02 5.1638E-04 -1.5481E-01
6
1 -1.5481E-01
2 3.3066E-01
MATRIX X
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 -45.0561
2 -4.6421 9.6216
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 1.7600E+01
2 2.9684E-01 -2.7363E-06
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 23.0074
2 -1.2491 7.4069
GROUP NUMBER: 2
G2: MZ twin pairs and sib
MATRIX A
This is a computed FULL matrix of order 2 by 2
[=X*X']
1 2
1 2030.0534 209.1538
2 209.1538 114.1243
MATRIX C
This is a computed FULL matrix of order 2 by 2
[=Y*Y']
1 2
1 309.7520 5.2243
2 5.2243 0.0881
MATRIX E
This is a computed FULL matrix of order 2 by 2
[=Z*Z']
1 2
1 529.3391 -28.7396
2 -28.7396 56.4224
MATRIX G
This is a FULL matrix of order 1 by 2
1 2
1 696.5500 62.4723
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX K
This is a computed FULL matrix of order 2 by 2
[=\STND(A)]
1 2
1 1.0000 0.4345
2 0.4345 1.0000
MATRIX L
This is a computed FULL matrix of order 2 by 2
[=\STND(C)]
1 2
1 1.0000 1.0000
2 1.0000 1.0000
MATRIX M
This is a computed FULL matrix of order 2 by 2
[=\STND(E)]
1 2
1 1.0000 -0.1663
2 -0.1663 1.0000
MATRIX P
This is a computed FULL matrix of order 2 by 6
[=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)]
1 2 3 4 5
1 7.0755E-01 1.1267E+00 1.0796E-01 2.8142E-02 1.8449E-01
2 1.1267E+00 6.6882E-01 2.8142E-02 5.1638E-04 -1.5481E-01
6
1 -1.5481E-01
2 3.3066E-01
MATRIX X
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 -45.0561
2 -4.6421 9.6216
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 1.7600E+01
2 2.9684E-01 -2.7363E-06
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 23.0074
2 -1.2491 7.4069
Vector of OBSERVED means
BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3
Mean 705.8689 62.3434 697.0530 61.3925 708.9617 61.8987
BBGMT4 WM4
Mean 0.0000 56.6969
Vector of EXPECTED means
BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3
Mean 696.5500 62.4723 696.5500 62.4723 696.5500 62.4723
BBGMT4 WM4
Mean 696.5500 62.4723
EXPECTED COVARIANCE MATRIX
BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3
BBGMT1 2869.1445
WM1 185.6385 170.6348
BBGMT2 2339.8054 214.3780 2869.1445
WM2 214.3780 114.2124 185.6385 170.6348
BBGMT3 1324.7787 109.8012 1324.7787 109.8012 2869.1445
WM3 109.8012 57.1503 109.8012 57.1503 185.6385 170.6348
BBGMT4 1324.7787 109.8012 1324.7787 109.8012 1324.7787 109.8012
WM4 109.8012 57.1503 109.8012 57.1503 109.8012 57.1503
BBGMT4 WM4
BBGMT1
WM1
BBGMT2
WM2
BBGMT3
WM3
BBGMT4 2869.1445
WM4 185.6385 170.6348
Function value of this group: 3422.2464
Where the fit function is -2 * Log-likelihood of raw data
GROUP NUMBER: 3
G3: DZ twin pairs and sib
MATRIX A
This is a computed FULL matrix of order 2 by 2
[=X*X']
1 2
1 2030.0534 209.1538
2 209.1538 114.1243
MATRIX C
This is a computed FULL matrix of order 2 by 2
[=Y*Y']
1 2
1 309.7520 5.2243
2 5.2243 0.0881
MATRIX E
This is a computed FULL matrix of order 2 by 2
[=Z*Z']
1 2
1 529.3391 -28.7396
2 -28.7396 56.4224
MATRIX G
This is a FULL matrix of order 1 by 2
1 2
1 696.5500 62.4723
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX K
This is a computed FULL matrix of order 2 by 2
[=\STND(A)]
1 2
1 1.0000 0.4345
2 0.4345 1.0000
MATRIX L
This is a computed FULL matrix of order 2 by 2
[=\STND(C)]
1 2
1 1.0000 1.0000
2 1.0000 1.0000
MATRIX M
This is a computed FULL matrix of order 2 by 2
[=\STND(E)]
1 2
1 1.0000 -0.1663
2 -0.1663 1.0000
MATRIX P
This is a computed FULL matrix of order 2 by 6
[=A%(A+C+E)|C%(A+C+E)|E%(A+C+E)]
1 2 3 4 5
1 7.0755E-01 1.1267E+00 1.0796E-01 2.8142E-02 1.8449E-01
2 1.1267E+00 6.6882E-01 2.8142E-02 5.1638E-04 -1.5481E-01
6
1 -1.5481E-01
2 3.3066E-01
MATRIX X
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 -45.0561
2 -4.6421 9.6216
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 1.7600E+01
2 2.9684E-01 -2.7363E-06
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 2 by 2
1 2
1 23.0074
2 -1.2491 7.4069
Vector of OBSERVED means
BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3
Mean 699.6267 62.3736 682.9521 61.8701 706.4447 66.1716
BBGMT4 WM4
Mean 694.2000 61.3105
Vector of EXPECTED means
BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3
Mean 696.5500 62.4723 696.5500 62.4723 696.5500 62.4723
BBGMT4 WM4
Mean 696.5500 62.4723
EXPECTED COVARIANCE MATRIX
BBGMT1 WM1 BBGMT2 WM2 BBGMT3 WM3
BBGMT1 2869.1445
WM1 185.6385 170.6348
BBGMT2 1324.7787 109.8012 2869.1445
WM2 109.8012 57.1503 185.6385 170.6348
BBGMT3 1324.7787 109.8012 1324.7787 109.8012 2869.1445
WM3 109.8012 57.1503 109.8012 57.1503 185.6385 170.6348
BBGMT4 1324.7787 109.8012 1324.7787 109.8012 1324.7787 109.8012
WM4 109.8012 57.1503 109.8012 57.1503 109.8012 57.1503
BBGMT4 WM4
BBGMT1
WM1
BBGMT2
WM2
BBGMT3
WM3
BBGMT4 2869.1445
WM4 185.6385 170.6348
Function value of this group: 4550.9224
Where the fit function is -2 * Log-likelihood of raw data
Your model has 11 estimated parameters and 933 Observed statistics
-2 times log-likelihood of data >>> 7973.169
Degrees of freedom >>>>>>>>>>>>>>>> 922
This problem used 5.9% of my workspace
Task Time elapsed (DD:HH:MM:SS)
Reading script & data 0: 0: 0: 0.87
Execution 0: 0: 0: 8.53
TOTAL 0: 0: 0: 9.40
Total number of warnings issued: 0
______________________________________________________________________________