!
! One relative pair per family so can test using Sham & Curtis' ERPA
! However 9 and 10 are inbred so these are problematic
!
! Fam  E     O
!  1  0.5   0.897045
!  2  0.5   0.25
!  3  0.25  0.32925
!  4  0.25  0.0
!  5  0.125 0.245435
!  6  0.125 0.145745
!  7  0.25  0.31399
!  8  0.25  0.281245
! 10  0.5   0.731725
!
set workdir d:\
set iter 10000
set locus trait aff
set locus marker mar
read pedigree pairs.ped
run
freq
set plevel 1
apm trait
! apm trait ibd
! write linkage pairs.pre
! write locus linkage pairs.loc
! 