If your datasets are small, then it is not difficult to "roll your own" using optim(). If you look at qimr.edu.au/davidD/sib-pair.R you will find such a routine at lines 1511-1562. This calls kinship.rel() that produces NRM etc. -- | David Duffy (MBBS PhD) ,-_|\ | email: davidD at qimr.edu.au ph: INT+61+7+3362-0217 fax: -0101 / * | Epidemiology Unit, Queensland Institute of Medical Research \_,-._/ | 300 Herston Rd, Brisbane, Queensland 4029, Australia GPG 4D0B994A v