!
! Wellcome IDDM genome scan chromosome 2 data Nature Genet 1994
! 96 families,  384 subjects
! markers: d2s131 d2s165 d2s177 d2s134 d2s139 d2s160 d2s114 d2s142 d2s326  
!          d2s152 d2s72 d2s157 d2s126 d2s206 d2s125
!
set data ../pedigrees
set ple -1
set imp off
set loc iddm aff
set loc D2S131 mar 18.856
set loc D2S165 mar 43.081
set loc D2S177 mar 52.663
set loc D2S134 mar 66.750
set loc D2S139 mar 94.069
set loc D2S160 mar 107.808
set loc D2S114 mar 143.788
set loc D2S142 mar 154.909
set loc D2S326 mar 172.976
set loc D2S152 mar 194.262
set loc D2S72  mar 211.695
set loc D2S157 mar 222.804
set loc D2S126 mar 236.259
set loc D2S206 mar 248.791
set loc D2S125 mar 253.992
read ped iddm.ped
run
lin
ckibs
bye
write loc tcl iddm.tcl
write tcl iddm.asp
bye
tdt iddm
set tdt first
keep iddm D2S134
tdt iddm
bye
ckibs
set iter 100
set plevel 3
share
keep $m
write loc lin iddm.loc
write lin iddm.pre
% makeped iddm.pre iddm.ppd n 
% echo 1                  > relatwin.in
% echo output.sib        >> relatwin.in
% echo iddm.loc iddm.ppd >> relatwin.in
% echo                   >> relatwin.in
% echo 0.90 0.08 0.02    >> relatwin.in
% echo 0.94 0.06         >> relatwin.in
% echo 0.99 0.01         >> relatwin.in
% /home/relative/relative relatwin.in
! % rm relatwin.in iddm.loc iddm.pre iddm.ppd

! write loc tcl iddm.tcl
! write tcl iddm.asp
! % sib_phase -f iddm.tcl iddm.asp
! % rm iddm.tcl iddm.asp
! asp iddm
! tdt iddm
