** Mx startup successful ** **MX-OSF version 1.52b** ! Simulate the data ! 30% additive genetic ! 20% common environment ! 50% nonshared environment The following MX script lines were read for group 1 #NGROUPS 3 Note: #NGroup set number of groups to 3 G1: MODEL PARAMETERS CALCULATION BEGIN MATRICES; X LOWER 1 1 FIXED Y LOWER 1 1 FIXED Z LOWER 1 1 FIXED END MATRICES; MATRIX X 0.5477 !UNSQUARED ESTIMATE FOR A OF .3 MATRIX Y 0.4472 !UNSQUARED ESTIMATE FOR C OF .2 MATRIX Z 0.7071 !UNSQUARED ESTIMATE FOR E OF .5 BEGIN ALGEBRA; A = X*X' ; C = Y*Y' ; E = Z*Z' ; END ALGEBRA; END The following MX script lines were read for group 2 G2: MZ TWIN PAIRS CALCULATION MATRICES = GROUP 1 COVARIANCES A+C+E | A+C _ A+C | A+C+E / OPTIONS MX%E=MZSIM.COV !STORES THE MZ COV. MATRIX IN MZSIM.COV END The following MX script lines were read for group 3 G3: DZ TWIN PAIRS CALCULATION MATRICES = GROUP 1 H FULL 1 1 COVARIANCES A+C+E | H@A+C _ H@A+C | A+C+E / MATRIX H 0.5 OPTIONS MX%E=DZSIM.COV !STORES THE DZ COV. MATRIX IN MZSIM.COV END PARAMETER SPECIFICATIONS No free parameters specified MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1: model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.3000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.2000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.5000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5477 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4472 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7071 GROUP NUMBER: 2 G2: MZ twin pairs Estimates of matrix % will be written to file: mzsim.cov MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.3000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.2000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.5000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5477 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4472 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7071 GROUP NUMBER: 3 G3: DZ twin pairs Estimates of matrix % will be written to file: dzsim.cov MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.3000 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.2000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.5000 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.5477 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.4472 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7071 Your model has 0 estimated parameters and 0 Observed statistics Chi-squared fit of model >>>>>>> 0.000 Degrees of freedom >>>>>>>>>>>>> 0 Probability incalculable Akaike's Information Criterion > 0.000 RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.000 This problem used 0.1% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.02 Execution 0: 0: 0: 0.01 TOTAL 0: 0: 0: 0.03 Total number of warnings issued: 0 ______________________________________________________________________________ ! Second part of script ! Fit the wrong model to the simulated data ! to calculate power The following MX script lines were read for group 1 #NGROUPS 3 Note: #NGroup set number of groups to 3 G1 : MODEL PARAMETERS CALCULATION BEGIN MATRICES; X LOWER 1 1 FREE Y LOWER 1 1 FIXED !DROPPED C FROM MODEL FOR TEST OF POWER Z LOWER 1 1 FREE END MATRICES; BEGIN ALGEBRA; A = X*X' ; C = Y*Y' ; E = Z*Z' ; END ALGEBRA; END The following MX script lines were read for group 2 G2 : MZ TWINS DATA NINPUT_VARS=2 NOBSERVATIONS=350 !NECESSARY TO SPECIFY NOBSERVATIONS FOR BASELINE CMATRIX FULL FILE=MZSIM.COV !SUPPLIES THE COVARIANCE MATRIX FROM ABOVE FOR MZ MATRICES= GROUP 1 COVARIANCES A+C+E | A+C _ A+C | A+C+E / OPTION RSIDUALS END The following MX script lines were read for group 3 G3 : DZ TWINS DATA NINPUT_VARS=2 NOBSERVATIONS=350 !NECESSARY TO SPECIFY NOBSERVATIONS FOR BASELINE CMATRIX FULL FILE=DZSIM.COV !SUPPLIES THE COVARIANCE MATRIX FROM ABOVE FOR MZ MATRICES= GROUP 1 H FULL 1 1 COVARIANCES A+C+E | H@A+C _ H@A+C | A+C+E / MATIX H 0.5 OPTION RSIDUALS ! Power for alpha = 0.05 and 1 df OPTION POWER= 0.05,1 !ASKS FOR POWER LEVELS WHERE ALPHA=.05 END PARAMETER SPECIFICATIONS GROUP NUMBER: 1 G1 : model parameters MATRIX A 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 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 GROUP NUMBER: 2 G2 : MZ twins MATRIX A 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 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 GROUP NUMBER: 3 G3 : DZ twins MATRIX A 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 1 by 1 It has no free parameters specified MATRIX E This is a computed FULL matrix of order 1 by 1 It has no free parameters specified MATRIX H This is a FULL matrix of order 1 by 1 It has no free parameters specified MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 1 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 It has no free parameters specified MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 2 Mx starting optimization; number of parameters = 2 MX PARAMETER ESTIMATES GROUP NUMBER: 1 G1 : model parameters MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5186 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.4765 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7201 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6903 GROUP NUMBER: 2 G2 : MZ twins MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5186 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.4765 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7201 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6903 OBSERVED COVARIANCE MATRIX 1 2 1 1.0000 2 0.5000 1.0000 EXPECTED COVARIANCE MATRIX 1 2 1 0.9950 2 0.5186 0.9950 RESIDUAL MATRIX 1 2 1 0.0049 2 -0.0186 0.0049 Function value of this group: 0.4263 Where the fit function is Maximum Likelihood GROUP NUMBER: 3 G3 : DZ twins MATRIX A This is a computed FULL matrix of order 1 by 1 [=X*X'] 1 1 0.5186 MATRIX C This is a computed FULL matrix of order 1 by 1 [=Y*Y'] 1 1 0.0000 MATRIX E This is a computed FULL matrix of order 1 by 1 [=Z*Z'] 1 1 0.4765 MATRIX H This is a FULL matrix of order 1 by 1 1 1 0.5000 MATRIX X This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.7201 MATRIX Y This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.0000 MATRIX Z This is a LOWER TRIANGULAR matrix of order 1 by 1 1 1 0.6903 OBSERVED COVARIANCE MATRIX 1 2 1 1.0000 2 0.3500 1.0000 EXPECTED COVARIANCE MATRIX 1 2 1 0.9950 2 0.2593 0.9950 RESIDUAL MATRIX 1 2 1 0.0049 2 0.0907 0.0049 Function value of this group: 3.5386 Where the fit function is Maximum Likelihood Your model has 2 estimated parameters and 6 Observed statistics Chi-squared fit of model >>>>>>> 3.965 Degrees of freedom >>>>>>>>>>>>> 4 Probability >>>>>>>>>>>>>>>>>>>> 0.411 Akaike's Information Criterion > -4.035 RMSEA >>>>>>>>>>>>>>>>>>>>>>>>>> 0.023 Power of this test, at the 0.0500 significance level with 1. df is 0.512506 Based on your combined observed sample size of 700. The following sample sizes would be required to reject the hypothesis: Power Total N .25 291. .50 678. .75 1225. .80 1386. .90 1855. .95 2294. .99 3244. This problem used 0.0% of my workspace Task Time elapsed (DD:HH:MM:SS) Reading script & data 0: 0: 0: 0.02 Execution 0: 0: 0: 0.06 TOTAL 0: 0: 0: 0.08 Total number of warnings issued: 0 ______________________________________________________________________________ ______________________________________________________________________________