require(ribd)
id<-c('gf','gm','b1','b2','m','c')
fid<-c(0,0,'gf','gf',0,'b1')
mid<-c(0,0,'gm','gm',0,'m')
sex<-c(1,2,1,1,2,0)
x<-ped(id,fid,mid,sex)
IBD<-pedigreeIBD(x,c('m','c','b1','b2'))
kappaIBD(x,c('m','c','b1','b2'))
Run the code above in your browser using DataLab