! My first Mx Model
! Univariate example testing ADE/ACE model
! Dutch 3 year old Males - Externalizing problems from CBCL
!
! A - additive genetic variance (VA or a^2)
! C - shared environmental variance (VC or c^2)
! E - unique environmental variance (VE or e^2)
! D - dominance genetic variance (VD or d^2)
!
! Phenotypic variance = A + C + E + D
! MZ covariance = A + C + D
! DZ covariance = 0.5 A + C + 0.25 D
!
! Path Coefficients Model
!
! Group 1 calculates variance components for females
! Group 2 fits model to MZ female data
! Group 3 fits model to DZ female data
The following MX script lines were read for group 1
#DEFINE NVAR 1
! Mx replaces every occurrence of nvar with 1
#NGROUP 3
Note: #NGroup set number of groups to 3
TITLE GROUP 1: PARAMETERS OF MODEL
CALCULATION
BEGIN MATRICES;
X FULL NVAR NVAR FREE ! A: ADDITIVE GENETIC PARAMETER
Y FULL NVAR NVAR FREE ! C: SHARED ENVIRONMENTAL PARAMETER
Z FULL NVAR NVAR FREE ! E: UNIQUE ENVIRONMENTAL PARAMETER
W FULL NVAR NVAR ! D: DOMINANCE PARAMETER
H FULL 1 1 ! SCALAR, .5
Q FULL 1 1 ! SCALAR, .25
END MATRICES;
MATRIX H .5
MATRIX Q .25
START .1 ALL
BEGIN ALGEBRA;
A=X*X'; ! A^2: ADDITIVE GENETIC VARIANCE
C=Y*Y'; ! C^2: SHARED ENVIRONMENTAL VARIANCE
E=Z*Z'; ! E^2: UNIQUE ENVIRONMENTAL VARIANCE
D=W*W'; ! D^2: DOMINANCE VARIANCE
V=A+C+E+D; ! TOTAL VARIANCE
P=A|C|E|D; ! PUT PARAMETER ESTIMATES IN ONE MATRIX
! the | operator concatenates matrices with same number of rows
S=P@V~; ! STANDARDIZED PARAMETER ESTIMATES
! the @ kronecker product multiplies every element of the first matrix
! with the second matrix, ~ takes the inverse a matrix
END ALGEBRA;
LABELS ROW X PAREST_A
LABELS ROW Y PAREST_C
LABELS ROW Z PAREST_E
LABELS ROW W PAREST_D
LABELS ROW A A^2
LABELS ROW C C^2
LABELS ROW E E^2
LABELS ROW D D^2
LABELS ROW V VARIANCE
LABELS ROW P ESTIMATE
LABELS COL P A C E D
LABELS ROW S STANDEST
LABELS COL S A^2 C^2 E^2 D^2
END
The following MX script lines were read for group 2
TITLE G2: MALE MZ TWIN PAIRS
Note: Opening #include file 1 extmzm.dat
DATA NINPUTVARS=2 NOBSERVATIONS=1029
LABELS EXT3T1 EXT3T2
CMATRIX FILE=EXTMZM.COV
Note: Closing #include file 1
MATRICES = GROUP 1
COVARIANCES A+C+E+D | A+C+D _
A+C+D | A+C+E+D ;
! the _ operator concatenates matrices with same number of columns
OPTION RSIDUALS ! PRINTS OBSERVED, EXPECTED & RESIDUAL MATRIX
END
The following MX script lines were read for group 3
TITLE G3: MALE DZ TWIN PAIRS
Note: Opening #include file 1 extdzm.dat
DATA NINPUT=2 NOBSERVATION=1024
LABELS EXT3T1 EXT3T2
CMATRIX FILE=EXTDZM.COV
Note: Closing #include file 1
MATRICES= GROUP 1
COVARIANCES A+C+E+D | H@A+C+Q@D _
H@A+C+Q@D | A+C+E+D ;
OPTION RSIDUALS
OPTION NDECIMALS=4
OPTION MULTIPLE ISSAT
END
PARAMETER SPECIFICATIONS
GROUP NUMBER: 1
Title Group 1: parameters of model
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 H
This is a 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 4
It has no free parameters specified
MATRIX Q
This is a 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 4
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 W
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 1
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 2
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 3
GROUP NUMBER: 2
Title G2: male 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 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 H
This is a 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 4
It has no free parameters specified
MATRIX Q
This is a 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 4
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 W
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 1
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 2
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 3
GROUP NUMBER: 3
Title G3: male 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 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 H
This is a 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 4
It has no free parameters specified
MATRIX Q
This is a 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 4
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 W
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 1
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 2
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 3
Mx starting optimization; number of parameters = 3
MX PARAMETER ESTIMATES
GROUP NUMBER: 1
Title Group 1: parameters of model
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
A^2 52.9935
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
C^2 29.9932
MATRIX D
This is a computed FULL matrix of order 1 by 1
[=W*W']
1
D^2 0.0000
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
E^2 22.1802
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX P
This is a computed FULL matrix of order 1 by 4
[=A|C|E|D]
A C E D
ESTIMATE 52.9935 29.9932 22.1802 0.0000
MATRIX Q
This is a FULL matrix of order 1 by 1
1
1 0.2500
MATRIX S
This is a computed FULL matrix of order 1 by 4
[=P@V~]
A^2 C^2 E^2 D^2
STANDEST 0.5039 0.2852 0.2109 0.0000
MATRIX V
This is a computed FULL matrix of order 1 by 1
[=A+C+E+D]
1
VARIANCE 105.1669
MATRIX W
This is a FULL matrix of order 1 by 1
1
PAREST_D 0.0000
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 7.2797
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 5.4766
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 4.7096
GROUP NUMBER: 2
Title G2: male MZ twin pairs
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
A^2 52.9935
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
C^2 29.9932
MATRIX D
This is a computed FULL matrix of order 1 by 1
[=W*W']
1
D^2 0.0000
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
E^2 22.1802
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX P
This is a computed FULL matrix of order 1 by 4
[=A|C|E|D]
A C E D
ESTIMATE 52.9935 29.9932 22.1802 0.0000
MATRIX Q
This is a FULL matrix of order 1 by 1
1
1 0.2500
MATRIX S
This is a computed FULL matrix of order 1 by 4
[=P@V~]
A^2 C^2 E^2 D^2
STANDEST 0.5039 0.2852 0.2109 0.0000
MATRIX V
This is a computed FULL matrix of order 1 by 1
[=A+C+E+D]
1
VARIANCE 105.1669
MATRIX W
This is a FULL matrix of order 1 by 1
1
PAREST_D 0.0000
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 7.2797
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 5.4766
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 4.7096
OBSERVED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 110.4300
EXT3T2 87.3080 108.7900
EXPECTED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 105.1669
EXT3T2 82.9867 105.1669
RESIDUAL MATRIX
EXT3T1 EXT3T2
EXT3T1 5.2631
EXT3T2 4.3213 3.6231
Function value of this group: 1.2547
Where the fit function is Maximum Likelihood
GROUP NUMBER: 3
Title G3: male DZ twin pairs
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
A^2 52.9935
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
C^2 29.9932
MATRIX D
This is a computed FULL matrix of order 1 by 1
[=W*W']
1
D^2 0.0000
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
E^2 22.1802
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX P
This is a computed FULL matrix of order 1 by 4
[=A|C|E|D]
A C E D
ESTIMATE 52.9935 29.9932 22.1802 0.0000
MATRIX Q
This is a FULL matrix of order 1 by 1
1
1 0.2500
MATRIX S
This is a computed FULL matrix of order 1 by 4
[=P@V~]
A^2 C^2 E^2 D^2
STANDEST 0.5039 0.2852 0.2109 0.0000
MATRIX V
This is a computed FULL matrix of order 1 by 1
[=A+C+E+D]
1
VARIANCE 105.1669
MATRIX W
This is a FULL matrix of order 1 by 1
1
PAREST_D 0.0000
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 7.2797
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 5.4766
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 4.7096
OBSERVED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 97.9530
EXT3T2 53.5340 105.2900
EXPECTED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 105.1669
EXT3T2 56.4899 105.1669
RESIDUAL MATRIX
EXT3T1 EXT3T2
EXT3T1 -7.2139
EXT3T2 -2.9559 0.1231
Function value of this group: 2.7727
Where the fit function is Maximum Likelihood
Your model has 3 estimated parameters and 6 Observed statistics
Chi-squared fit of model >>>>>>> 4.027
Degrees of freedom >>>>>>>>>>>>> 3
Probability >>>>>>>>>>>>>>>>>>>> 0.259
Akaike's Information Criterion > -1.973
RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.014
This problem used 0.5% of my workspace
Task Time elapsed (DD:HH:MM:SS)
Reading script & data 0: 0: 0: 0.00
Execution 0: 0: 0: 0.17
TOTAL 0: 0: 0: 0.17
Total number of warnings issued: 0
______________________________________________________________________________
Multiple fit option in effect.
The following MX script lines have been read:
SAVE EXT3MACE.MXS
! Test significance of additive genetic variance
DROP @0 X 1 1 1
END
PARAMETER SPECIFICATIONS
GROUP NUMBER: 1
Title Group 1: parameters of model
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 H
This is a 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 4
It has no free parameters specified
MATRIX Q
This is a 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 4
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 W
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX X
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 2
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 3
GROUP NUMBER: 2
Title G2: male 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 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 H
This is a 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 4
It has no free parameters specified
MATRIX Q
This is a 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 4
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 W
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX X
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 2
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 3
GROUP NUMBER: 3
Title G3: male 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 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 H
This is a 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 4
It has no free parameters specified
MATRIX Q
This is a 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 4
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 W
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX X
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 2
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 3
Mx starting optimization; number of parameters = 2
MX PARAMETER ESTIMATES
GROUP NUMBER: 1
Title Group 1: parameters of model
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
A^2 0.0000
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
C^2 70.4622
MATRIX D
This is a computed FULL matrix of order 1 by 1
[=W*W']
1
D^2 0.0000
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
E^2 35.1633
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX P
This is a computed FULL matrix of order 1 by 4
[=A|C|E|D]
A C E D
ESTIMATE 0.0000 70.4622 35.1633 0.0000
MATRIX Q
This is a FULL matrix of order 1 by 1
1
1 0.2500
MATRIX S
This is a computed FULL matrix of order 1 by 4
[=P@V~]
A^2 C^2 E^2 D^2
STANDEST 0.0000 0.6671 0.3329 0.0000
MATRIX V
This is a computed FULL matrix of order 1 by 1
[=A+C+E+D]
1
VARIANCE 105.6255
MATRIX W
This is a FULL matrix of order 1 by 1
1
PAREST_D 0.0000
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 0.0000
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 8.3942
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 5.9299
GROUP NUMBER: 2
Title G2: male MZ twin pairs
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
A^2 0.0000
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
C^2 70.4622
MATRIX D
This is a computed FULL matrix of order 1 by 1
[=W*W']
1
D^2 0.0000
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
E^2 35.1633
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX P
This is a computed FULL matrix of order 1 by 4
[=A|C|E|D]
A C E D
ESTIMATE 0.0000 70.4622 35.1633 0.0000
MATRIX Q
This is a FULL matrix of order 1 by 1
1
1 0.2500
MATRIX S
This is a computed FULL matrix of order 1 by 4
[=P@V~]
A^2 C^2 E^2 D^2
STANDEST 0.0000 0.6671 0.3329 0.0000
MATRIX V
This is a computed FULL matrix of order 1 by 1
[=A+C+E+D]
1
VARIANCE 105.6255
MATRIX W
This is a FULL matrix of order 1 by 1
1
PAREST_D 0.0000
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 0.0000
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 8.3942
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 5.9299
OBSERVED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 110.4300
EXT3T2 87.3080 108.7900
EXPECTED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 105.6255
EXT3T2 70.4622 105.6255
RESIDUAL MATRIX
EXT3T1 EXT3T2
EXT3T1 4.8045
EXT3T2 16.8458 3.1645
Function value of this group: 98.9044
Where the fit function is Maximum Likelihood
GROUP NUMBER: 3
Title G3: male DZ twin pairs
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
A^2 0.0000
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
C^2 70.4622
MATRIX D
This is a computed FULL matrix of order 1 by 1
[=W*W']
1
D^2 0.0000
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
E^2 35.1633
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX P
This is a computed FULL matrix of order 1 by 4
[=A|C|E|D]
A C E D
ESTIMATE 0.0000 70.4622 35.1633 0.0000
MATRIX Q
This is a FULL matrix of order 1 by 1
1
1 0.2500
MATRIX S
This is a computed FULL matrix of order 1 by 4
[=P@V~]
A^2 C^2 E^2 D^2
STANDEST 0.0000 0.6671 0.3329 0.0000
MATRIX V
This is a computed FULL matrix of order 1 by 1
[=A+C+E+D]
1
VARIANCE 105.6255
MATRIX W
This is a FULL matrix of order 1 by 1
1
PAREST_D 0.0000
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 0.0000
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 8.3942
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 5.9299
OBSERVED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 97.9530
EXT3T2 53.5340 105.2900
EXPECTED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 105.6255
EXT3T2 70.4622 105.6255
RESIDUAL MATRIX
EXT3T1 EXT3T2
EXT3T1 -7.6725
EXT3T2 -16.9282 -0.3355
Function value of this group: 65.4867
Where the fit function is Maximum Likelihood
Your model has 2 estimated parameters and 6 Observed statistics
Chi-squared fit of model >>>>>>> 164.391
Degrees of freedom >>>>>>>>>>>>> 4
Probability >>>>>>>>>>>>>>>>>>>> 0.000
Akaike's Information Criterion > 156.391
RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.197
Saturated model fit* >>>>>>>>>>> 4.027
Saturated model df* >>>>>>>>>>> 3
Difference Chi-squared >>>>>>>> 160.364
Difference d.f. >>>>>>>>>>>>>>> 1
Probability >>>>>>>>>>>>>>>>>>>> 0.000
Akaike's Information Criterion > 158.364
* Saturated model statistic computed earlier in this job
This problem used 0.5% of my workspace
Task Time elapsed (DD:HH:MM:SS)
Reading script & data 0: 0: 0: 0.00
Execution 0: 0: 0: 0.00
TOTAL 0: 0: 0: 0.00
Total number of warnings issued: 0
______________________________________________________________________________
! Test significance of shared environmental variance
Multiple fit option in effect.
The following MX script lines have been read:
GET EXT3MACE.MXS
Note: 1351 binary records read
DROP Y 1 1 1
END
PARAMETER SPECIFICATIONS
GROUP NUMBER: 1
Title Group 1: parameters of model
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 H
This is a 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 4
It has no free parameters specified
MATRIX Q
This is a 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 4
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 W
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 1
MATRIX Y
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 3
GROUP NUMBER: 2
Title G2: male 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 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 H
This is a 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 4
It has no free parameters specified
MATRIX Q
This is a 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 4
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 W
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 1
MATRIX Y
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 3
GROUP NUMBER: 3
Title G3: male 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 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 H
This is a 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 4
It has no free parameters specified
MATRIX Q
This is a 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 4
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 W
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 1
MATRIX Y
This is a FULL matrix of order 1 by 1
It has no free parameters specified
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 3
Mx starting optimization; number of parameters = 2
MX PARAMETER ESTIMATES
GROUP NUMBER: 1
Title Group 1: parameters of model
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
A^2 80.8553
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
C^2 0.0000
MATRIX D
This is a computed FULL matrix of order 1 by 1
[=W*W']
1
D^2 0.0000
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
E^2 21.1793
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX P
This is a computed FULL matrix of order 1 by 4
[=A|C|E|D]
A C E D
ESTIMATE 80.8553 0.0000 21.1793 0.0000
MATRIX Q
This is a FULL matrix of order 1 by 1
1
1 0.2500
MATRIX S
This is a computed FULL matrix of order 1 by 4
[=P@V~]
A^2 C^2 E^2 D^2
STANDEST 0.7924 0.0000 0.2076 0.0000
MATRIX V
This is a computed FULL matrix of order 1 by 1
[=A+C+E+D]
1
VARIANCE 102.0346
MATRIX W
This is a FULL matrix of order 1 by 1
1
PAREST_D 0.0000
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 8.9920
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 0.0000
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 4.6021
GROUP NUMBER: 2
Title G2: male MZ twin pairs
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
A^2 80.8553
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
C^2 0.0000
MATRIX D
This is a computed FULL matrix of order 1 by 1
[=W*W']
1
D^2 0.0000
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
E^2 21.1793
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX P
This is a computed FULL matrix of order 1 by 4
[=A|C|E|D]
A C E D
ESTIMATE 80.8553 0.0000 21.1793 0.0000
MATRIX Q
This is a FULL matrix of order 1 by 1
1
1 0.2500
MATRIX S
This is a computed FULL matrix of order 1 by 4
[=P@V~]
A^2 C^2 E^2 D^2
STANDEST 0.7924 0.0000 0.2076 0.0000
MATRIX V
This is a computed FULL matrix of order 1 by 1
[=A+C+E+D]
1
VARIANCE 102.0346
MATRIX W
This is a FULL matrix of order 1 by 1
1
PAREST_D 0.0000
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 8.9920
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 0.0000
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 4.6021
OBSERVED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 110.4300
EXT3T2 87.3080 108.7900
EXPECTED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 102.0346
EXT3T2 80.8553 102.0346
RESIDUAL MATRIX
EXT3T1 EXT3T2
EXT3T1 8.3954
EXT3T2 6.4527 6.7554
Function value of this group: 4.4305
Where the fit function is Maximum Likelihood
GROUP NUMBER: 3
Title G3: male DZ twin pairs
MATRIX A
This is a computed FULL matrix of order 1 by 1
[=X*X']
1
A^2 80.8553
MATRIX C
This is a computed FULL matrix of order 1 by 1
[=Y*Y']
1
C^2 0.0000
MATRIX D
This is a computed FULL matrix of order 1 by 1
[=W*W']
1
D^2 0.0000
MATRIX E
This is a computed FULL matrix of order 1 by 1
[=Z*Z']
1
E^2 21.1793
MATRIX H
This is a FULL matrix of order 1 by 1
1
1 0.5000
MATRIX P
This is a computed FULL matrix of order 1 by 4
[=A|C|E|D]
A C E D
ESTIMATE 80.8553 0.0000 21.1793 0.0000
MATRIX Q
This is a FULL matrix of order 1 by 1
1
1 0.2500
MATRIX S
This is a computed FULL matrix of order 1 by 4
[=P@V~]
A^2 C^2 E^2 D^2
STANDEST 0.7924 0.0000 0.2076 0.0000
MATRIX V
This is a computed FULL matrix of order 1 by 1
[=A+C+E+D]
1
VARIANCE 102.0346
MATRIX W
This is a FULL matrix of order 1 by 1
1
PAREST_D 0.0000
MATRIX X
This is a FULL matrix of order 1 by 1
1
PAREST_A 8.9920
MATRIX Y
This is a FULL matrix of order 1 by 1
1
PAREST_C 0.0000
MATRIX Z
This is a FULL matrix of order 1 by 1
1
PAREST_E 4.6021
OBSERVED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 97.9530
EXT3T2 53.5340 105.2900
EXPECTED COVARIANCE MATRIX
EXT3T1 EXT3T2
EXT3T1 102.0346
EXT3T2 40.4277 102.0346
RESIDUAL MATRIX
EXT3T1 EXT3T2
EXT3T1 -4.0816
EXT3T2 13.1063 3.2554
Function value of this group: 34.6379
Where the fit function is Maximum Likelihood
Your model has 2 estimated parameters and 6 Observed statistics
Chi-squared fit of model >>>>>>> 39.068
Degrees of freedom >>>>>>>>>>>>> 4
Probability >>>>>>>>>>>>>>>>>>>> 0.000
Akaike's Information Criterion > 31.068
RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.080
Saturated model fit* >>>>>>>>>>> 4.027
Saturated model df* >>>>>>>>>>> 3
Difference Chi-squared >>>>>>>> 35.041
Difference d.f. >>>>>>>>>>>>>>> 1
Probability >>>>>>>>>>>>>>>>>>>> 0.000
Akaike's Information Criterion > 33.041
* Saturated model statistic computed earlier in this job
This problem used 0.5% of my workspace
Task Time elapsed (DD:HH:MM:SS)
Reading script & data 0: 0: 0: 0.00
Execution 0: 0: 0: 0.05
TOTAL 0: 0: 0: 0.05
Total number of warnings issued: 0
______________________________________________________________________________