** Mx startup successful ** **MX-OSF version 1.50c** ! Script to estimate haplotype frequencies The following MX script lines were read for group 1 #NGROUPS 2 Note: #NGroup set number of groups to 2 TITLE ESTIMATE FREQS DATA NI=1 BEGIN MATRICES; I UNIT 1 1 P FULL 3 3 FREE ! EXPECTED PROPORTIONS O FULL 3 3 ! OBSERVED GENOTYPE FREQS T FULL 1 1 END MATRICES; MATRIX O ! OBSERVED FREQS 10 15 5 10 50 13 3 13 10 START .111 ALL MATRIX T 2 BEGIN ALGEBRA; END ALGEBRA; BOUND 0 1 P 1 1 TO P 3 3 COMPUTE -T. \SUM(O.\LN(P)); OPTION USER END GROUP The following MX script lines were read for group 2 CONSTRAIN HAPLOTYPE FREQS TO SUM TO 1 CONSTRAINT BEGIN MATRICES = GROUP 1; END MATRICES; CONSTRAIN I=\SUM(P); OPTION RS END GROUP; PARAMETER SPECIFICATIONS GROUP NUMBER: 1 Title estimate freqs MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX O This is a FULL matrix of order 3 by 3 It has no free parameters specified MATRIX P This is a FULL matrix of order 3 by 3 1 2 3 1 1 2 3 2 4 5 6 3 7 8 9 MATRIX T This is a FULL matrix of order 1 by 1 It has no free parameters specified GROUP NUMBER: 2 Constrain haplotype freqs to sum to 1 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX O This is a FULL matrix of order 3 by 3 It has no free parameters specified MATRIX P This is a FULL matrix of order 3 by 3 1 2 3 1 1 2 3 2 4 5 6 3 7 8 9 MATRIX T This is a FULL matrix of order 1 by 1 It has no free parameters specified MX PARAMETER ESTIMATES GROUP NUMBER: 1 Title estimate freqs MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX O This is a FULL matrix of order 3 by 3 1 2 3 1 10.0000 15.0000 5.0000 2 10.0000 50.0000 13.0000 3 3.0000 13.0000 10.0000 MATRIX P This is a FULL matrix of order 3 by 3 1 2 3 1 0.0775 0.1163 0.0388 2 0.0775 0.3876 0.1008 3 0.0233 0.1008 0.0775 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.0000 GROUP NUMBER: 2 Constrain haplotype freqs to sum to 1 MATRIX I This is a UNIT matrix of order 1 by 1 MATRIX O This is a FULL matrix of order 3 by 3 1 2 3 1 10.0000 15.0000 5.0000 2 10.0000 50.0000 13.0000 3 3.0000 13.0000 10.0000 MATRIX P This is a FULL matrix of order 3 by 3 1 2 3 1 0.0775 0.1163 0.0388 2 0.0775 0.3876 0.1008 3 0.0233 0.1008 0.0775 MATRIX T This is a FULL matrix of order 1 by 1 1 1 2.0000 CONSTRAINT VALUES (should be near zero) 1 1 -1.3136E-07 Your model has 9 estimated parameters and 1 Observed statistics Observed statistics include 1 constraints. User defined function value = 487.169 'Degrees of freedom' >>>>>>>>>>>>>>>> -8 This problem used 0.1% 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 ______________________________________________________________________________ ______________________________________________________________________________