** 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 No free parameters specified MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: PARAMETER ESTIMATES MATRIX A This is a FULL matrix of order 1 by 1 1 1 0.000000000 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.000000000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.000000000 MATRIX H This is a computed FULL matrix of order 1 by 1 [=A*A] 1 1 0.000000000 MATRIX I This is a computed FULL matrix of order 1 by 1 [=B*B] 1 1 0.000000000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=C*C] 1 1 0.000000000 MATRIX M This is a computed FULL matrix of order 1 by 1 [=A*B] 1 1 0.000000000 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A*C] 1 1 0.000000000 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=B*C] 1 1 0.000000000 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.000000000 MATRIX B This is a FULL matrix of order 1 by 1 1 1 0.000000000 MATRIX C This is a FULL matrix of order 1 by 1 1 1 0.000000000 MATRIX H This is a computed FULL matrix of order 1 by 1 [=A*A] 1 1 0.000000000 MATRIX I This is a computed FULL matrix of order 1 by 1 [=B*B] 1 1 0.000000000 MATRIX J This is a computed FULL matrix of order 1 by 1 [=C*C] 1 1 0.000000000 MATRIX M This is a computed FULL matrix of order 1 by 1 [=A*B] 1 1 0.000000000 MATRIX N This is a computed FULL matrix of order 1 by 1 [=A*C] 1 1 0.000000000 MATRIX Q This is a computed FULL matrix of order 1 by 1 [=B*C] 1 1 0.000000000 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 22.642702709 Your model has 0 estimated parameters and 6 Observed statistics -2 times log-likelihood of data >>> 22.643 Degrees of freedom >>>>>>>>>>>>>>>> 6 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.00 TOTAL 0: 0: 0: 0.00 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 ______________________________________________________________________________