! unim.mx ! ! Univariate Model fitting for ADE, ACE, AE, CE and E ! Raw data ! Child Behavior Checklist, NTR MZM and DZM twins ! ! Mx homepage: http://views.vcu.edu/mx #NGroups 3 #define nvar 1 ! variable AGG or RB, mother rating, at age 10, males #define nsib 2 ! twin-1 & twin-2 G1: Model Parameters Calculation Begin Matrices; X Lower nvar nvar Free ! additive genetic structure Y Lower nvar nvar Free ! common environmental structure Z Lower nvar nvar Free ! specific environmental structure W Lower nvar nvar ! dominance structure H Full 1 1 ! scalar fixed @ .5 Q Full 1 1 ! scalar fixed @ .25 End Matrices; Matrix H .5 Matrix Q .25 Start .5 all ! starting values for all free parameters Begin Algebra; A= X*X' ; ! additive genetic variance C= Y*Y' ; ! common environmental variance E= Z*Z' ; ! specific environmental variance D= W*W' ; ! dominance variance V= A+C+E+D; ! total variance P= A|C|E|D; ! put parameter estimates in one matrix S= P@V~; ! standardized parameter estimates End Algebra; !Interval @95 S 1 1-S 1 3 ! confidence intervals End G2: male MZ twin pairs Data NInput_vars=5 Missing=-1.00 Rectangular File=Agg10.rec Labels ZYG RB10A AGG10A RB10B AGG10B Select if zyg =1 ; ! select MZM twins Select AGG10A AGG10B ; Begin Matrices = Group 1; M Full nsib nvar ! means End Matrices; Specify M 100 100 Means M ; ! model for means Covariances ! model for covariances A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; Options RSiduals End G3: male DZ twin pairs Data NInput_vars=5 Missing=-1.00 Rectangular File=Agg10.rec Labels ZYG RB10A AGG10A RB10B AGG10B Select if zyg =2 ; ! select DZM twins Select AGG10A AGG10B ; Begin Matrices = Group 1; M Full nsib nvar =M2 ! means End Matrices; Specify M 100 100 Means M ; ! model for means Covariances ! model for covariances A+C+D+E | H@A+C+Q@D _ H@A+C+Q@D | A+C+D+E ; Option RSiduals Option Multiple Issat End Save unim.mxs !AE model Drop Y 1 1 1 End Get unim.mxs !CE model Drop X 1 1 1 End !E model Drop Y 1 1 1 End