** Mx startup successful ** **MX-VMS 1.27** The following MX script lines were read for group 1 USE MOMENTS OF TRUNCATED NORMAL FUNCTION TO CORRECT FOR ASCERTAINMENT #DEFINE NXY 2 DATA NG=2 NI=NXY NO=489 CMATRIX 0.3375284321 0.1399967791 0.8317113181 MEANS 0.791434 0.374071 MATRICES A LOWR 2 2 FREE ! COVARIANCE OF X,Y M FULL 1 2 FREE ! MEAN OF X,Y T FULL 1 NXY !THRESHOLDS S FULL 1 NXY !SELECTION VECTOR N FULL 1 NXY !# OF ABSCISSAE X FULL 1 4 ! TO PARTITION RESULT OF MOMNOR Y FULL 1 4 ! TO PARTITION RESULT OF MOMNOR END MATRICES MATRIX X 1 1 2 2 ! TO PULL THE 2X2 MATRIX AT TOP OF MOMENT MATRIX MATRIX Y 3 1 3 2 ! TO PULL THE BOTTOM ROW OF THE MOMENT MATRIX MATRIX M 0.5 0.5 ! STARTING VALUES OF MEANS OF VAR1 & VAR2 MATRIX T 0 0 ! THRESHOLDS FOR VAR 1 & 2 MATRIX S 1 0 ! INDICATES SELECTED ABOVE THRESHOLD ON VAR1, NO SEL ON VAR2 MATRIX N 24 24 BEGIN ALGEBRA; B = A*A' ; ! COVARIANCE MATRIX IN POPULATION BEFORE SELECTION C = \MOMNOR(B_M_T_S_N) ; ! COMPUTE THE MOMENTS OF THE TRUNCATED NORMAL ! it will return 3x2 matrix: ! newvar1 newcov12 ! newcov12 newvar2 ! newmean1 newmean2 END ALGEBRA; MEANS \PART(C, Y) / COVARIANCE \PART(C, X) / MATRIX A 1 .1 1 BOUND -2 2 M 1 1 M 1 2 BOUND .1 2 A 1 1 A 2 2 BOUND -1 1 A 2 1 OPTION RS END The following MX script lines were read for group 2 GROUP2: Y IS NOT OBSERVED BECAUSE X<0 #DEFINE NX 1 DATA NI=NX NO=511 CMATRIX 0.3772516322 MEANS -0.797566 MATRICES A LOWR 1 1 ! COVARIANCE OF X,Y M FULL 1 1 ! MEAN OF X,Y T FULL 1 NX !THRESHOLDS S FULL 1 NX !SELECTION VECTOR N FULL 1 NX !# OF ABSCISSAE X FULL 1 4 ! TO PARTITION RESULT OF MOMNOR Y FULL 1 4 ! TO PARTITION RESULT OF MOMNOR END MATRICES MATRIX X 1 1 1 1 MATRIX Y 2 1 2 1 MATRIX M 0 MATRIX T 0 MATRIX S -1 MATRIX N 24 BEGIN ALGEBRA; B = A*A' ; ! COVARIANCE MATRIX IN POPULATION BEFORE SELECTION C = \MOMNOR(B_M_T_S_N) ; ! COMPUTE THE MOMENTS OF THE TRUNCATED NORMAL ! it will return 2x1 matrix: ! newvar1 ! newmean1 END ALGEBRA; MEANS \PART(C, Y) / COVARIANCE \PART(C, X) / SPECIFY A 1 ! EQUALITY CONSTRAINT TO VAR OF X IN GROUP 1 SPECIFY M 4 ! EQUALITY CONSTRAINT TO MEAN OF X IN GRP 1 OPTION NAG=10 DB=10 OPTION MU RS END PARAMETER SPECIFICATIONS GROUP NUMBER: 1 USE MOMENTS OF TRUNCATED NORMAL FUNCTION TO CORRECT FOR ASCERTAINMENT MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1 2 2 3 MATRIX B This is a computed FULL matrix of order 2 by 2 It has no free parameters specified MATRIX C This is a computed FULL matrix of order 3 by 2 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 4 5 MATRIX N This is a FULL matrix of order 1 by 2 It has no free parameters specified MATRIX S This is a FULL matrix of order 1 by 2 It has no free parameters specified MATRIX T This is a FULL matrix of order 1 by 2 It has no free parameters specified MATRIX X This is a FULL matrix of order 1 by 4 It has no free parameters specified MATRIX Y This is a FULL matrix of order 1 by 4 It has no free parameters specified GROUP NUMBER: 2 GROUP2: Y IS NOT OBSERVED BECAUSE X<0 MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX B 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 2 by 1 It has no free parameters specified MATRIX M This is a FULL matrix of order 1 by 1 1 1 4 MATRIX N This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX S This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX T 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 4 It has no free parameters specified MATRIX Y This is a FULL matrix of order 1 by 4 It has no free parameters specified MX PARAMETER ESTIMATES GROUP NUMBER: 1 USE MOMENTS OF TRUNCATED NORMAL FUNCTION TO CORRECT FOR ASCERTAINMENT MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 0.9931 2 0.4162 0.8786 MATRIX B This is a computed FULL matrix of order 2 by 2 1 2 1 0.9863 0.4134 2 0.4134 0.9451 MATRIX C This is a computed FULL matrix of order 3 by 2 1 2 1 0.3508 0.1455 2 0.1455 0.8340 3 0.7798 0.3693 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 -0.0360 0.0273 MATRIX N This is a FULL matrix of order 1 by 2 1 2 1 24.0000 24.0000 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 1.0000 0.0000 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX X This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 2.0000 2.0000 MATRIX Y This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 2.0000 Vector of EXPECTED means 1 2 Mean 0.7798 0.3693 OBSERVED COVARIANCE MATRIX 1 2 1 0.3375 2 0.1400 0.8317 EXPECTED COVARIANCE MATRIX 1 2 1 0.3508 2 0.1455 0.8340 RESIDUAL MATRIX 1 2 1 -0.0132 2 -0.0055 -0.0023 Function value of this group: 0.5444 Where the fit function is Maximum Likelihood GROUP NUMBER: 2 GROUP2: Y IS NOT OBSERVED BECAUSE X<0 MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.9931 MATRIX B This is a computed FULL matrix of order 1 by 1 1 1 0.9863 MATRIX C This is a computed FULL matrix of order 2 by 1 1 1 0.3661 2 -0.8056 MATRIX M This is a FULL matrix of order 1 by 1 1 1 -0.0360 MATRIX N This is a FULL matrix of order 1 by 1 1 1 24.0000 MATRIX S This is a FULL matrix of order 1 by 1 1 1 -1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 1.0000 1.0000 MATRIX Y This is a FULL matrix of order 1 by 4 1 2 3 4 1 2.0000 1.0000 2.0000 1.0000 Vector of EXPECTED means 1 Mean -0.8056 OBSERVED COVARIANCE MATRIX 1 1 0.3773 EXPECTED COVARIANCE MATRIX 1 1 0.3661 RESIDUAL MATRIX 1 1 0.0111 Function value of this group: 0.3193 Where the fit function is Maximum Likelihood Your model has 5 estimated parameters and 7 Observed statistics Chi-squared fit of model >>>>>>> 0.864 Degrees of freedom >>>>>>>>>>>>> 2 Probability >>>>>>>>>>>>>>>>>>>> 0.649 Akaike's Information Criterion > -3.136 This problem used 0.1% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.14 Execution 0: 0: 0: 6.95 TOTAL 0: 0: 0: 7.09 ______________________________________________________________________________ Multiple fit option in effect. The following MX script lines have been read: MATRIX 1 A 1 .5 .866025 MATRIX 1 M 0 0 MATRIX 2 A 1 MATRIX 2 M 0 FIX ALL END PARAMETER SPECIFICATIONS No free parameters specified MX PARAMETER ESTIMATES GROUP NUMBER: 1 USE MOMENTS OF TRUNCATED NORMAL FUNCTION TO CORRECT FOR ASCERTAINMENT MATRIX A This is a LOWER TRIANGULAR matrix of order 2 by 2 1 2 1 1.0000 2 0.5000 0.8660 MATRIX B This is a computed FULL matrix of order 2 by 2 1 2 1 1.0000 0.5000 2 0.5000 1.0000 MATRIX C This is a computed FULL matrix of order 3 by 2 1 2 1 0.3620 0.1810 2 0.1810 0.8405 3 0.7980 0.3990 MATRIX M This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX N This is a FULL matrix of order 1 by 2 1 2 1 24.0000 24.0000 MATRIX S This is a FULL matrix of order 1 by 2 1 2 1 1.0000 0.0000 MATRIX T This is a FULL matrix of order 1 by 2 1 2 1 0.0000 0.0000 MATRIX X This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 2.0000 2.0000 MATRIX Y This is a FULL matrix of order 1 by 4 1 2 3 4 1 3.0000 1.0000 3.0000 2.0000 Vector of EXPECTED means 1 2 Mean 0.7980 0.3990 OBSERVED COVARIANCE MATRIX 1 2 1 0.3375 2 0.1400 0.8317 EXPECTED COVARIANCE MATRIX 1 2 1 0.3620 2 0.1810 0.8405 RESIDUAL MATRIX 1 2 1 -0.0245 2 -0.0410 -0.0088 Function value of this group: 3.3659 Where the fit function is Maximum Likelihood GROUP NUMBER: 2 GROUP2: Y IS NOT OBSERVED BECAUSE X<0 MATRIX A This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1.0000 MATRIX B This is a computed FULL matrix of order 1 by 1 1 1 1.0000 MATRIX C This is a computed FULL matrix of order 2 by 1 1 1 0.3620 2 -0.7980 MATRIX M This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX N This is a FULL matrix of order 1 by 1 1 1 24.0000 MATRIX S This is a FULL matrix of order 1 by 1 1 1 -1.0000 MATRIX T This is a FULL matrix of order 1 by 1 1 1 0.0000 MATRIX X This is a FULL matrix of order 1 by 4 1 2 3 4 1 1.0000 1.0000 1.0000 1.0000 MATRIX Y This is a FULL matrix of order 1 by 4 1 2 3 4 1 2.0000 1.0000 2.0000 1.0000 Vector of EXPECTED means 1 Mean -0.7980 OBSERVED COVARIANCE MATRIX 1 1 0.3773 EXPECTED COVARIANCE MATRIX 1 1 0.3620 RESIDUAL MATRIX 1 1 0.0152 Function value of this group: 0.4392 Where the fit function is Maximum Likelihood Your model has 0 estimated parameters and 7 Observed statistics Chi-squared fit of model >>>>>>> 3.805 Degrees of freedom >>>>>>>>>>>>> 7 Probability >>>>>>>>>>>>>>>>>>>> 0.802 Akaike's Information Criterion > -10.195 This problem used 0.1% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.07 Execution 0: 0: 0: 0.17 TOTAL 0: 0: 0: 0.24 ______________________________________________________________________________ ______________________________________________________________________________