! E4. Genetic Cholesky Model
! Schieken Skinfold Measurements of Male Twins

G1: genetic structure
Data Calc NGroups=4
Matrices
A Lower 4 4 Free  ! genetic cholesky (lower triangular)
Compute A*A' /
End

G2: specific environmental structure
Data Calc 
Matrices
E Lower 4 4 Free  ! specific environmental cholesky
Compute E*E' /
End

G3: MZ twin pairs
Data NInput_vars=8 NObservations=84
Labels
Bicep1 Subsca1 Supra1 Tricep1 Bicep2 Subsca2 Supra2 Tricep2
CMatrix
.1285 
.1270  .1759 
.1704  .2156  .3031 
.1035  .1101  .1469  .1041 
.0982  .1069  .1491  .0824  .1233 
.0999  .1411  .1848  .0880  .1295  .1894
.1256  .1654  .2417  .1095  .1616  .2185  .2842 
.0836  .0907  .1341  .0836  .1010  .1134  .1436  .1068
Matrices
A Symm 4 4 = %E1  ! expected matrix of group 1
E Symm 4 4 = %E2  ! expected matrix of group 2
Covariances A+E | A _
            A   | A+E /
End 

G4: DZ twin pairs
Data NInput_vars=8 NObservations=33
Labels
Bicep1 Subsca1 Supra1 Tricep1 Bicep2 Subsca2 Supra2 Tricep2
CMatrix
.1538 
.1999  .3007 
.2266  .3298  .3795 
.1285  .1739  .2007  .1271 
.0435  .0336  .0354  .0376  .1782 
.0646  .0817  .0741  .0543  .2095  .3081 
.0812  .0901  .0972  .0666  .2334  .3241  .3899 
.0431  .0388  .0376  .0373  .1437  .1842  .2108  .1415
Matrices
A Symm 4 4 = %E1
E Symm 4 4 = %E2
H Full 1 1
Covariances A+E   | H@A _
            H@A   | A+E /
Matrix H .5
Start .5 All
Options Multiple NDecimals=4
End