! Boulder 2003 ! saturated model ! using raw data input file ! MZM, DZM, MZF, DZF, DOSMF, DOSFM twins ! Mx homepage: http://views.vcu.edu/mx #define nvar 2 ! For every instance of the word nvar in the script, ! Mx reads 1. This corresponds to one phenotype to be analysed !#Define twonvar 2 ! twice nvar. will be used to calculate means G1: MZM twins, datagroup Data NInput_vars= 14 NGroups=1 ! nr of inputvars per family Missing=-1.00 ! missing values = -1.00 REctangular file=all.dat ! read in raw datafile, must be in same directory LABELS ! tell Mx what is found in the datafile nih_id country zygot6 rear ageA sexA ldllnA apoblnA apoelnA ageB sexB ldllnB apoblnB apoelnB Select if country = 1 ; ! select Dutch twins, younger sample (Dorret) Select if zygot6 = 1 ; ! select MZM's Select ldllnA ldllnB ; ! select one phenotype per twin member Matrices = group 1 G Full 1 nvar free ! mean X Lower nvar nvar Free ! covariance end matrices; Begin Algebra; Y = \stnd(X*X'); ! twin correlation End Algebra; Means G ; ! model for means Covariances X*X' ; ! model for MZM variance/covariances st 15 G 1 1 - G 1 nvar st 2.5 X 1 1 X 2 2 st 0.5 X 2 1 Options RSiduals End G2: DZM twins, datagroup Data NInput_vars= 14 NGroups=1 ! nr of inputvars per family Missing=-1.00 ! missing values = -1.00 REctangular file=all.dat ! read in raw datafile, must be in same directory LABELS ! tell Mx what is found in the datafile nih_id country zygot6 rear ageA sexA ldllnA apoblnA apoelnA ageB sexB ldllnB apoblnB apoelnB Select if country = 1 ; ! select Dutch twins, younger sample (Dorret) Select if zygot6 = 2 ; ! select MZM's Select ldllnA ldllnB ; ! select one phenotype per twin member Matrices = group 1 G Full 1 nvar free ! mean X Lower nvar nvar Free ! covariance end matrices; Begin Algebra; Y = \stnd(X*X'); ! twin correlation End Algebra; Means G ; ! model for means Covariances X*X' ; ! model for MZM variance/covariances st 15 G 1 1 - G 1 nvar st 2.5 X 1 1 X 2 2 st 0.5 X 2 1 Options RSiduals End G3: MZF twins, datagroup Data NInput_vars= 14 NGroups=1 ! nr of inputvars per family Missing=-1.00 ! missing values = -1.00 REctangular file=all.dat ! read in raw datafile, must be in same directory LABELS ! tell Mx what is found in the datafile nih_id country zygot6 rear ageA sexA ldllnA apoblnA apoelnA ageB sexB ldllnB apoblnB apoelnB Select if country = 1 ; ! select Dutch twins, younger sample (Dorret) Select if zygot6 = 3 ; ! select MZM's Select ldllnA ldllnB ; ! select one phenotype per twin member Matrices = group 1 G Full 1 nvar free ! mean X Lower nvar nvar Free ! covariance end matrices; Begin Algebra; Y = \stnd(X*X'); ! twin correlation End Algebra; Means G ; ! model for means Covariances X*X' ; ! model for MZM variance/covariances st 15 G 1 1 - G 1 nvar st 2.5 X 1 1 X 2 2 st 0.5 X 2 1 Options RSiduals End G4: DZF twins, datagroup Data NInput_vars= 14 NGroups=1 ! nr of inputvars per family Missing=-1.00 ! missing values = -1.00 REctangular file=all.dat ! read in raw datafile, must be in same directory LABELS ! tell Mx what is found in the datafile nih_id country zygot6 rear ageA sexA ldllnA apoblnA apoelnA ageB sexB ldllnB apoblnB apoelnB Select if country = 1 ; ! select Dutch twins, younger sample (Dorret) Select if zygot6 = 4 ; ! select MZM's Select ldllnA ldllnB ; ! select one phenotype per twin member Matrices = group 1 G Full 1 nvar free ! mean X Lower nvar nvar Free ! covariance end matrices; Begin Algebra; Y = \stnd(X*X'); ! twin correlation End Algebra; Means G ; ! model for means Covariances X*X' ; ! model for MZM variance/covariances st 15 G 1 1 - G 1 nvar st 2.5 X 1 1 X 2 2 st 0.5 X 2 1 Options RSiduals End G5: DOSMF twins, datagroup Data NInput_vars= 14 NGroups=1 ! nr of inputvars per family Missing=-1.00 ! missing values = -1.00 REctangular file=all.dat ! read in raw datafile, must be in same directory LABELS ! tell Mx what is found in the datafile nih_id country zygot6 rear ageA sexA ldllnA apoblnA apoelnA ageB sexB ldllnB apoblnB apoelnB Select if country = 1 ; ! select Dutch twins, younger sample (Dorret) Select if zygot6 = 5 ; ! select MZM's Select ldllnA ldllnB ; ! select one phenotype per twin member Matrices = group 1 G Full 1 nvar free ! mean X Lower nvar nvar Free ! covariance end matrices; Begin Algebra; Y = \stnd(X*X'); ! twin correlation End Algebra; Means G ; ! model for means Covariances X*X' ; ! model for MZM variance/covariances st 15 G 1 1 - G 1 nvar st 2.5 X 1 1 X 2 2 st 0.5 X 2 1 Options RSiduals End G6: DOSFM twins, datagroup Data NInput_vars= 14 NGroups=1 ! nr of inputvars per family Missing=-1.00 ! missing values = -1.00 REctangular file=all.dat ! read in raw datafile, must be in same directory LABELS ! tell Mx what is found in the datafile nih_id country zygot6 rear ageA sexA ldllnA apoblnA apoelnA ageB sexB ldllnB apoblnB apoelnB Select if country = 1 ; ! select Dutch twins, younger sample (Dorret) Select if zygot6 = 6 ; ! select MZM's Select ldllnA ldllnB ; ! select one phenotype per twin member Matrices = group 1 G Full 1 nvar free ! mean X Lower nvar nvar Free ! covariance end matrices; Begin Algebra; Y = \stnd(X*X'); ! twin correlation End Algebra; Means G ; ! model for means Covariances X*X' ; ! model for MZM variance/covariances st 15 G 1 1 - G 1 nvar st 2.5 X 1 1 X 2 2 st 0.5 X 2 1 Options RSiduals End