!
! Mooser et al  Am J Hum Genet 61: 402-417.  
! Contents of Table 1:
! Plasma Apoa genotypes, Apoa isoforms, and Lpa levels, 
! in 52 African-American families.
!
! Genehunter (run on nuclear families):
! lod=9.06 m_m=36.85 m_f=38.96 VA=61.36 VQ=862.38 VE=209.49
! (AE model: VA=1122.65, VE=102.21)
!
set imp -1
set iter 1000
set impute on 
set data ../pedigrees
set loc age qua
set loc Lpa qua
set loc Apoa mar
read ped apoa.ped
run
set loc llpa qua
llpa=log(Lpa+1)
nuc
set ple 1
des llpa
regress Lpa on Apoa
qtl llpa
sib llpa
vis llpa
he1 llpa
bye
des Lpa
regress Lpa on age 
regress Lpa on age Apoa
hwe
sib Lpa sim
he1 Lpa sim
keep Lpa Apoa
write locus gh apoa.loc dummy
nuclear 
write gh apoa.pre dummy
bye
all Lpa
asp Lpa over 40
asp Lpa under 40
bye
desc Lpa
assoc Lpa
breed Lpa 0.95
sib Lpa
keep Lpa
write ped x.ped
bye
all Lpa ano
freq
tdt Lpa over 40
asp Lpa over 40
asp Lpa under 40
set iter 5000
apm Lpa over 40
apm Lpa under 40
sib Lpa
trans Lpa 1 -5 0
sib Lpa

