** Mx startup successful **
**MX-PC 1.50** Job started on 03/06/01 at 13:53:19
!Correlations for 3 DZ IBD and MZ groups
!'Partitioned Twin Analysis'
The following MX script lines were read for group 1
#DEFINE NVAR 1
#NGROUPS 4
Note: #NGroup set number of groups to 4
G1: DZ IBD2 TWINS
DATA NINPUT=13
RECTANGULAR FILE=MARKER.DAT
NOTE: Rectangular file contained 519 records with data
LABELS FAMID INDID1 INDID2 PIBD0 PIBD1 PIBD2 ZYGOSITY PLT11 PLT12 PLT13 PLT21 PLT22 PLT23
SELECT IF ZYGOSITY >2;
NOTE: Select if yields 232 data vectors for analysis
NOTE: Vectors contain a total of 2662 observations
SELECT IF PIBD2 >.80;
NOTE: Select if yields 43 data vectors for analysis
NOTE: Vectors contain a total of 500 observations
SELECT PLT11 PLT21;
MATRICES = GROUP 1
NOTE: Selection yields 38 data vectors for analysis
NOTE: Vectors contain a total of 75 observations
S DIAG 2 2 FREE
R STND 2 2 FREE
M FULL 2 1 FREE
END MATRICES;
MATRIX S 1 1
MATRIX R .5
MATRIX M 3 3
MEANS M;
COVARIANCE S*R*S';
END
The following MX script lines were read for group 2
G2: DZ IBD1 TWINS
DATA NINPUT=13
RECTANGULAR FILE=MARKER.DAT
NOTE: Rectangular file contained 519 records with data
LABELS FAMID INDID1 INDID2 PIBD0 PIBD1 PIBD2 ZYGOSITY PLT11 PLT12 PLT13 PLT21 PLT22 PLT23
SELECT IF ZYGOSITY >2;
NOTE: Select if yields 232 data vectors for analysis
NOTE: Vectors contain a total of 2662 observations
SELECT IF PIBD1 >.80;
NOTE: Select if yields 67 data vectors for analysis
NOTE: Vectors contain a total of 771 observations
SELECT PLT11 PLT21;
MATRICES
NOTE: Selection yields 58 data vectors for analysis
NOTE: Vectors contain a total of 114 observations
S DIAG 2 2 FREE
R STND 2 2 FREE
M FULL 2 1 FREE
END MATRICES;
MATRIX S 1 1
MATRIX R .5
MATRIX M 3 3
MEANS M;
COVARIANCE S*R*S';
END
The following MX script lines were read for group 3
G3: DZ IBD0 TWINS
DATA NINPUT=13
RECTANGULAR FILE=MARKER.DAT
NOTE: Rectangular file contained 519 records with data
LABELS FAMID INDID1 INDID2 PIBD0 PIBD1 PIBD2 ZYGOSITY PLT11 PLT12 PLT13 PLT21 PLT22 PLT23
SELECT IF ZYGOSITY >2;
NOTE: Select if yields 232 data vectors for analysis
NOTE: Vectors contain a total of 2662 observations
SELECT IF PIBD0 >.80;
NOTE: Select if yields 37 data vectors for analysis
NOTE: Vectors contain a total of 419 observations
SELECT PLT11 PLT21;
MATRICES
NOTE: Selection yields 32 data vectors for analysis
NOTE: Vectors contain a total of 64 observations
S DIAG 2 2 FREE
R STND 2 2 FREE
M FULL 2 1 FREE
END MATRICES;
MATRIX S 1 1
MATRIX R .5
MATRIX M 3 3
MEANS M;
COVARIANCE S*R*S';
END
The following MX script lines were read for group 4
G4: MZ TWINS
DATA NINPUT=13
RECTANGULAR FILE=MARKER.DAT
NOTE: Rectangular file contained 519 records with data
LABELS FAMID INDID1 INDID2 PIBD0 PIBD1 PIBD2 ZYGOSITY PLT11 PLT12 PLT13 PLT21 PLT22 PLT23
SELECT IF ZYGOSITY <3;
NOTE: Select if yields 287 data vectors for analysis
NOTE: Vectors contain a total of 2962 observations
SELECT PLT11 PLT21;
MATRICES
NOTE: Selection yields 166 data vectors for analysis
NOTE: Vectors contain a total of 328 observations
S DIAG 2 2 FREE
R STND 2 2 FREE
M FULL 2 1 FREE
END MATRICES;
MATRIX S 1 1
MATRIX R .5
MATRIX M 3 3
MEANS M;
COVARIANCE S*R*S';
OPTION MULTIPLE ISSAT
END
Summary of VL file data for group 1
PLT11 PLT21
Code 1.0000 2.0000
Number 38.0000 37.0000
Mean 2.7950 2.7276
Variance 0.2711 0.2431
Minimum 1.9200 1.7200
Maximum 4.2600 4.0200
Summary of VL file data for group 2
PLT11 PLT21
Code 1.0000 2.0000
Number 58.0000 56.0000
Mean 2.9200 2.8020
Variance 0.2255 0.2195
Minimum 1.9800 1.7700
Maximum 3.9400 3.9100
Summary of VL file data for group 3
PLT11 PLT21
Code 1.0000 2.0000
Number 32.0000 32.0000
Mean 2.8994 2.7906
Variance 0.2188 0.3680
Minimum 2.0000 1.6600
Maximum 4.2900 4.5900
Summary of VL file data for group 4
PLT11 PLT21
Code 1.0000 2.0000
Number 165.0000 163.0000
Mean 2.8489 2.8618
Variance 0.2925 0.3018
Minimum 1.8000 1.8000
Maximum 4.4800 4.5800
PARAMETER SPECIFICATIONS
GROUP NUMBER: 1
G1: DZ IBD2 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 4
2 5
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 3 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 1
2 0 2
GROUP NUMBER: 2
G2: DZ IBD1 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 9
2 10
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 8 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 6
2 0 7
GROUP NUMBER: 3
G3: DZ IBD0 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 14
2 15
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 13 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 11
2 0 12
GROUP NUMBER: 4
G4: MZ twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 19
2 20
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 18 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 16
2 0 17
MX PARAMETER ESTIMATES
GROUP NUMBER: 1
G1: DZ IBD2 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.7950
2 2.7285
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.6811 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.5206
2 0.0000 0.4900
GROUP NUMBER: 2
G2: DZ IBD1 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.9200
2 2.8113
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.3424 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.4748
2 0.0000 0.4705
GROUP NUMBER: 3
G3: DZ IBD0 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.8994
2 2.7906
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.1797 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.4677
2 0.0000 0.6066
GROUP NUMBER: 4
G4: MZ twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.8495
2 2.8677
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.8508 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.5397
2 0.0000 0.5504
Your model has 20 estimated parameters and 581 Observed statistics
-2 times log-likelihood of data >>> 657.202
Degrees of freedom >>>>>>>>>>>>>>>> 561
This problem used 26.8% of my workspace
Task Time elapsed (DD:HH:MM:SS)
Reading script & data 0: 0: 0: 2.58
Execution 0: 0: 0: 1.97
TOTAL 0: 0: 0: 4.55
Total number of warnings issued: 0
______________________________________________________________________________
Multiple fit option in effect.
The following MX script lines have been read:
SAVE COR.MXS
! Test for linkage
! Set 3 DZ IBD correlations equal
EQUATE R 1 2 1 R 2 2 1 R 3 2 1
END
Summary of VL file data for group 1
PLT11 PLT21
Code 1.0000 2.0000
Number 38.0000 37.0000
Mean 2.7950 2.7276
Variance 0.2711 0.2431
Minimum 1.9200 1.7200
Maximum 4.2600 4.0200
Summary of VL file data for group 2
PLT11 PLT21
Code 1.0000 2.0000
Number 58.0000 56.0000
Mean 2.9200 2.8020
Variance 0.2255 0.2195
Minimum 1.9800 1.7700
Maximum 3.9400 3.9100
Summary of VL file data for group 3
PLT11 PLT21
Code 1.0000 2.0000
Number 32.0000 32.0000
Mean 2.8994 2.7906
Variance 0.2188 0.3680
Minimum 2.0000 1.6600
Maximum 4.2900 4.5900
Summary of VL file data for group 4
PLT11 PLT21
Code 1.0000 2.0000
Number 165.0000 163.0000
Mean 2.8489 2.8618
Variance 0.2925 0.3018
Minimum 1.8000 1.8000
Maximum 4.4800 4.5800
PARAMETER SPECIFICATIONS
GROUP NUMBER: 1
G1: DZ IBD2 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 4
2 5
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 3 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 1
2 0 2
GROUP NUMBER: 2
G2: DZ IBD1 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 9
2 10
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 3 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 6
2 0 7
GROUP NUMBER: 3
G3: DZ IBD0 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 14
2 15
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 3 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 11
2 0 12
GROUP NUMBER: 4
G4: MZ twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 19
2 20
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 18 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 16
2 0 17
MX PARAMETER ESTIMATES
GROUP NUMBER: 1
G1: DZ IBD2 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.7950
2 2.7282
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.4249 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.4851
2 0.0000 0.4576
GROUP NUMBER: 2
G2: DZ IBD1 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.9200
2 2.8136
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.4249 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.4840
2 0.0000 0.4809
GROUP NUMBER: 3
G3: DZ IBD0 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.8994
2 2.7906
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.4249 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.4966
2 0.0000 0.6441
GROUP NUMBER: 4
G4: MZ twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.8495
2 2.8677
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.8508 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.5397
2 0.0000 0.5504
Your model has 18 estimated parameters and 581 Observed statistics
-2 times log-likelihood of data >>> 665.288
Degrees of freedom >>>>>>>>>>>>>>>> 563
Saturated model fit* >>>>>>>>>>> 657.202
Saturated model df* >>>>>>>>>>> 561
Difference Chi-squared >>>>>>>> 8.086
Difference d.f. >>>>>>>>>>>>>>> 2
Probability >>>>>>>>>>>>>>>>>>>> 0.018
Akaike's Information Criterion > 4.086
* Saturated model statistic computed earlier in this job
This problem used 26.8% of my workspace
Task Time elapsed (DD:HH:MM:SS)
Reading script & data 0: 0: 0: 0.33
Execution 0: 0: 0: 0.99
TOTAL 0: 0: 0: 1.32
Total number of warnings issued: 0
______________________________________________________________________________
Multiple fit option in effect.
The following MX script lines have been read:
GET COR.MXS
Note: 53224 binary records read
! Test for residual polygenic variance
! Set DZ IBD2 correlation equal to MZ correlation
EQUATE R 1 2 1 R 4 2 1
END
Summary of VL file data for group 1
PLT11 PLT21
Code 1.0000 2.0000
Number 38.0000 37.0000
Mean 2.7950 2.7276
Variance 0.2711 0.2431
Minimum 1.9200 1.7200
Maximum 4.2600 4.0200
Summary of VL file data for group 2
PLT11 PLT21
Code 1.0000 2.0000
Number 58.0000 56.0000
Mean 2.9200 2.8020
Variance 0.2255 0.2195
Minimum 1.9800 1.7700
Maximum 3.9400 3.9100
Summary of VL file data for group 3
PLT11 PLT21
Code 1.0000 2.0000
Number 32.0000 32.0000
Mean 2.8994 2.7906
Variance 0.2188 0.3680
Minimum 2.0000 1.6600
Maximum 4.2900 4.5900
Summary of VL file data for group 4
PLT11 PLT21
Code 1.0000 2.0000
Number 165.0000 163.0000
Mean 2.8489 2.8618
Variance 0.2925 0.3018
Minimum 1.8000 1.8000
Maximum 4.4800 4.5800
PARAMETER SPECIFICATIONS
GROUP NUMBER: 1
G1: DZ IBD2 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 4
2 5
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 3 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 1
2 0 2
GROUP NUMBER: 2
G2: DZ IBD1 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 9
2 10
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 8 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 6
2 0 7
GROUP NUMBER: 3
G3: DZ IBD0 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 14
2 15
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 13 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 11
2 0 12
GROUP NUMBER: 4
G4: MZ twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 19
2 20
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 0
2 3 0
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 16
2 0 17
MX PARAMETER ESTIMATES
GROUP NUMBER: 1
G1: DZ IBD2 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.7950
2 2.7287
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.8277 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.6123
2 0.0000 0.5756
GROUP NUMBER: 2
G2: DZ IBD1 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.9200
2 2.8113
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.3424 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.4748
2 0.0000 0.4705
GROUP NUMBER: 3
G3: DZ IBD0 twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.8994
2 2.7906
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.1797 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.4677
2 0.0000 0.6066
GROUP NUMBER: 4
G4: MZ twins
MATRIX M
This is a FULL matrix of order 2 by 1
1
1 2.8495
2 2.8676
MATRIX R
This is a STANDARDISED matrix of order 2 by 2
1 2
1 1.0000
2 0.8277 1.0000
MATRIX S
This is a DIAGONAL matrix of order 2 by 2
1 2
1 0.5234
2 0.0000 0.5335
Your model has 19 estimated parameters and 581 Observed statistics
-2 times log-likelihood of data >>> 662.730
Degrees of freedom >>>>>>>>>>>>>>>> 562
Saturated model fit* >>>>>>>>>>> 657.202
Saturated model df* >>>>>>>>>>> 561
Difference Chi-squared >>>>>>>> 5.528
Difference d.f. >>>>>>>>>>>>>>> 1
Probability >>>>>>>>>>>>>>>>>>>> 0.019
Akaike's Information Criterion > 3.528
* Saturated model statistic computed earlier in this job
This problem used 26.8% of my workspace
Task Time elapsed (DD:HH:MM:SS)
Reading script & data 0: 0: 0: 0.06
Execution 0: 0: 0: 0.99
TOTAL 0: 0: 0: 1.05
Total number of warnings issued: 0
______________________________________________________________________________