! Estimate genetic and environmental components - ACED model
! OZ BMI data - younger females

#NGroups 4
#define nvar 1
#define nvar2 2

Title G1: Model Parameters
 Calculation
 Begin Matrices;
  X Lower nvar nvar Free  ! additive genetic path, a
  Y Lower nvar nvar       ! common environmental path, c
  Z Lower nvar nvar Free  ! specific environmental path, e
  W Lower nvar nvar Free  ! dominance genetic path, d
  H Full 1 1              ! scalar, 0.5
  Q Full 1 1              ! scalar, 0.25
 End Matrices;
  Label Row X add_gen
  Label Row Y comm_env
  Label Row Z spec_env
  Label Row W dom_gen
  Matrix H .5
  Matrix Q .25
 Begin Algebra;
 End Algebra;
End Group

Title G2: MZ data
 #include ozbmi.dat
 Select if zyg=1
 Select bmi1 bmi2 ;
 Begin Matrices = Group 1;
  M Full 1 nvar2 Free
 Means M;
 Covariance 
 Option RSiduals;
End

Title G3: 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
 Option RSiduals
End

Title G4: Standardization
 Calculation
 Begin Matrices = Group 1;
 End Matrices;
 Begin Algebra;
  V=A+C+E+D;              ! total variance
  P=A|C|E|D;              ! concatenate parameter estimates
 End Algebra;
!ADE model
! Interval S 1 1 - S 1 4
 Option NDecimals=4
End