Example Mx Output
Multivariate ACE Twin Model
** Mx startup successful **
**MX-Linux version 1.47c**
! Mx ACE script for twin data
The following MX script lines were read for group 1
#DEFINE NVAR 3
GROUP1: DEFINES MATRICES
DATA CALC NGROUPS=4
BEGIN MATRICES;
X LOWER NVAR NVAR FREE ! GENETIC STRUCTURE
Y LOWER NVAR NVAR FREE ! SHARED ENVIRONMENT
Z LOWER NVAR NVAR FREE ! NON-SHARED PATH CO-EFFICIENTS
H FULL 1 1
I IDEN 4 4
END MATRICES;
MATRIX H .5
BEGIN ALGEBRA;
A= X*X'; ! GENETIC COVARIANCE MATRIX
C= Y*Y'; ! ENVIRONMENTAL COVARIANCE MATRIX
E= Z*Z'; ! NONSHARED ENVIRONMENTAL COVARIANCE MATRIX
END ALGEBRA;
START 0 X 1 1 TO X NVAR NVAR
START 0 Y 1 1 TO Y NVAR NVAR
START 1.7 Z 1 1 TO Z NVAR NVAR
END
The following MX script lines were read for group 2
GROUP2: MZ TWIN PAIRS
DATA NINPUT_VARS=6 NOBSERVATIONS=1500
CMATRIX
4.24756
-.280441 3.39532
2.01213 .05206 2.95854
3.2865 -.243047 2.01504 4.39321
-.210246 2.38388 .096846 -.202754 3.42946
2.02091 .050704 2.44859 2.02949 .114946 2.86705
LABELS X1 Y1 Z1 X2 Y2 Z2
MATRICES= GROUP 1
COVARIANCES A + C + E | A + C _
A + C | A + C + E /
OPTION RS
END
The following MX script lines were read for group 3
GROUP3: DZ TWIN PAIRS
DATA NINPUT_VARS=6 NOBSERVATIONS=2000
CMATRIX
4.35496
-.354591 3.68418
2.04503 .097915 3.17886
2.12497 .329648 1.41837 4.15677
.241837 1.64101 .430632 -.216949 3.43144
1.34426 .501343 1.70592 2.01647 .139561 3.12658
LABELS X1 Y1 Z1 X2 Y2 Z2
MATRICES= GROUP 1
COVARIANCES A + C + E | H@A + C _
H@A + C | A + C + E /
OPTION RS
END
The following MX script lines were read for group 4
GROUP4: STANDARDISED SOLUTION
CALCULATION
MATRICES = GROUP 1
BEGIN ALGEBRA;
! PHENOTYPIC COVARIANCE MATRIX
P = A + C + E;
! DIAGONAL MATRIX OF PHENOTYPIC STANDARD DEVIATIONS
D = \SQRT(\V2D(\D2V(P)));
! PHENOTYPICALLY STANDARDISED GENETIC COVARIANCE MATRIX
T = D~ * A * D~;
! PHENOTYPICALLY STANDARDISED SHARED ENVIRONMENTAL COVARIANCE MATRIX
U = D~ * C * D~;
! PHENOTYPICALLY STANDARDISED NONSHARED ENVIRONMENTAL COVARIANCE MATRIX
V = D~ * E * D~;
! GENETIC CORRELATION MATRIX
G = \STND(T);
How am I supposed to take the square root of 0.?
Diagonal elements are:
0. 0.
Matrix is
0.0000E+00 0.0000E+00 0.0000E+00
0.0000E+00 0.0000E+00 0.0000E+00
0.0000E+00 0.0000E+00 0.0000E+00
Matrix not fully standardized
! SHARED ENVIRONMENTAL CORRELATION MATRIX
S = \STND(U);
How am I supposed to take the square root of 0.?
Diagonal elements are:
0. 0.
Matrix is
0.0000E+00 0.0000E+00 0.0000E+00
0.0000E+00 0.0000E+00 0.0000E+00
0.0000E+00 0.0000E+00 0.0000E+00
Matrix not fully standardized
! NONSHARED ENVIRONMENTAL CORRELATION MATRIX
N = \STND(V);
How am I supposed to take the square root of 0.?
Diagonal elements are:
0. 0.
Matrix is
0.0000E+00 0.0000E+00 0.0000E+00
0.0000E+00 0.0000E+00 0.0000E+00
0.0000E+00 0.0000E+00 0.0000E+00
Matrix not fully standardized
END ALGEBRA;
END
How am I supposed to take the square root of 0.?
Diagonal elements are:
0. 0.
Matrix is
0.0000E+00 0.0000E+00 0.0000E+00
0.0000E+00 0.0000E+00 0.0000E+00
0.0000E+00 0.0000E+00 0.0000E+00
Matrix not fully standardized
How am I supposed to take the square root of 0.?
Diagonal elements are:
0. 0.
Matrix is
0.0000E+00 0.0000E+00 0.0000E+00
0.0000E+00 0.0000E+00 0.0000E+00
0.0000E+00 0.0000E+00 0.0000E+00
Matrix not fully standardized
PARAMETER SPECIFICATIONS
GROUP NUMBER: 1
Group1: Defines Matrices
MATRIX A
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX C
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX E
This is a computed FULL matrix of order 3 by 3
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 I
This is an IDENTITY matrix of order 4 by 4
MATRIX X
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1
2 2 3
3 4 5 6
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 7
2 8 9
3 10 11 12
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 13
2 14 15
3 16 17 18
GROUP NUMBER: 2
Group2: MZ twin pairs
MATRIX A
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX C
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX E
This is a computed FULL matrix of order 3 by 3
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 I
This is an IDENTITY matrix of order 4 by 4
MATRIX X
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1
2 2 3
3 4 5 6
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 7
2 8 9
3 10 11 12
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 13
2 14 15
3 16 17 18
GROUP NUMBER: 3
Group3: DZ twin pairs
MATRIX A
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX C
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX E
This is a computed FULL matrix of order 3 by 3
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 I
This is an IDENTITY matrix of order 4 by 4
MATRIX X
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1
2 2 3
3 4 5 6
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 7
2 8 9
3 10 11 12
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 13
2 14 15
3 16 17 18
GROUP NUMBER: 4
Group4: Standardised solution
MATRIX A
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX C
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX D
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX E
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX G
This is a computed FULL matrix of order 3 by 3
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 I
This is an IDENTITY matrix of order 4 by 4
MATRIX N
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX P
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX S
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX T
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX U
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX V
This is a computed FULL matrix of order 3 by 3
It has no free parameters specified
MATRIX X
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1
2 2 3
3 4 5 6
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 7
2 8 9
3 10 11 12
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 13
2 14 15
3 16 17 18
*** WARNING! ***
I am not sure I have found a solution that satisfies
Kuhn-Tucker conditions for a minimum.
NAG's IFAIL parameter is 1
We probably have a minimum here, but you might consider trying different
starting values. You can randomize these with TH=n on the OU line, where
n is the number of times you wish to do this.
I STRONGLY recommend BOundaries to be set if you use TH
MX PARAMETER ESTIMATES
GROUP NUMBER: 1
Group1: Defines Matrices
MATRIX A
This is a computed FULL matrix of order 3 by 3
[=X*X']
1 2 3
1 2.1760 -1.0631 1.2854
2 -1.0631 1.6775 -0.6641
3 1.2854 -0.6641 1.8815
MATRIX C
This is a computed FULL matrix of order 3 by 3
[=Y*Y']
1 2 3
1 1.0759 0.8058 0.7510
2 0.8058 0.7914 0.7540
3 0.7510 0.7540 0.7195
MATRIX E
This is a computed FULL matrix of order 3 by 3
[=Z*Z']
1 2 3
1 1.0355 -0.0177 -0.0010
2 -0.0177 1.0384 0.0083
3 -0.0010 0.0083 0.4697
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX I
This is an IDENTITY matrix of order 4 by 4
MATRIX X
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.4751
2 -0.7207 1.0762
3 0.8714 -0.0336 1.0588
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.0373E+00
2 7.7683E-01 4.3352E-01
3 7.2405E-01 4.4189E-01 -9.1639E-07
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.0176
2 -0.0174 1.0189
3 -0.0010 0.0081 0.6853
GROUP NUMBER: 2
Group2: MZ twin pairs
MATRIX A
This is a computed FULL matrix of order 3 by 3
[=X*X']
1 2 3
1 2.1760 -1.0631 1.2854
2 -1.0631 1.6775 -0.6641
3 1.2854 -0.6641 1.8815
MATRIX C
This is a computed FULL matrix of order 3 by 3
[=Y*Y']
1 2 3
1 1.0759 0.8058 0.7510
2 0.8058 0.7914 0.7540
3 0.7510 0.7540 0.7195
MATRIX E
This is a computed FULL matrix of order 3 by 3
[=Z*Z']
1 2 3
1 1.0355 -0.0177 -0.0010
2 -0.0177 1.0384 0.0083
3 -0.0010 0.0083 0.4697
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX I
This is an IDENTITY matrix of order 4 by 4
MATRIX X
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.4751
2 -0.7207 1.0762
3 0.8714 -0.0336 1.0588
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.0373E+00
2 7.7683E-01 4.3352E-01
3 7.2405E-01 4.4189E-01 -9.1639E-07
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.0176
2 -0.0174 1.0189
3 -0.0010 0.0081 0.6853
OBSERVED COVARIANCE MATRIX
X1 Y1 Z1 X2 Y2 Z2
X1 4.2476
Y1 -0.2804 3.3953
Z1 2.0121 0.0521 2.9585
X2 3.2865 -0.2430 2.0150 4.3932
Y2 -0.2102 2.3839 0.0968 -0.2028 3.4295
Z2 2.0209 0.0507 2.4486 2.0295 0.1149 2.8670
EXPECTED COVARIANCE MATRIX
X1 Y1 Z1 X2 Y2 Z2
X1 4.2874
Y1 -0.2751 3.5073
Z1 2.0353 0.0982 3.0708
X2 3.2519 -0.2574 2.0364 4.2874
Y2 -0.2574 2.4690 0.0899 -0.2751 3.5073
Z2 2.0364 0.0899 2.6010 2.0353 0.0982 3.0708
RESIDUAL MATRIX
X1 Y1 Z1 X2 Y2 Z2
X1 -0.0398
Y1 -0.0053 -0.1120
Z1 -0.0232 -0.0462 -0.1122
X2 0.0346 0.0143 -0.0213 0.1058
Y2 0.0471 -0.0851 0.0069 0.0723 -0.0779
Z2 -0.0155 -0.0392 -0.1524 -0.0059 0.0167 -0.2037
Function value of this group: 11.3182
Where the fit function is Maximum Likelihood
GROUP NUMBER: 3
Group3: DZ twin pairs
MATRIX A
This is a computed FULL matrix of order 3 by 3
[=X*X']
1 2 3
1 2.1760 -1.0631 1.2854
2 -1.0631 1.6775 -0.6641
3 1.2854 -0.6641 1.8815
MATRIX C
This is a computed FULL matrix of order 3 by 3
[=Y*Y']
1 2 3
1 1.0759 0.8058 0.7510
2 0.8058 0.7914 0.7540
3 0.7510 0.7540 0.7195
MATRIX E
This is a computed FULL matrix of order 3 by 3
[=Z*Z']
1 2 3
1 1.0355 -0.0177 -0.0010
2 -0.0177 1.0384 0.0083
3 -0.0010 0.0083 0.4697
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX I
This is an IDENTITY matrix of order 4 by 4
MATRIX X
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.4751
2 -0.7207 1.0762
3 0.8714 -0.0336 1.0588
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.0373E+00
2 7.7683E-01 4.3352E-01
3 7.2405E-01 4.4189E-01 -9.1639E-07
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.0176
2 -0.0174 1.0189
3 -0.0010 0.0081 0.6853
OBSERVED COVARIANCE MATRIX
X1 Y1 Z1 X2 Y2 Z2
X1 4.3550
Y1 -0.3546 3.6842
Z1 2.0450 0.0979 3.1789
X2 2.1250 0.3296 1.4184 4.1568
Y2 0.2418 1.6410 0.4306 -0.2169 3.4314
Z2 1.3443 0.5013 1.7059 2.0165 0.1396 3.1266
EXPECTED COVARIANCE MATRIX
X1 Y1 Z1 X2 Y2 Z2
X1 4.2874
Y1 -0.2751 3.5073
Z1 2.0353 0.0982 3.0708
X2 2.1639 0.2742 1.3937 4.2874
Y2 0.2742 1.6302 0.4220 -0.2751 3.5073
Z2 1.3937 0.4220 1.6603 2.0353 0.0982 3.0708
RESIDUAL MATRIX
X1 Y1 Z1 X2 Y2
X1 6.7568E-02
Y1 -7.9488E-02 1.7684E-01
Z1 9.6869E-03 -3.0825E-04 1.0809E-01
X2 -3.8908E-02 5.5443E-02 2.4665E-02 -1.3062E-01
Y2 -3.2368E-02 1.0833E-02 8.6615E-03 5.8154E-02 -7.5898E-02
Z2 -4.9445E-02 7.9372E-02 4.5648E-02 -1.8873E-02 4.1338E-02
Z2
X1
Y1
Z1
X2
Y2
Z2 5.5810E-02
Function value of this group: 11.0339
Where the fit function is Maximum Likelihood
GROUP NUMBER: 4
Group4: Standardised solution
MATRIX A
This is a computed FULL matrix of order 3 by 3
[=X*X']
1 2 3
1 2.1760 -1.0631 1.2854
2 -1.0631 1.6775 -0.6641
3 1.2854 -0.6641 1.8815
MATRIX C
This is a computed FULL matrix of order 3 by 3
[=Y*Y']
1 2 3
1 1.0759 0.8058 0.7510
2 0.8058 0.7914 0.7540
3 0.7510 0.7540 0.7195
MATRIX D
This is a computed FULL matrix of order 3 by 3
[=\SQRT(\V2D(\D2V(P)))]
1 2 3
1 2.0706 0.0000 0.0000
2 0.0000 1.8728 0.0000
3 0.0000 0.0000 1.7524
MATRIX E
This is a computed FULL matrix of order 3 by 3
[=Z*Z']
1 2 3
1 1.0355 -0.0177 -0.0010
2 -0.0177 1.0384 0.0083
3 -0.0010 0.0083 0.4697
MATRIX G
This is a computed FULL matrix of order 3 by 3
[=\STND(T)]
1 2 3
1 1.0000 -0.5564 0.6353
2 -0.5564 1.0000 -0.3738
3 0.6353 -0.3738 1.0000
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX I
This is an IDENTITY matrix of order 4 by 4
MATRIX N
This is a computed FULL matrix of order 3 by 3
[=\STND(V)]
1 2 3
1 1.0000 -0.0171 -0.0015
2 -0.0171 1.0000 0.0119
3 -0.0015 0.0119 1.0000
MATRIX P
This is a computed FULL matrix of order 3 by 3
[=A+C+E]
1 2 3
1 4.2874 -0.2751 2.0353
2 -0.2751 3.5073 0.0982
3 2.0353 0.0982 3.0708
MATRIX S
This is a computed FULL matrix of order 3 by 3
[=\STND(U)]
1 2 3
1 1.0000 0.8732 0.8536
2 0.8732 1.0000 0.9992
3 0.8536 0.9992 1.0000
MATRIX T
This is a computed FULL matrix of order 3 by 3
[=D~*A*D~]
1 2 3
1 0.5075 -0.2742 0.3542
2 -0.2742 0.4783 -0.2024
3 0.3542 -0.2024 0.6127
MATRIX U
This is a computed FULL matrix of order 3 by 3
[=D~*C*D~]
1 2 3
1 0.2509 0.2078 0.2070
2 0.2078 0.2256 0.2298
3 0.2070 0.2298 0.2343
MATRIX V
This is a computed FULL matrix of order 3 by 3
[=D~*E*D~]
1 2 3
1 2.4153E-01 -4.5754E-03 -2.8605E-04
2 -4.5754E-03 2.9606E-01 2.5342E-03
3 -2.8605E-04 2.5342E-03 1.5297E-01
MATRIX X
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.4751
2 -0.7207 1.0762
3 0.8714 -0.0336 1.0588
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.0373E+00
2 7.7683E-01 4.3352E-01
3 7.2405E-01 4.4189E-01 -9.1639E-07
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 3 by 3
1 2 3
1 1.0176
2 -0.0174 1.0189
3 -0.0010 0.0081 0.6853
*** WARNING! ***
Minimization may not be successful. See above
CODE GREEN - it probably was OK
Your model has 18 estimated parameters and 42 Observed statistics
Chi-squared fit of model >>>>>>> 22.352
Degrees of freedom >>>>>>>>>>>>> 24
Probability >>>>>>>>>>>>>>>>>>>> 0.558
Akaike's Information Criterion > -25.648
RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.000
This problem used 0.2% of my workspace
Task Time elapsed (DD:HH:MM:SS)
Reading script & data 0: 0: 0: 0.04
Execution 0: 0: 0: 0.25
TOTAL 0: 0: 0: 0.29
Total number of warnings issued: 2
______________________________________________________________________________
______________________________________________________________________________