Example Mx Output
Univariate 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 1
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.5 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=2 NOBSERVATIONS=1500
CMATRIX
2.12
1.67 2.08
LABELS TRAIT1 TRAIT2
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=2 NOBSERVATIONS=2000
CMATRIX
1.98
1.25 2.06
LABELS TRAIT1 TRAIT2
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
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
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
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
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
Matrix not fully standardized
PARAMETER SPECIFICATIONS
GROUP NUMBER: 1
Group1: Defines Matrices
MATRIX A
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX C
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX E
This is a computed FULL matrix of order 1 by 1
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 1 by 1
1
1 1
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 2
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 3
GROUP NUMBER: 2
Group2: MZ twin pairs
MATRIX A
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX C
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX E
This is a computed FULL matrix of order 1 by 1
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 1 by 1
1
1 1
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 2
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 3
GROUP NUMBER: 3
Group3: DZ twin pairs
MATRIX A
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX C
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX E
This is a computed FULL matrix of order 1 by 1
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 1 by 1
1
1 1
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 2
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 3
GROUP NUMBER: 4
Group4: Standardised solution
MATRIX A
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX C
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX D
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX E
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX G
This is a computed FULL matrix of order 1 by 1
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 1 by 1
It has no free parameters specified
MATRIX P
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX S
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX T
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX U
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX V
This is a computed FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX X
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 1
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 2
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 3
MX PARAMETER ESTIMATES
GROUP NUMBER: 1
Group1: Defines Matrices
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
1 0.6890
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
1 0.9335
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
1 0.4287
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 1 by 1
1
1 0.8301
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 0.9662
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 0.6547
GROUP NUMBER: 2
Group2: MZ twin pairs
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
1 0.6890
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
1 0.9335
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
1 0.4287
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 1 by 1
1
1 0.8301
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 0.9662
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 0.6547
OBSERVED COVARIANCE MATRIX
TRAIT1 TRAIT2
TRAIT1 2.1200
TRAIT2 1.6700 2.0800
EXPECTED COVARIANCE MATRIX
TRAIT1 TRAIT2
TRAIT1 2.0512
TRAIT2 1.6225 2.0512
RESIDUAL MATRIX
TRAIT1 TRAIT2
TRAIT1 0.0688
TRAIT2 0.0475 0.0288
Function value of this group: 0.8823
Where the fit function is Maximum Likelihood
GROUP NUMBER: 3
Group3: DZ twin pairs
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
1 0.6890
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
1 0.9335
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
1 0.4287
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 1 by 1
1
1 0.8301
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 0.9662
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 0.6547
OBSERVED COVARIANCE MATRIX
TRAIT1 TRAIT2
TRAIT1 1.9800
TRAIT2 1.2500 2.0600
EXPECTED COVARIANCE MATRIX
TRAIT1 TRAIT2
TRAIT1 2.0512
TRAIT2 1.2780 2.0512
RESIDUAL MATRIX
TRAIT1 TRAIT2
TRAIT1 -0.0712
TRAIT2 -0.0280 0.0088
Function value of this group: 1.6083
Where the fit function is Maximum Likelihood
GROUP NUMBER: 4
Group4: Standardised solution
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
1 0.6890
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
1 0.9335
MATRIX D
This is a computed FULL matrix of order 1 by 1
[=\SQRT(\V2D(\D2V(P)))]
1
1 1.4322
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
1 0.4287
MATRIX G
This is a computed FULL matrix of order 1 by 1
[=\STND(T)]
1
1 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 1 by 1
[=\STND(V)]
1
1 1.0000
MATRIX P
This is a computed FULL matrix of order 1 by 1
[=A+C+E]
1
1 2.0512
MATRIX S
This is a computed FULL matrix of order 1 by 1
[=\STND(U)]
1
1 1.0000
MATRIX T
This is a computed FULL matrix of order 1 by 1
[=D~*A*D~]
1
1 0.3359
MATRIX U
This is a computed FULL matrix of order 1 by 1
[=D~*C*D~]
1
1 0.4551
MATRIX V
This is a computed FULL matrix of order 1 by 1
[=D~*E*D~]
1
1 0.2090
MATRIX X
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 0.8301
MATRIX Y
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 0.9662
MATRIX Z
This is a LOWER TRIANGULAR matrix of order 1 by 1
1
1 0.6547
Your model has 3 estimated parameters and 6 Observed statistics
Chi-squared fit of model >>>>>>> 2.491
Degrees of freedom >>>>>>>>>>>>> 3
Probability >>>>>>>>>>>>>>>>>>>> 0.477
Akaike's Information Criterion > -3.509
RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.003
This problem used 0.1% of my workspace
Task Time elapsed (DD:HH:MM:SS)
Reading script & data 0: 0: 0: 0.04
Execution 0: 0: 0: 0.01
TOTAL 0: 0: 0: 0.05
Total number of warnings issued: 0
______________________________________________________________________________
______________________________________________________________________________