** Mx startup successful ** **MX-Linux version 1.52b** ! This is example script for k=3 trait test The following MX script lines were read for group 1 #DEFINE TOTAL 6 #NGROUPS 3 Note: #NGroup set number of groups to 3 G1: PARAMETER ESTIMATES ! PARAMETER SPECIFICATION CALCULATION BEGIN MATRICES; A FULL 1 1 FREE ! A IS FREE PARAMETER WHITH DIMENSION 1 BY 1 B FULL 1 1 FREE C FULL 1 1 FREE X FULL TOTAL TOTAL FIXED ! X IS 6 BY 6 MATRIX WITH FIXED VALUES ASSIGNED IN THE FOLLOWING LINES END MATRICES; ! (THIS IS SUBMATRIX OF INVERTED FISHER INFORMATION MATRIXTHAT CORRESPOND TO SIX PARAMETER OF INTERESTS VALUE 7.710474640438282 X 1 1 ! ASSIGN THE VALUE TO EACH ELEMENT OF X FROM THE CALCULATED SUBMATRIX OF INVERTED FISHER INFORMATION VALUE 7.709260663212016 X 2 2 VALUE 7.711691774397430 X 3 3 VALUE 4.663702915655440 X 4 4 VALUE 4.620667030846704 X 5 5 VALUE 4.679200672505286 X 6 6 VALUE 1.617538218457437 X 1 2 X 2 1 VALUE 1.530250893440179 X 1 3 X 3 1 VALUE 3.546715188981548 X 1 4 X 4 1 VALUE -3.449854515258192 X 1 5 X 5 1 VALUE -1.573289443908528 X 1 6 X 6 1 VALUE 1.647925282478313 X 2 3 X 3 2 VALUE 3.546421968829633 X 2 4 X 4 2 VALUE -1.632671504988344 X 2 5 X 5 2 VALUE -3.578418676335090 X 2 6 X 6 2 VALUE 1.588011722661225 X 3 4 X 4 3 VALUE -3.450140420881972 X 3 5 X 5 3 VALUE -3.579009967331238 X 3 6 X 6 3 VALUE -2.576001700429044 X 4 5 X 5 4 VALUE -2.541120409739163 X 4 6 X 6 4 VALUE 2.567510428234816 X 5 6 X 6 5 BEGIN ALGEBRA; ! THEN SPECIFY THE CONSTRAINTS AMONG THE PARAMETERS OF INTERSTS (THETA1, THETA2,..., THETA6) H=A*A; ! THIS IS THETA 1 I=B*B; ! THIS IS THETA 2 J=C*C; ! THIS IS THETA 3 M=A*B; ! THIS IS THETA 4 N=A*C; ! THIS IS THETA 5 Q=B*C; ! THIS IS THETA 6 END ALGEBRA; END The following MX script lines were read for group 2 G2: MLE DATA NINPUT=6 NOBSERBATIONS=1 ! READ THE DATA FILE "MYDATA2.DAT", WHICH HAS ONE OBSERVATION Z, WHICH HAS 6 COLUMNS AND ONE ROW. RECTANGULAR FILE=MYDATA2.DAT Rectangular continuous data read initiated NOTE: Rectangular file contained 1 records with data that contained a total of 6 observations LABELS X1 X2 X3 X4 X5 X6 BEGIN MATRICES = GROUP 1; MEANS H_I_J_M_N_Q; ! THIS FORM MEAN PARAMETERS OF Z COVARIANCES X; ! X IS COVARIANCE OF Z, BUT SINCE WE FIXED THE VALUES IN X, IT WILL NOT ESTIMATE THE COVARIANCE. START 1 ALL; ! STARTING VALUES FOR ALL PARAMETERS OPTION NDECIMALS=9 ! DECIMAL PLACES !DROP 1 2 3 ! Use this script twice both for the full model and the null model. When you estimate the null model remove "!" END ! SIGN AT THE BEGINIG OF THIS LINE, THEN IT WILL DROP PARAMETER A, B, C. IF YOU ESTIMATE THE FULL MODEL, KEEP "!" SIGN The following MX script lines were read for group 3 ! then it won't drop any parameters and thus estimate the full model G3: GET FUNCTION ! THIS IS TO GET THE MAXIMUM LIKELIHOOD FUNCTION VALUE CALCULATION BEGIN MATRICES; Z FULL 1 1 = %F2 END MATRICES; END Summary of VL file data for group 2 X1 X2 X3 X4 Code 1.000000000 2.000000000 3.000000000 4.000000000 Number 1.000000000 1.000000000 1.000000000 1.000000000 Mean -0.753932012 0.539824235 -0.628481290 -1.243107783 Variance 0.000000000 0.000000000 0.000000000 0.000000000 Minimum -0.753932012 0.539824235 -0.628481290 -1.243107783 Maximum -0.753932012 0.539824235 -0.628481290 -1.243107783 X5 X6 Code 5.000000000 6.000000000 Number 1.000000000 1.000000000 Mean -0.883479696 -1.597953348 Variance 0.000000000 0.000000000 Minimum -0.883479696 -1.597953348 Maximum -0.883479696 -1.597953348 PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1: PARAMETER ESTIMATES MATRIX A This is a FULL matrix of order 1 by 1 1 1 1 MATRIX B This is a FULL matrix of order 1 by 1 1 1 2 MATRIX C This is a FULL matrix of order 1 by 1 1 1 3 MATRIX H This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX N This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Q This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 6 by 6 It has no free parameters specified GROUP NUMBER: 2 G2: MLE MATRIX A This is a FULL matrix of order 1 by 1 1 1 1 MATRIX B This is a FULL matrix of order 1 by 1 1 1 2 MATRIX C This is a FULL matrix of order 1 by 1 1 1 3 MATRIX H This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX I This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX J This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX M This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX N This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX Q This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a FULL matrix of order 6 by 6 It has no free parameters specified GROUP NUMBER: 3 G3: GET FUNCTION MATRIX Z This is a constrained FULL matrix of order 1 by 1 It has no free parameters specified Mx starting optimization; number of parameters = 3 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: PARAMETER ESTIMATES MATRIX A This is a FULL matrix of order 1 by 1 1 1 -0.321476627 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.180574145 MATRIX C This is a FULL matrix of order 1 by 1 1 1 -0.522928479 MATRIX H This is a computed FULL matrix of order 1 by 1 [=A*A] 1 1 0.103347222 MATRIX I This is a computed FULL matrix of order 1 by 1 [=B*B] 1 1 1.393755312 MATRIX J This is a computed FULL matrix of order 1 by 1 [=C*C] 1 1 0.273454194 MATRIX M This is a computed FULL matrix of order 1 by 1 [=A*B] 1 1 -0.379526994 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A*C] 1 1 0.168109283 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=B*C] 1 1 -0.617355842 MATRIX X This is a FULL matrix of order 6 by 6 1 2 3 4 5 1 7.710474640 1.617538218 1.530250893 3.546715189 -3.449854505 2 1.617538218 7.709260663 1.647925282 3.546421969 -1.632671495 3 1.530250893 1.647925282 7.711691774 1.588011723 -3.450140419 4 3.546715189 3.546421969 1.588011723 4.663702916 -2.576001700 5 -3.449854505 -1.632671495 -3.450140419 -2.576001700 4.620667031 6 -1.573289436 -3.578418664 -3.579009953 -2.541120390 2.567510428 6 1 -1.573289436 2 -3.578418664 3 -3.579009953 4 -2.541120390 5 2.567510428 6 4.679200673 GROUP NUMBER: 2 G2: MLE MATRIX A This is a FULL matrix of order 1 by 1 1 1 -0.321476627 MATRIX B This is a FULL matrix of order 1 by 1 1 1 1.180574145 MATRIX C This is a FULL matrix of order 1 by 1 1 1 -0.522928479 MATRIX H This is a computed FULL matrix of order 1 by 1 [=A*A] 1 1 0.103347222 MATRIX I This is a computed FULL matrix of order 1 by 1 [=B*B] 1 1 1.393755312 MATRIX J This is a computed FULL matrix of order 1 by 1 [=C*C] 1 1 0.273454194 MATRIX M This is a computed FULL matrix of order 1 by 1 [=A*B] 1 1 -0.379526994 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A*C] 1 1 0.168109283 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=B*C] 1 1 -0.617355842 MATRIX X This is a FULL matrix of order 6 by 6 1 2 3 4 5 1 7.710474640 1.617538218 1.530250893 3.546715189 -3.449854505 2 1.617538218 7.709260663 1.647925282 3.546421969 -1.632671495 3 1.530250893 1.647925282 7.711691774 1.588011723 -3.450140419 4 3.546715189 3.546421969 1.588011723 4.663702916 -2.576001700 5 -3.449854505 -1.632671495 -3.450140419 -2.576001700 4.620667031 6 -1.573289436 -3.578418664 -3.579009953 -2.541120390 2.567510428 6 1 -1.573289436 2 -3.578418664 3 -3.579009953 4 -2.541120390 5 2.567510428 6 4.679200673 GROUP NUMBER: 3 G3: GET FUNCTION MATRIX Z This is a constrained FULL matrix of order 1 by 1 1 1 21.882200914 Your model has 3 estimated parameters and 6 Observed statistics -2 times log-likelihood of data >>> 21.882 Degrees of freedom >>>>>>>>>>>>>>>> 3 This problem used 0.1% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.00 Execution 0: 0: 0: 0.01 TOTAL 0: 0: 0: 0.01 Total number of warnings issued: 0 ______________________________________________________________________________ ! to run this scirpt in unix, use the command "mx < myscriptF.mx > outfile.FULL.mxo" (outfile.FULL.mxo can be any output file name) ! Note you also need input data file, "myData2.dat" to run this script ______________________________________________________________________________