! Univariate Model fitting for ADE, ACE, AE, CE and E ! NTR Child Behavior Checklist ! #define nvar 1 ! variable AGG or RB, mother rating, at age 10, males #define nsib 2 ! twin-1 & twin-2 #NGroups 6 G1: Model Parameters Boys 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 !Free ! 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 I= P@V~; ! standardized parameter estimates End Algebra; Interval @95 I 1 1-I 1 3 ! confidence intervals End G2: Model Parameters Girls 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 !Free ! 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 J= P@V~; ! standardized parameter estimates End Algebra; Interval @95 J 1 1-J 1 3 ! confidence intervals End G3: 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 free ! means End Matrices; Means M ; ! model for means Covariances ! model for covariances A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; st 1.9 M 1 1 - M 2 1 Options RSiduals End G4: male DZ twin pairs Data NInput_vars=5 Missing=-1.00 ectangular 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 free ! means End Matrices; 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 ; st 1.9 M 1 1 - M 2 1 EQ M 3 1 1 M 3 2 1 M 4 1 1 M 4 2 1 Option RSiduals Options Multiple End G5: female MZ twin pairs Data NInput_vars=5 Missing=-1.00 Rectangular File=Agg10.rec Labels ZYG RB10A AGG10A RB10B AGG10B Select if zyg =3 ; ! select MZF twins Select AGG10A AGG10B ; Begin Matrices = Group 2; M Full nsib nvar free ! means End Matrices; !Specify M 200 200 Means M ; ! model for means Covariances ! model for covariances A+C+D+E | A+C+D _ A+C+D | A+C+D+E ; st 1.9 M 1 1 - M 2 1 Options RSiduals End G6: female DZ twin pairs Data NInput_vars=5 Missing=-1.00 Rectangular File=Agg10.rec Labels ZYG RB10A AGG10A RB10B AGG10B Select if zyg =4 ; ! select DZM twins Select AGG10A AGG10B ; Begin Matrices = Group 2; M Full nsib nvar free ! means End Matrices; 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 ; st 1.9 M 1 1 - M 2 1 EQ M 5 1 1 M 5 2 1 M 6 1 1 M 6 2 1 Option RSiduals Options Multiple issat End !no sex differences EQ X 1 1 1 X 2 1 1 EQ Y 1 1 1 Y 2 1 1 EQ Z 1 1 1 Z 2 1 1 end !Save aced.mxs !Get aced.mxs ! ADE model ! Drop Y 1 1 1 !End ! AE model ! Drop Y 1 1 1 ! Drop Y 2 1 1 !End !Get aced.mxs ! ACE model ! Drop W 1 1 1 !End ! CE mmodel ! Drop X 1 1 1 ! Drop X 2 1 1 !End ! E model ! Drop Y 1 1 1 !End