**MX-PC 1.68e** Job started on 02/21/14 at 13:51:20 ! F.V. Rijsdijk ! Polychoric Cor Script for CTables, of categorical variables (twin1/twin2) ! Note: the threshold model is different then when using RAW ORD: ! TWIN 2 ! ! 0 1 2 ! T 0 cell1 | cell2 | cell3 ! W ___________________________ row threshold 1 ! I 1 cell4 | cell5 | cell6 ! N ____________________________ row threshold 2 ! 1 2 cell7 | cell8 | cell9 ! ^ ^ ! col thr 1 col thr2 ! THRESHOLD MODEL : ! (1,1) (1,2) = row thresholds = TH of twin1 =T ! (2,1) (2,2) = col thresholds = TH of twin2 =U ! So, need to do T_U #ngroups 2 The following MX script lines were read for group 1 #NGROUPS 2 Note: #NGroup set number of groups to 2 G1: Correlation matrix for MZ pairs Data NInputv=2 CTable 3 3 385 23 6 28 37 4 3 3 12 Begin Matrices; R STAN 2 2 FREE ! MZ cor T FUL 1 2 FREE ! TH of twin1 U FUL 1 2 FREE ! TH of twin2 End matrices; COvariance R/ !Model for the MZ correlation Thresholds T_U / !Model for the thresholds MA T 1.8 2.4 !starting values for the Twin 1 thresholds MA U 1.8 2.4 !starting values for the Twin 2 thresholds MA R .6 !starting value Cor Interval R 2 1 options RS end The following MX script lines were read for group 2 G2: Correlation matrix DZ Data NInputv=2 CTable 3 3 334 37 19 56 32 1 12 2 10 Begin Matrices; R STAN 2 2 FREE ! MZ cor T FUL 1 2 FREE ! TH of twin1 U FUL 1 2 FREE ! TH of twin2 End matrices; COvariance R/ !Model for the MZ correlation Thresholds T_U / !Model for the thresholds MA T 1.8 2.4 !starting values for the Twin 1 thresholds MA U 1.8 2.4 !starting values for the Twin 2 thresholds MA R .6 !starting value Cor Interval R 2 1 options RS end PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: Correlation matrix for MZ pairs MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 1 0 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 2 3 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 4 5 GROUP NUMBER: 2 G2: Correlation matrix DZ MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 0 2 6 0 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 7 8 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 9 10 Mx starting optimization; number of parameters = 10 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: Correlation matrix for MZ pairs MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.8147 1.0000 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 0.9429 1.7850 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 0.9611 1.6894 Expected ROW thresholds 1 2 t 0.9429 1.7850 Expected COLUMN thresholds 1 2 t 0.9611 1.6894 Observed frequencies 1 2 3 1 385.0000 23.0000 6.0000 2 28.0000 37.0000 4.0000 3 3.0000 3.0000 12.0000 Expected frequencies 1 2 3 1 384.6495 27.0920 2.6509 2 30.1982 27.8129 9.9950 3 1.8645 6.5525 10.1846 Residual frequencies 1 2 3 1 0.3505 -4.0920 3.3491 2 -2.1982 9.1871 -5.9950 3 1.1355 -3.5525 1.8154 EXPECTED COVARIANCE MATRIX 1 2 1 1.0000 0.8147 2 0.8147 1.0000 Function value of this group: 14.6370 Where the fit function is ML fit function for counts GROUP NUMBER: 2 G2: Correlation matrix DZ MATRIX R This is a STANDARDISED matrix of order 2 by 2 1 2 1 1.0000 2 0.4426 1.0000 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 0.7544 1.6649 MATRIX U This is a FULL matrix of order 1 by 2 1 2 1 0.8406 1.5543 Expected ROW thresholds 1 2 t 0.7544 1.6649 Expected COLUMN thresholds 1 2 t 0.8406 1.5543 Observed frequencies 1 2 3 1 334.0000 37.0000 19.0000 2 56.0000 32.0000 1.0000 3 12.0000 2.0000 10.0000 Expected frequencies 1 2 3 1 333.3188 42.8612 13.4860 2 57.7075 20.6451 10.8564 3 11.2312 7.0266 5.8672 Residual frequencies 1 2 3 1 0.6812 -5.8612 5.5140 2 -1.7075 11.3549 -9.8564 3 0.7688 -5.0266 4.1328 EXPECTED COVARIANCE MATRIX 1 2 1 1.0000 0.4426 2 0.4426 1.0000 Function value of this group: 30.6500 Where the fit function is ML fit function for counts Your model has 10 estimated parameters and 16 Observed statistics Chi-squared fit of model >>>>>>> 45.287 Degrees of freedom >>>>>>>>>>>>> 6 Probability >>>>>>>>>>>>>>>>>>>> 0.000 Akaike's Information Criterion > 33.287 RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.112 1 Confidence intervals requested in group 1 Matrix Element Int. Estimate Lower Upper Lfail Ufail R 1 2 1 95.0 0.8147 0.7344 0.8746 0 0 0 0 1 Confidence intervals requested in group 2 Matrix Element Int. Estimate Lower Upper Lfail Ufail R 2 2 1 95.0 0.4426 0.3012 0.5682 0 0 0 0 This problem used 0.2% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.01 Execution 0: 0: 0: 0.01 TOTAL 0: 0: 0: 0.02 Total number of warnings issued: 0 ______________________________________________________________________________