! satmf.mx ! ! Univariate saturated model for one continuous phenotype ! Raw data ! Child Behavior Checklist, NTR MZM DZM, MZF and DZF twins ! ! Mx homepage: http://views.vcu.edu/mx #NGroups 5 #define nvar 1 #define nsib 2 G1: 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; X Symm nsib nsib Free ! covariances I Iden nsib nsib M Full nvar nsib Free ! means End Matrices; Start 2 X 1 1 X 2 2 ! starting values for variances Start 0.5 M 1 1 M 1 2 ! starting values for means Begin Algebra; O= \sqrt(I.X)~&X; ! MZM correlation End Algebra; Means M; ! model for MZM means Covariances X; ! model for MZM variance/covariances ! Interval @95 O 2 1 Option RSiduals End G2: 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; Y Symm nsib nsib Free ! covariances I Iden nsib nsib N Full nvar nsib Free ! means End Matrices; Start 2 Y 1 1 Y 2 2 ! starting values for variances Start 0.5 N 1 1 N 1 2 ! starting values for means Begin Algebra; P= \sqrt(I.Y)~&Y; ! DZM correlation End Algebra; Means N; ! model for DZM means Covariances Y; ! model for DZM variance/covariances ! Interval @95 P 2 1 Option RSiduals Option Multiple Issat End G3: 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; X Symm nsib nsib Free ! covariances I Iden nsib nsib M Full nvar nsib Free ! means End Matrices; Start 2 X 1 1 X 2 2 ! starting values for variances Start 0.5 M 1 1 M 1 2 ! starting values for means Begin Algebra; O= \sqrt(I.X)~&X; ! MZF correlation End Algebra; Means M; ! model for MZF means Covariances X; ! model for MZF variance/covariances ! Interval @95 O 2 1 Option RSiduals End G4: 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; Y Symm nsib nsib Free ! covariances I Iden nsib nsib N Full nvar nsib Free ! means End Matrices; Start 2 Y 1 1 Y 2 2 ! starting values for variances Start 0.5 N 1 1 N 1 2 ! starting values for means Begin Algebra; P= \sqrt(I.Y)~&Y; ! DZF correlation End Algebra; Means N; ! model for DZF means Covariances Y; ! model for DZF variance/covariances ! Interval @95 P 2 1 Option RSiduals Option Multiple Issat End G5: male-female DZ twin pairs Data NInput_vars=5 Missing=-1.00 Rectangular File=Agg10.rec Labels ZYG RB10A AGG10A RB10B AGG10B Select if zyg =5 ; ! select DZMF twins Select AGG10A AGG10B ; Begin Matrices; Y Symm nsib nsib Free ! covariances I Iden nsib nsib N Full nvar nsib Free ! means End Matrices; Start 2 Y 1 1 Y 2 2 ! starting values for variances Start 0.5 N 1 1 N 1 2 ! starting values for means Begin Algebra; P= \sqrt(I.Y)~&Y; ! DZMF correlation End Algebra; Means N; ! model for DZMF means Covariances Y; ! model for DZMF variance/covariances ! Interval @95 P 2 1 Option RSiduals Option Multiple Issat End Save satmf5.mxs ! equate means same-sex male Equate M 1 1 1 M 1 1 2 N 2 1 1 N 2 1 2 End ! equate means same-sex and opposite-sex male Equate M 1 1 1 N 5 1 1 End ! equate means same-sex female Equate M 3 1 1 M 3 1 2 N 4 1 1 N 4 1 2 End ! equate means same-sex and opposite-sex female Equate M 3 1 1 N 5 1 2 End Get satmf5.mxs ! equate means same-sex female Equate M 3 1 1 M 3 1 2 N 4 1 1 N 4 1 2 End ! equate means same-sex and opposite-sex female Equate M 3 1 1 N 5 1 2 End Get satmf5.mxs ! equate variances same-sex male Equate X 1 1 1 X 1 2 2 Y 2 1 1 Y 2 2 2 End ! equate variances same-sex female Equate X 3 1 1 X 3 2 2 Y 4 1 1 Y 4 2 2 End ! equate variances same-sex and opposite-sex male Equate X 1 1 1 Y 5 1 1 End ! equate variances same-sex and opposite-sex female Equate X 3 1 1 Y 5 2 2 End ! equate all variances Equate X 1 1 1 X 3 1 1 End