# reg.out <- ibdreg(formula=~1, ped.id=ped.id, person.id=person.id,
# status=status, status.method="AA", data=data,
# ibd.dat=ibd.dat.obj, ibd.var=ibd.var.obj)
# plot.linkage.tests(reg.out$AA.linkage, title="AA, formula: ~1")
Run the code above in your browser using DataLab