Univariate AnalysisTestsEquality TestsSpecific Equality TestsMultiple FitUnivariate Genetic AnalysisACE ModelACE Model + MeansSlide 9! Estimate variance components - ACED model ! OZ BMI data - younger females! Estimate variance components - ACED model ! OZ BMI data - younger females II! Estimate variance components - ACED model ! OZ BMI data - younger females III! Estimate variance components - ACED model ! OZ BMI data - younger females IVSubmodels: ozbmifyaces.mxSlide 15Submodels: ozbmifyaces2.mxGoodness-of-FitParameter EstimatesUnivariate AnalysisHGEN619 class 2006TestsSaturated modelEquality of meansIs m1 = m2 ?Is m1MZ = m1DZ = m2MZ = m2DZ ?Equality of variancesIs v1 = v2 ?Is v1MZ = v1DZ = v2MZ = v2DZ ?Equality TestsMain ScriptLast Group.... Option Multiple IssatEndSave ozbmisat.mxs! equate means and variances Equate M 1 1 1 M 1 1 2 M 2 1 1 M 2 1 2 Equate X 1 1 1 X 1 2 2 X 2 1 1 X 2 2 2EndSpecific Equality TestsGet ozbmisat.mxs! equate means within zygosity groups Equate M 1 1 1 M 1 1 2 Equate M 2 1 1 M 2 1 2End! equate means across zygosity groups Equate M 1 1 1 M 1 1 2 M 2 1 1 M 2 1 2End! equate variances within zygosity groups Equate X 1 1 1 X 1 2 2 Equate X 2 1 1 X 2 2 2End! equate variances across zygosity groups Equate X 1 1 1 X 1 2 2 X 2 1 1 X 2 2 2EndMultiple FitMZ (group 1) DZ (group 2) parm1 m2 v1 cov v2 m3 m4 v3 cov v4Full 1 2 3 4 5 6 7 8 9 10 10Save filename.mxs: 10 free parameters 10I 1 1 3 4 3 1 1 3 9 3 4Get filename.mxs: back to 10 free parameters 10II 1 1 3 4 5 6 6 8 9 10 8III 1 1 3 4 5 1 1 8 9 10 7IV 1 1 3 4 3 1 1 8 9 8 5V 1 1 3 4 3 1 1 3 9 3 4Univariate Genetic AnalysisSaturated ModelsFree variances, covariancesFree meansUnivariate ModelsVariances partitioned in a, c/d and eFree means (or not)ACE ModelT2AEA Ca ac eT11 111Ee1Cc11 or .51ACE Model + MeansT2AEA Ca ac eT11 111Ee1Cc11 or .511m1 or m3 m2 or m4TestsACE modelIs a significant ? -> CE modelIs c significant ? -> AE modelIs there significant family resemblance ? -> E modelADE model! Estimate variance components - ACED model! OZ BMI data - younger females#NGroups 4#define nvar 1#define nvar2 2Title 1: Model Parameters Calculation Begin Matrices; X Lower nvar nvar Free ! a Y Lower nvar nvar ! c Z Lower nvar nvar Free ! e W Lower nvar nvar Free ! d H Full 1 1 ! 0.5 Q Full 1 1 ! 0.25 End Matrices; Matrix H .5 Matrix Q .25 Label Row X add_gen Label Row Y com_env Label Row Z spec_env Label Row W dom_genBegin Algebra; A= X*X'; ! a^2 C= Y*Y'; ! c^2 E= Z*Z'; ! e^2 D= W*W'; ! d^2 End Algebra;Endozbmifyace.mx! Estimate variance components - ACED model! OZ BMI data - younger females II Title 2: MZ data #include ozbmi.dat Select if zyg =1 Select bmi1 bmi2 ; Begin Matrices = Group 1; M Full 1 nvar2 Free Means M; Covariance A+C+E+D | A+C+D _ A+C+D | A+C+E+D; Option RSiduals;EndTitle 3: DZ data #include ozbmi.dat Select if zyg =3 Select bmi1 bmi2 ; Begin Matrices = Group 1; M Full 1 nvar2 Free End Matrices; Means M; Covariance A+C+E+D | H@A+C+Q@D _ H@A+C+Q@D | A+C+E+D; Option RSidualsEndozbmifyace.mx! Estimate variance components - ACED model! OZ BMI data - younger females IIITitle 4: Standardization Calculation Begin Matrices = Group 1; End Matrices; Start .6 all Start 20 M 2 1 1 - M 2 1 nvar2 Start 20 M 3 1 1 - M 3 1 nvar2 Begin Algebra; V=A+C+E+D; ! total variance P=A|C|E|D; ! concatenate parameter estimates S=P@V~; ! standardized parameter estimates End Algebra;!ADE model Interval S 1 1 - S 1 4 Option NDecimals=4 Option Sat=4055.935,1767Endozbmifyace.mx! Estimate variance components - ACED model! OZ BMI data - younger females IVTitle 4: Standardization Calculation Begin Matrices = Group 1; End Matrices; Start .6 all Start 20 M 2 1 1 - M 2 1 2 Start 20 M 3 1 1 - M 3 1 2 Begin Algebra; V=A+C+E+D; P=A|C|E|D; S=P@V~; End Algebra;!ADE model Interval S 1 1 - S 1 4 Option NDecimals=4 Option Sat=4055.935,1767 Option MultipleEnd!AE model Drop W 1 1 1End!ACE model Free Y 1 1 1End!CE model Drop X 1 1 1End!E model Drop Y 1 1 1Endozbmifyaces.mxSubmodels: ozbmifyaces.mxMatrix / ModelX (a) Y (c) Z (e) W (d) Cov NPMean NPNP DFSat 6 4 10ADE Free Free Free 4 7 3AE Free Free Drop 4 6 4ACE Free Free Free 4 7 3CE Drop Free Free 4 6 4E Drop Free 4 5 5! Estimate variance components - ACED model! OZ BMI data - younger females IV.....!ADE model Interval S 1 1 - S 1 4 Option NDecimals=4 Option Sat=4055.935,1767 Option MultipleEnd!Save ozbmify.mxs Option IssatEnd!AE model Drop W 1 1 1End!ACE model Free Y 1 1 1 Option Sat=4055.935,1767End Option IssatEnd!CE model Drop X 1 1 1End!E model Drop Y 1 1 1Endozbmifyaces2.mxSubmodels: ozbmifyaces2.mxMatrix / ModelX (a) Y (c) Z (e) W (d) Cov NPMean NPNP Sat DFSat # 6 4 10ADE $ Free Free Free 4 7 10# 3AE Free Free Drop 4 6 7 $ 1ACE & Free Free Free 4 7 10# 3CE Drop Free Free 4 6 7 & 1E Drop Free 4 5 7 & 2Goodness-of-Fit-2LL df2df p AIC2df pSatADE 3AE 4 1ACE 3CE 4 1E 5 2Parameter Estimatesa c e d
View Full Document