data(hapRes)
res <- hapRes$res
sPF <- hapRes$sPF
annot <- hapRes$annot
nnL <- length(Z(res)[1,])
labelsA <- cbind(as.character(1:nnL),
as.character(1:nnL),as.character(1:nnL),
as.character(1:nnL))
resIBDsegmentList <- extractIBDsegments(res=res,
sPF=sPF,annot=annot,chrom="1",labelsA=labelsA,
ps=0.9,psZ=0.8,inteA=50,thresA=6,mintagSNVs=6,
off=0,procMinIndivids=0.1,thresPrune=1e-3)
summary(resIBDsegmentList)
print("Position of the first IBD segment:")
print(IBDsegmentPos(resIBDsegmentList[[1]]))
print("Length of the first IBD segment:")
print(IBDsegmentLength(resIBDsegmentList[[1]]))
Run the code above in your browser using DataLab