!
! Pedigrees from Keavney et al
! Measured haplotype analysis of the angiotensin-I converting enzyme gene
! Hum Mol Genet 1998 7: 1745-1751
!
! as included in MENDEL 4 example pedigree files: ACE I/D polymorphism and
! sex-standardised sACE.
! 
set data ../pedigrees
set loc aceins mar
set loc sace qua
read ped ace.ped
run
des
hist sace

set ple 1
var sace
qtl sace
ass sace
tdt sace over 0
set ple 0
!
! Does the I/D polymorphism association explain all the linkage signal?
!
set loc adjace qua
set loc g11 aff
set loc g12 aff
set loc g22 aff
set loc lev qua
factor aceins lev
tab aceins lev
g11=(lev==1)
g12=(lev==2)
g22=(lev==4)
tab aceins g11 g12 g22
adjace=sace
residuals adjace on g12 g22
set ple 1
des adjace
var adjace
qtl adjace
